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AB ST RACT 


As part of a continuing program of Quaternary stratigraphic studies 
in the Edmonton area, Alberta, till fabric determinations were made in 
an attempt to distinguish the till units. As the local stratigraphy be- 
came established, emphasis of the till fabric work was changed to con- 
centrate on (1) examination of possible models for the statistical 
treatment of till fabrics, (2) investigation of the problems of deter- 
mining the pebble fabric at an exposure and of interpreting the fabric 
pattern, (3) evaluation of till fabrics as an aid to local stratigraphic 
studies, and (4) establishment of local ice-movement directions. 

Two mathematical models were examined to determine their value as 
bases for calculating descriptive statistics to represent till fabric 
data. A previously used model, that of a single axially symmetric 
spherical normal (Fisher) distribution, was examined and found unsuitable 
since most of the fabrics encountered were not unimodal. The testing of 
a model consisting of several Fisher distributions, one corresponding to 
each mode of the fabric, revealed that most of the fabric modes lacked 
axial symmetry. 

The orientation of 550 elongate pebbles from a 300-foot exposure 
of a single till unit (six feet thick) revealed that: (1) The fabric 
varied both laterally and vertically, but when samples were taken to- 
gether they showed a preferred trend parallel to groove molds at the 
base of the till. (2) Larger pebbles exhibited less scatter in their 
orientations than did smaller ones. (3) The orientation of the a-axis 
of a pebble was in part determined by the elongation of the pebble. 


Pebble fabrics appear to be of l'ittle help in the recognition and 
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differentiation of Quaternary stratigraphic units in the Edmonton area. 

Upper till fabrics are consistent with ice movement from the north- 
east as indicated by sole markings and surface features. Transverse 
fabrics were found at four out of 10 locations. 

The direction of movement of the glacier that deposited the lower 
till cannot be reliably determined from the fabric evidence presented 
here. Interpretation of lower till fabrics is made difficult by the fact 
that deformation of this deposit occurred when the area was overridden 
by the glacier that laid down the upper till. Distinction between pri- 
mary and wholly or partially altered lower till fabrics is a major 


problem in the study of lower till fabrics. 
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INTRODUCTION 


Historical Review 


Till Fabrics 


On March 8, 1856 Henry Youle Hind presented a paper to the 
Canadian Institute entitled '"'0n the blue clay of Toronto" in which 
he described the orientation of the stones and suggested a glacial 
origin for the deposit rather than deposition from floating ice or 
by the sea (Elson, 1966). Miller (1884) in describing ''pavement 
boulders'' in the till near Edinburgh, Scotland, stated (p.167): ''The 
longer axis of the stone is often directed in the line of glaciation, 
and the pointed end is frequently, but not always, toward the ice.'' 

He also suggested that the orientation of the Stones might reflect 
their original orientation in the ice. Bell (1888) observed a tendency 
for boulders on Swiss glaciers to assume a longitudinal position. 
"Upham (1891) described the characteristic position of 'oblong' stones 
in sub-glacial till as having their long axes parallel to contiguous 
striae. Flat stones were said to lie parallel with the surface of 
deposition.'' (Holmes, 1941, p.1302). 

The first systematic investigation of till fabric appears to 
have been carried out by Richter in the early 1930's. According to 
Holmes (1941, p.1302) Richter's conclusions confirmed the earlier 
observations of Miller and Bell and were supported by quantitative 
statistical data. Richter was apparently the first to note a tendency 


to an orientation transverse to glacier flow. 
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Holmes (1941) investigated the relationship between the orienta- 
tions of pebbles in a till and their physical properties. Using a total 
of 1180 pebbles from till of central New York he attempted to demonstrate 
that groups of pebbles having different combinations of physical proper- 
ties (shape, rounding, axial ratios, size) had significantly different 
preferred orientations. His work drew attention to many important 
questions, and led him to suggest the ''plastering-on'' process of till 
deposition. 

In an attempt to answer some of these questions, Harftison (1957) 
studied the orientation of disk- and blade-shaped pebbles having b:c 
ratios greater than two. He felt that "particles of this shape are 
better indicators than any others of dynamic forces responsible for 
their orientation." (Harrison, 1957, p»2/6)4 His work showed that disk- 
and blade-shaped particles tend to lie in a near-horizontal plane and 
to be imbricated upstream to ice movement. Harrison favoured the theory 
that the orientations of the pebbles in a till reflect the orientations 
that they achieved in the ice. 

Glen, Donner and West (1957) showed that the observed orientations 
of till pebbles could be entirely the result of flow in moving ice. 
Collisions between particles and the. mode of deposition could change 
the distribution, however. 

More recently, attention has been given to till fabric variability, 
the reliability of fabric determinations, and problems of sampling. For 
example, Kauranne (1960), Andrews and Smith. (1966), Andrews and King. - 
(1968), Johansson (1968) and Harris (1969) have investigated till fabric 
variability. These studies have shown that till fabrics may be highly. 


variable, in both preferred orientation and degree of preferred orientas 
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tion. Thus any till exposure must be regarded as a unique problem in 
sampling. 

The problem of determining sediment origin from pebble fabric 
patterns has apparently received little attention. Lindsay (1968), 
in relation to simulation of mudflow pebble fabrics by computer, 
briefly considered the problem of differentiating between till pebble 
fabrics and mudflow pebble fabrics, but the work was inconclusive. 
Recently, with more promising results, he compared the pebble fabrics 
of two diamicton deposits with known till fabrics to test the theory 
that the deposits were tillites (Lindsay, et al., 1970). 

Previous work on the Quaternary sediments of the Edmonton area 
(fig. 1) has been summarized by Westgate (1969). Although the basic 
Stratigraphy has been established, little work has been done on the 
till fabrics of the area. Till fabric data obtained by other workers 


in recent years has been incorporated in this report. 


Statistical Treatment of Geological Orientation Data 


Much work has been done on the statistical treatment of orienta- 
tions in two dimensions. One important method involving vector sum- 
mation was apparently first used by Reiche (1938) and was applied to 
the orientations of pebbles in sediments, including till, by Krumbein 
(1939). A second powerful method developed by Tukey (1954) using a 
combination of the chi-square technique and vector summation was 
applied to the orientations of till pebbles by Harrison (1957). A 
simple chi-square test was used on till fabric data by Kauranne (1960) 
and forms the basis of a method described by Harris (1969) for deter- 


mining minimum satisfactory sample sizes. 
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The present study deals with methods of treating the orientations 
of pebble axes in three dimensions. For the treatment of directions in 
three-dimensional space, spherical probability distributions were 
constructed by Arnold (1941) but were not published until included by 
Pincus in his comprehensive review of methods of analyzing orientation 
data in 1953. In the same year Fisher (1953) developed the theory of 
the spherical normal model as applied to the axially symmetric case 
for the treatment of geomagnetic vectors. Further statistical pro- 
cedures based on Fisher's work were developed by Watson (1956a, 1956b, 
1960), Watson and Williams (1956), and Watson and Irving (1957). A 
somewhat different model was suggested by Scheidegger (1965) designed 
for the treatment of non-directed axes. 

The procedure described by Watson and Irving (1957) for the 
treatment of paleomagnetic vectors was adopted by Steinmetz (1962) 
for representing measurements of direction and angle of dip of cross- 


beds, and by Andrews and Shimizu (1966) for till fabric data. 


Objectives of Present Study 


The initial objective of this study was to map exposures of 
glacial deposits in the Edmonton area so as to get a better under- 
standing of the stratigraphy. Till fabric determinations were made 
in an attempt to distinguish the till units. As the local stratigraphy 
became established, emphasis of the work was changed to concentrate on 
(1) evaluating the reliability of the till fabric measurements, (2) 
assessing the usefulness of till fabrics in the local stratigraphic 


studies, and (3) establishing local ice movement directions. 
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Later a fourth major objective was added: the application and 
evaluation of Andrews and Shimizu's (1966) method of calculating descrip- 
tive statistics to represent till fabric samples and serve as a basis 


for comparison of samples. 


The Approach 


Since the principal objective initially was to map as many new 
sections as possible, only one fabric measurement was made on each till 
unit at each site, except in the case of certain important exposures 
where it was considered worthwhile obtaining more detailed fabric infor- 
mation. 

The reliability of the till fabric determinations was investi- 
gated by taking duplicate samples from certain sites to test the 
adequacy of 50-pebble samples, and by measuring both lateral and 
vertical sequences of fabrics to test the within-site variability and 
the value of single samples. At one location, the orientations and 
dimensions of about 200 pebbles were measured to see if any effects of 
pebble shape or size could be detected in the fabric patterns. 

Based on a FORTRAN IV program developed at the University of 
Alberta for producing stereographic projections of point density, a 
modified program was developed that would take account of measurement 
errors and calculate the probable density of true axes rather than the 
actual density of recorded axes. 

Andrews and Shimizu's (1966) adaptation of Watson and Irving's 
(1957) procedure for calculating descriptive statistics for three- 


dimensional orientations was employed. Statistics calculated by the 
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method were compared with point density diagrams of the data in order 
to evaluate the success of the method. An alternative method devised 
by the writer based on a different assumption about the distribution 


of pebble orientations was also used and evaluated. 
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PROCEDURES OF DATA COLLECTION 


AND GRAPHICAL REPRESENTATION 


Data Collection 


Definitions 


The dimensions and shape of a pebble are described by the lengths 
of three imaginary mutually perpendicular ''axes''. The c-axis is per- 
pendicular to the maximum projection plane (Krumbein, 1939, p.677) of 
the pebble. The a-axis lies in this plane, and is parallel to the 
longest dimension of the pebble in this plane. The b-axis is perpen- 
dicular to both the a- and c-axes. The a-, b-, and c-axes are also 
referred to as the long, intermediate, and short axes respectively. 
The planes perpendicular to these axes are referred to as the bc-plane, 
the ac-plane, and the ab-plane. 

The lengths of the three axes are indicated by the lower case 
letters a, b, and c. Pebble shape is often described by the ratios 
of these lengths, for example, the length of the a-axis divided by 
the length of the b-axis is called the a:b ratio or simply a/b. The 
a:b ratio is here called the ''elongation'' of a pebble, and the b:c 
ratio the "flatness". 

The definitions of the a-, b-, and c-axes used here are the same 
as those described by Krumbein (1939). However, it should be noted 
that Krumbein's statement that ''the a-axis is normal to the minimum 


projection plane! (Krumbein, 1939, p.678) is not necessarily true. 
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Selection and Measurement of Pebbles 


At most sites only the orientation of the a-axis of each pebble 
was measured. At location F, the lengths a, b, and c were recorded 
for some of the pebbles, and for some of these the attitude of the 
ab-plane was recorded in addition to the orientation of the a-axis. 

Particles smaller than about one quarter inch were not measured 
due to the difficulty of obtaining accurate measurements. Pebbles 
smaller than about one inch were measured only when their use was 
necessitated by the scarcity of larger pebbles. Pebbles with an a:b 
ratio of less than about 1.2 were not measured, except at location F 
as described in Chapter 6. 

Most samples consisted of measurements of 50 pebbles. The 
volume of till from which a sample was taken varied from about one 
cubic foot to about six cubic feet, depending on the abundance of 
pebbles and the condition of the till. 

The number of samples taken from one till or till-like unit at 
a locality varied, the minimum being one and the maximum 11. When 
two or more samples were taken from a unit at a single location, they 
were located so as to yield the greatest amount of information about 
both vertical and lateral fabric variability within the limits of the 


exposure. 


Field Technique and Sources of Error 


Care was taken to select sample sites in undisturbed till, 


avoiding slumped material. After a sample site had been selected, 
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the existing surface of the outcrop was removed to a depth of about 
six inches, or as much as was necessary to expose fresh, undisturbed, 
workable material. 
The procedure of exposing and measuring pebbles was as follows. 
The till was removed using a small pointed trowel until a pebble was 
encountered. The pebble was held firmly in position while the sur- 
rounding matrix was carefully removed until about half of the pebble 
was exposed. The pebble was then removed, while care was taken to 
preserve its mold in the remaining matrix. At this point the pebble 
was discarded if it had too small an elongation to be usable. If the 
pebble was usable, the position of its long axis was estimated, but 
not marked on the pebble. The pebble was then carefully replaced in 
its mold and the trend and plunge of its long axis were measured to 
the nearest degree using a Brunton compass. In some early samples 
the trend and plunge were measured to the nearest five degrees. The 
sources of error in the measurements thus obtained are thought to be 
as follows. 
1. Error involved in replacing pebble in its original position. 
This error was probably less than two degrees in most cases. 
2. Error involved in estimating the position of the long axis 
after replacement of the pebble. This error was highly 
variable, depending on the shape of the pebble, how much of 
the pebble was exposed, which portion of the pebble was ex- 
posed, and the attitude of the pebble in relation to the 
outcrop surface (that is the angle from which the pebble had 
to be viewed by the observer). The maximum probable value 


of this error ranged up to about five degrees. 
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3. Error in aligning Brunton compass with estimated position 
of long axis. This error was probably less than two de- 
grees in all cases. 

4, Error in reading scale: one half degree. 

Thus the discrepancy between the recorded orientation of an 

axis and the true orientation of the axis was considerably less than 


10 degrees in all cases. 


Graphical Representation: 


Contoured Point Density Diagram 


Introduction 


Initially, the writer obtained equal-area projections of point 
density using a FORTRAN IV program developed at the University of 
Alberta (Cruden, 1966, p.84). This program determines point density 
by counting axes or the poles to planes directly on the reference 
hemisphere. The axis of a circular cone whose apex is at the centre 
of the hemisphere is placed successively through each of 333 counting 
locations, and the number of observations falling within the cone at 
each location is recorded as a percentage of the total. The solid 
angle contained within the cone can be set to any desired fraction of 
the total solid angle of the hemisphere. This cone defines a circle 
on the surface of the hemisphere whose area is the same fraction of 
the total surface area of the hemisphere. The program uses the 
computer's line printer to produce a 10-inch radius equal-area pro- 


jection of the counting locations, each location being represented by 
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the appropriate point density value. To give a simple line-printer 
format, the counting locations were chosen so that their projections 
form a rectangular grid with a spacing of one tenth the radius of the 
reference hemisphere. Sixteen locations near the periphery of the 
projection were added to make a total of 333. 

Suppose that a sample of N axes is taken from a population of 
axes, and that the program is used to prepare a density diagram using 
a P per cent counting cone (or circle). The contoured diagram then 
represents a ''density surface", the elevation of which at any point 
shows the percentage of axes in the sample falling within a P per cent 
counting circle centered at that point, and gives the probability of 
an axis selected randomly from the sample falling within a P per cent 
counting circle centered at that point. It is taken as an estimate of 
the probability that an axis selected randomly from the population will 
fall within a P per cent counting circle centered at that point. 

This method of determining point densities is a great improvement 
over previous methods of counting points on a projection. Previously, 
accurate contouring on an equal area projection required an elliptical 
counter whose eccentricity increased progressively with distance from 
the centre of the projection, so that use of a circular counter intro- 
duced inaccuracy that was greatest at the periphery and decreased to 
zero at the centre. Since the present method counts the points directly 
on the reference hemisphere using a circular cone, all such inaccuracy 
is eliminated. Use of a computer also eliminates plotting and counting 
errors and allows diagrams to be prepared in a small fraction of the 


time required to prepare them manually. 
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Modified Procedure 


In the program just described, a counter is incremented by one 
for each observation that falls within the counting circle, regardless 
of where in the counting circle it falls. Thus an observation falling 
just inside the counting circle, for example at point A in Figure 2, 
is counted as one, while an observation falling just outside the counting 
circle, for example at point B in Figure 2 is not counted at all. 

With relatively small samples, this results in erratic contours. 
Also, in interpreting the diagram in terms of the distribution of axes 
in the population, it is assumed that the axes in the sample are identi- 
cal to the axes that they represent in the population. In other words, 
errors of measurement are ignored, as are inherent variations in the 
attitude of the feature being measured. For example, when the attitude 
of a joint that is not perfectly planar is represented éy a single 
measurement, that measurement is assumed to indicate the mean orienta- 
tion of the joint: both measurement errors and the variation in attitude 
or "roughness" of the joint are ignored. It was felt worthwhile to 
develop a modification of the program that would take into account the 
distribution of errors in the measurements, thereby smoothing the contours 
and producing a diagram that is more likely to approximate the actual 


distribution of axis density in the population. 
Method Used in the Program 


Following a basic idea originally used by Dr. D. M. Cruden 


(Cruden, 1966), the writer developed a program that would take measurement 
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Counting Location 


Counting Circle 


Surface of 


Reference Sphere 


Figure 2. Section through centre of reference sphere, two sample axes 
A and B, and one counting axis, illustrating the difference 
between the methods of incrementing the counter in the mod- 
ified and unmodified programs for preparing point density 
diagrams. 
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errors into account. The problem was to design a counting procedure 
that, rather than counting 1 for each observed axis that fell within 
the counting circle, would increment the counter by a quantity Pr, 
where Pr is the probability that the true axis represented by the 
observation lies within the counting circle. Errors in the measurement 
of an axis were assumed to conform to a spherical normal or Fisher 
distribution (Fisher, 1953). Each observation was then no longer 
thought of as a point on the surface of the reference hemisphere, but 
rather as a spherical normal probability density surface centred at the 
observed point, as illustrated in Figure 2. The shape of the Fisher 
distribution is described by a parameter K, called the precision para- 
meter. K must be estimated by repeated measurements of selected pebbles. 
After a value for K has been established, it is assumed to be constant 
for all measurements. 

If the probability of occurrence of a point in a two-dimensional 
sample space is represented by a probability density surface, then the 
probability of that point occurring within a specified region of the 
Sample space is given by the volume beneath the probability surface 
over the specified region. For the modified point density program, 
the two-dimensional sample space is the surface of the reference hemi- 
sphere; the point is the true axis represented by an observation, and 
its probable position is described by a spherical normal probability 
density surface centred at the observed position; the specified region 
of the sample space is the counting circle. Hence the probability of 
the point (true axis) falling within the counting circle is given by 
the volume beneath that part of the probability surface that lies 


within the counting circle. 
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To facilitate computation of this volume, the program approximates 
the probability density function used to represent each true axis by a 
series of concentric discs arranged in a vertical stack, with the radii 
of the discs decreasing from bottom to top (fig. 2). Progressively 
larger discs are added to the bottom of the stack until the approximation 
contains more than 99.9 per cent of the volume beneath the probability 
surface, that is, more than 99.9 per cent of the probability. 

For a given observation and a given position of the counting 
circle the required volume is determined by calculating how much of 
the area of each disc is contained within the counting circle and 
multiplying this area by the height of the disc, then adding the 
results for all the discs that overlap the counting circle. The 
result is a number between 0 and 1, representing the fraction of 
the total volume of the distribution that overlaps the counting circle. 
For a given position of the counting circle, this fraction is computed 
for each observation, and the fractions for all observations are added. 
The result is expressed as a percentage of the total number of observa- 
tions, and recorded as the density value for that position of the 


counting circle. 


Statistical Significance of Modified Procedure 


It will now be shown that the counting procedure just described 
determines, for a specified position of the counting circle, the 
probability that the true axis represented by an observation selected 


randomly from a sample of N observations lies within the counting circle. 
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The selection from N observations of an observation whose true 
axis lies within the counting circle can happen in N different ways: 
(1) selection of the first observation, and its true axis lying within 
the counting circle, (2) selection of the second observation, and its 
true axis lying with the counting circle, and so on up to: (N) selection 
of the Nth observation, and its true axis lying within the counting 
circle. The occurrence of any one of these N events itself consists of 
the simultaneous occurrence of two events: (1) the selection of the 
particular observation, and (2) its true axis lying within the counting 
circle. The probability of the simultaneous occurrence of two events 
is given by the product of the probabilities of the individual events 
occurring alone. Hence the probability of obtaining a true axis within 
the counting circle in any one of the above N ways is (1/N) x Pres 
where 1/N is the probability that the ith observation will be selected, 
and Pr. is the probability that its true axis lies within the counting 
circle. Furthermore, the N different ways of obtaining an axis within 
the counting circle are mutually exclusive, and the probability of the 
occurrence of at least one of a number of mutually exclusive events 
is given by the sum of the probabilities of the individual events. 
Hence the probability of obtaining an axis within the counting circle 
in at least one of the N mutually exclusive ways is S x Fe, SF" Sinee 


L=| 
, oS 1 : 
N is constant, this is the same as —2 Pr; . But 2 PY is the 
© re 


N 


quantity determined for each counting circle by the counting procedure 
described above, and hence when this is divided by N the result is the 
probability that an observation selected randomly from the N observations 


represents a true axis lying within the counting circle. 
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Determination of K 


The meaning of K in terms of the tightness of grouping of the 
population it describes is illustrated by the following table. Each 
value of K is given with its corresponding value of bo5? the significance 
of which is as follows. If an axis selected randomly from the population 
makes an angle ¢g with the central direction of the distribution, then Gor 


is an angle such that Pr(¢< or )=0.95. 


K 95 (degrees) 
10 Lo .6 
20 31.9 
50 19.9 
100 ba] 
200 9.9 
500 623 
1000 4.5 


So that the program can approximate the probability distribution 
representing the true axis of a pebble, the degree of scattering of 
repeated measurements, must be specified by means of the precision 
parameter K. Two sets of repeated measurements on single pebbles were 
made by the writer to test the assumption of a spherical normal distribu- 
tion and to estimate K. These sets of measurements are labeled 991 and 
992, and consist of 23 and 20 measurements respectively. For set 991, k, 
the best estimate of K, is 695. For set 992, k is 119. A value of 
K=500 was more or less arbitrarily chosen and used to prepare modified 
point density diagrams. The value of 6.3 degrees for Boo implied by 
this value of K appears to be reasonable in view of the probable 


magnitude of errors discussed above. 
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Validity of Assumptions 


It is thought likely that K varies from pebble to pebble, 
depending chiefly on the shape of the pebble. Other factors mentioned 
in the section on errors in the measurements may also have an effect. 
The fit to Fisher distributions with the estimated parameters was 
tested using the modified mean vector program described in Chapter 3 
and given in Appendix D. Both sets 991 and 992 contained too few 
observations to allow a chi-square test of the distribution of the 
angles between the observations and the estimated mean direction. A 
chi-square value was calculated for each axial distribution about the 
mean, however. The chi-square value was 5.00 with one degree of free- 
dom for set 991, indicating that there is a probability of less than 
five per cent that this set has a Fisher distribution. For set 992, 
the chi-square value was 13.60 with one degree of freedom, indicating 
a probability of considerably less than one half of one per cent that 
this set has a Fisher distribution. Hence the assumption that the 
measurement errors have a Fisher distribution is at best an approxima- 
tion. 

It is evident from the foregoing discussion that the two major 
assumptions involved in the modified procedure are both invalid when 
the technique is used on field measurements of long axes of till stones. 
However, reasons will be advanced for using the modified diagrams 


in spite of the invalidity of these assumptions. 
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Reasons for Using the Modified Program 


Modified diagrams were prepared using a three per cent area, and 
were compared with diagrams obtained using the unmodified point density 
program using a three per cent area. Figure 3 shows the three per cent 
diagrams produced by both programs for samples 1, 2, 6 and 34-1. For 
all four samples, the modified program produces smoother, more regular 
contours. In some cases the unmodified program produces areas contain- 
ing numberous small highs and lows, whereas the modified program simpli- 
fies the contours greatly, eliminating many of the small highs and lows. 
This can be seen, for example, in sample 34-1 in the southwesterly and 
southeasterly parts of the diagrams, and in the western part of the 
diagram for sample 6. If it is assumed that the actual distribution 
of axes in the population has a simple pattern, then the simplified 
contours must be regarded as a better approximation to the actual 
distribution of density in the population. 

In certain situations the modified program eliminated maxima 
produced by the unmodified program that could be shown to be spurious. 
These may appear when the counting circle is located between two group- 
ings of observations so that it partially overlaps both. One effect of 
the modified counting procedure is to reduce the effect of points lying 
near the periphery of the counting circle, hence largely avoiding the 
above difficulty. This can be seen in sample 1, at the periphery of 
the projection in the northeastern and southwestern quadrants. 

In view of the above-mentioned advantages of the modified program, 
it was decided to use it to illustrate the fabrics rather than the un- 


modified program. It must be remembered, however, that in view of the 
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invalidity of the assumptions concerning error distribution in the 
measurements, the numbers appearing on the diagrams cannot be taken as 
accurate estimates of probability. Nevertheless, the numbers probably 
do represent useful approximations of the probabilities involved. 
Similarity of the unmodified and the modified diagrams in their major 
features and relative magnitudes in different areas of the diagram 
bears this out, and shows that no loss of information is incurred by 
use of the modified diagrams. Furthermore, comparison of the maximum 
value appearing on diagrams for different samples using the same size 
counting circle is believed to give a realistic indication of the rela- 
tive strength of the modes, that is, of the relative strength of pre- 
ferred orientation in the data, taking into account a certain amount of 
expected random variation in the samples. The question of interpreta- 


tion of the diagrams is discussed more fully in Chapter 4. 


Selection of Size of Counting Circle 


Following experience gained with the unmodified point density 
program, modified diagrams were prepared using a three per cent counting 
circle. Later, modified diagrams were also obtained using a seven per 
cent circle in order to achieve more of an ''averaging'' effect. It was 
felt this would be helpful in the interpretation of the actual direction 
of preferred orientation represented by each fabric sample. 

Figure 4 shows modified diagrams prepared for samples 2, 43, 44 
and 45 using counting circles of one, three, five, seven and nine per 
cent. All diagrams were contoured using a contour interval of e, where 


e is the "'expected'' density in the case of a random sample from a 


rot em +tett nt = betttbom etl? bre bot? thom: ant to en 
 mexestb eft te sor Inevettib ni eebutinpem evizeiey bos 2zowwiss? = 
yd ber suont et not Temorat eo. 2eor on ted? ewore bis eto eid vi cael y | 
munitxem ond Yo Host reqmos .eromiediiuy .amerpsth beittbom edt 1 ‘ 
Ste emee oft pereweetipiee” fne19t3 th 10% zmeripeib no. oneal on 
-efor ont to rio toes tbat ote ise & ovie es bevet ied ef etovis a 
| of 


4 


Basie, to dagnonte aviteter att Yo yet enh <2ebom srit to dipnevte s 

to inudme ovetiso 6 Swosss orn? BREA BSeb sit rit notsesnetw 
-sisigqrsiar to elena edt .eglqmee SA nt nottetisy mobney eee) 
a aetqed9 nt ytful rom bezewsetb ei amaipsib odd remot 


<i 


“etartd pnttnucd to ssi2 to cena 
ry i 
. i ‘a 
ytfaneb tniog bettibomy sit dtiw benisp sansitteqxe priwol lo? - Ad”. ae 
pntsnues lin Nase fs B pnizy borsge7q svow emerges tb bottibom a i 
oe 
oq sores a Chheu tebaddy ozs sei eieiperb 'bettibom etei aia ae 
re 
@ow. 31 . toot te i leaeacll radi Siem svstros ot eatad nt efovig 
nor sei saaeie ort hell 


a" 


22 


uniform population. For example, for a diagram based on a three per 

cent counting circle, e is three per cent; when a seven per cent counting 
circle is used, e is seven per cent. The contours thus provide an 
indication of the degree of deviation from uniformity in the sample. and 
facilitate comparison of diagrams prepared using counting circles of 
different size. 

It can be seen that the contours become more and more generalized 
as the counting circle increases in size. The three per cent diagrams 
give a detailed picture of the sample, indicating unmistakably all 
major and even minor groupings. The double peak in sample 43 is clearly 
evident in the three per cent diagram. The seven per cent diagrams do 
not differ markedly from the three per cent diagrams, but smaller highs 
are usually assimilated into the larger ones. The two maxima about 
60 degrees apart in the northwestern and northeastern quadrants of 
sample 2 remain separate in the seven per cent diagram. Even the two 
close maxima in sample 43 remain distinct in the seven per cent plot, 
but merge into a single elongate high in the nine per cent plot. 

The position of a maximum on a point density diagram will be 
indicated by giving the trend and plunge of the line passing through 
the density peak. In sample 2, the position of the largest maximum 
moves from 032 25 to 026 16 as the counting circle increases from three 
to nine per cent. In sample 43, as counting circle size increases from 
one to seven per cent, the north-northwesterly-trending maximum moves 
from 344 32 to 356 29, and the north-northeasterly-trending maximum 
from 008 12 to 019 07. With a counting circle of nine per cent, a 
single maximum appears at 005 20. The trend of the principal maximum 


in sample 44 migrates from 348 in the one per cent plot to 353 in the 
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five per cent plot, then back to 349 in the nine per cent plot. Plunge 
of the principal maximum in sample 45 decreases continuously from 15 
degrees to six degrees in going from the one per cent to the nine per 


cent plots, while the trend changes from 303 to 288 and back to 296. 


Various rules have been used for determining the optimum size of 
the counting circle. For instance, Flinn (1958) used a circle of per 
cent area 100/N where N is the number of observations, so that the 
expected areas between the contours would be the same regardless of 
sample size. 

Kamb (1959), in constructing ice crystal fabric diagrams of the Blue 
Glacier, used a different approach. He observed that, in the case of 
a random sample of N points from a uniform population, the actual 
number of points, n, falling within the counting circle has a binomial 
distribution, with mean NA and standard deviation o=/NA(1-A) where A 
is the area of the counting circle assuming the area of the projection 
to be unity. In order to avoid wildly fluctuating counts as a result of 
random sampling, he sets the arbitrary condition that the mean of n be 
three times its standard deviation. This gives A = 9/(N + 9), or 
A = 900/(N + 9) per cent. Thus the expected density in a random sample 
from a uniform population is n or 36 . Kamb then used contours at 0, 
20, 40, and so on, so that an indication of the degree of deviation 
from uniformity was obtained. Kamb stated (1959, p.1909): 'With the 
choice of A used here, the difference between the Schmidt and Mellis 
contouring methods (Flinn, 1958) can lead to no statistically signi- 
Ficant differences in the positions of the contours and may therefore 


be disregarded."' 
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In the case of a sample of size 50, Flinn's rule gives a counting 
circle size of 100/50 per cent or two per cent. Kamb's rule gives 
900/59 per cent or 15.3 per cent. This difference is not surprising. 
The two workers subjected their diagrams to different analytical pro- 
cedures, and therefore had different requirements of the diagrams. 
Clearly, no universal rule can be devised for determining the ''best'' 
size of the counting circle; in any given study, the use to which the 
diagrams are to be put and the information required from them are the 
ultimate controlling factors in determining the size of counting circle 
to be used. For example, a small counting circle emphasizes the 
small-scale variations in the sample. A larger circle will tend to 
"average out'' these local variations while emphasizing the large-scale 
variations in density within the sample. Thus use of a large counting 
circle is one way of separating the random and non-random elements in 
the sample, assuming that the local density variations in the sample 
are likely to be random while the larger density variations are less 
likely to be random and more likely to reflect characteristics of the 
population. This principle of minimizing the effect of the random 
element is the basis for Kamb's rule for determining counting circle 


size. 
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PROCEDURES OF NUMERICAL REPRESENTATION 


AND ANALYSIS 


Review of Methods for Three-Dimensional 


Orientation Data 


Introduction 

Flinn (1958, p.527) pointed out the great variety of models that 
must be considered in analyzing three-dimensional fabric data: "In 
rock fabrics a single preferred direction is a special case -- more 
commonly several preferred directions (groups of concentrations) occur; 
or alternatively a preferred plane or planes (girdles); or concentra- 
tions in which the direction is more dispersed in one plane than in 
that normal to it (partial girdles and elongate concentrations).'' 

Curray (1956, p.118), in discussing the analysis of two-dimensional 
orientation data, listed the requirements of a method of analyzing 
orientation data as: (1) A measure of central tendency or preferred 
orientation which is independent of origin. (2) A measure of disper- 
sion which is independent of origin. (3) A test of the statistical 
significance of the results against a model of randomness. (4) A 
model distribution from which deviations can be tested. Although Curray 
was concerned only with two-dimensional orientations, the above list 
is equally applicable to the study of three-dimensional orientation data. 

Various statistical procedures have been devised to assist in 
the performance of one or more of the above steps. Two excellent review 
articles dealing with many of these procedures are those by Flinn (1958) 


and Steinmetz (1962). Many of the concepts discussed by Curray (1956) 
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in relation to two-dimensional orientation data are also relevant. 

A brief review of methods that have been used in the treatment 
of three-dimensional orientation data is given in the following sections. 
A subdivision of the methods into the following two groups was found 
convenient. (1) Direct methods: those that treat the data directly 
in a three-dimensional framework. (2) Indirect methods: those that 
treat two-dimensional projections of the data. This subdivision is 


similar to those employed by Flinn (1958) and Steinmetz (1962). 
Directed and Non-directed Orientations 


This section will define some important concepts and terms to 
facilitate subsequent discussion of linear quantities. To begin with, 
any straight line has an orientation in three-dimensional space, usually 
referred to a conventional system of three mutual ly perpendicular axes. 
An orientation may be either directed or non-directed. The orientation 
of the a-axis of a till pebble is an example of a non-directed orienta- 


tion. A directed orientation has a sense. Consider, for example, the 


wind direction indicated by a simple weather vane. This is a directed 
orientation. In this report the word direction will always denote a 
directed orientation. In some cases the term ''directed orientation" 
will be used for emphasis instead of ''direction''. It should be noted 
that a non-directed orientation may be made into a directed orientation 
by assigning a sense to it, either arbitrarily or according to some 
convention. 

A vector is a quantity that has an orientation, a sense and a magni- 


tude. If a magnitude is added to the wind direction mentioned above, 
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the result is the wind velocity, a vector. 

Table 1 gives examples of geological quantities having different 
combinations of the above properties, in both two and three dimensions. 
Observations of types 2 and 3, that is directed orientations, require 
methods of treatment that take sense into account, while such methods 
are not suitable for treating observations of type 1. Similarly, 
methods suitable for treating non-directed orientations (type 1) are 
not suitable for the treatment of directed orientations (types 2 and 3). 
Thus the chief distinction to be made for the purpose of determining 
suitable methods of treatment is between directed and non-directed 


quantities. 


Indirect Methods 


Introduction 


Steinmetz (1962, p.805) argued that the tests that analyze 
three-dimensional data indirectly by means of projections are weaker 
than tests that analyze the data directly. This is probably true if 
a suitable direct method is available. However, several years earlier 
Flinn (1958) had pointed out that all available direct (vectorial) 
methods were designed to analyze only distributions having a single 
grouping; as previously noted, it was not until 1964 that Selby dealt 
analytically with girdle distributions. Accordingly Flinn observed 
that the vectorial methods were of limited application to rock fabrics, 
which more commonly showed several modes, or whole or partial girdles, 
or elongate groupings, to which the vector methods were not applicable. 
He pointed out that the type of test in which a stereographic diagram 


is compared by means of the chi-square test with a model random diagram 
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Table |. Combinations of Properties Observed when Measuring 
Oriented Quantities in Geology, with Examples. 


Number of Dimensions 
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is better able to handle the complexities of a multi-modal, possibly 
non-normal, distribution (Flinn, 1958, p.527). 

Flinn (1958) was concerned only with tests of significance of 
preferred orientation, and not with the calculation of descriptive 
statistics. He described several tests of significance, all of which 
used a stereographic projection, either contoured or uncontoured. 
Based on these various tests, Flinn derived a procedure for testing 
a fabric diagram for preferred orientation. Since this procedure 
incorporates the most useful of the previous tests, only Flinn's 


method will be described here. 


Flinn's Procedure for Detecting Preferred Orientation 


The method is based on comparison of the fabric diagiamto be 
tested with random diagrams, that is, diagrams obtained by plotting 
points drawn at random from a uniform population. The first step 
is to plot the data on an equal-area projection, and then to contour 
the point density by the Mellis method by drawing around each point 
a circle whose area is the area of the projection divided by the number 
of points plotted. For example, a 100-point plot would be contoured 
uSing a one per cent circle. Since a circle is used for contouring, 
the densities are inaccurate in the outer part of the diagram. 

The basis of the tests is the predictability of certain 
characteristics of random diagrams. Two characteristics were used 
by Flinn: (1) The total areas contained between the various con- 
tours. (2) The frequencies of occurrence of different concentra- 


tions, that is, the number of independent and separate areas 
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of one point per unit area, two points per unit area, and so on. 
Both of these characteristics of a random diagram may be predicted 
by means of the Poisson distribution. Flinn also produced 20 random 
diagrams and verified the predicted values. One advantage of con- 
touring with a circle of area A/N, where N is the number of points 
plotted and A the area of the diagram, is that the expected density 
is always the same-regardless of the number of points plotted, and 
therefore the Poisson probabilities of different density values will 
also be constant, resulting in the same expected areas between given 
contours. The frequencies of occurrence of different concentrations 
are proportional to the number of points plotted. 

Briefly, application of the test is as follows. The above 
characteristics are obtained for the fabric diagram under test, and 
the areas and frequencies compared with those expected in random 
diagrams by means of chi-square. Any significant deviation is re- 
garded as an indication of a oreferred orientation, that is, a de- 
parture from complete randomness. 

Flinn (1958, p.531) defined three independent ways in which a 
fabric diagram may show preferred orientation: (1) Grouping: the 
points may be more concentrated than is to be expected in a diagram 


drawn at random from a uniform population. (2) Spatial distributions: 
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independently of whether the grouping of the points departs from random 


expectancy or not, the spatial distribution of points may or may not 
conform to that expected in diagrams drawn at random from a uniform 
population. (3) Structural relations: a single concentration or a 
pattern of concentrations may have a Special geometric relation to a 


structural direction or plane in the rock. This can occur even in a 
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diagram whose points are both grouped and spatially distributed in a 
manner to be expected in a diagram drawn at random from a uniform popu- 
lation. 

The "'area test'' is used to test the grouping of the points. The 
concentration frequencies are used in the ''frequency test'' to test the 
Lae ry distribution of the points against randomness. One final tech- 
nique that Flinn calls the ''two-cell zonal test'' completes the procedure. 
This consists of drawing a circle concentric with the projection that 
divides the projection into two equal areas or cells. In a random fab- 
ric, the same number of points would be expected in each cell, and the 
observed numbers are tested against this expectancy using chi-square. 
Thus if the fabric is first rotated to an appropriate position, the 
test may be used to detect either a girdle or a point concentration. 

In Flinn's procedure, this test is used as one possible test for parti- 
cular patterns after a non-random grouping or spatial distribution has 
been detected. It is also used to dectect a point concentration or 
girdle related to a particular structural direction in the rock, by 
first rotating the fabric until the structural direction to be tested 


lies at the center of the projection. 


The chief purpose of Flinn's procedure is the detection of a 
significant deviation from randomness in a fabric diagram. The two- 
cell zonal test can be used to test for certain patterns. The procedure 
is not concerned with the calculation of descriptive statistics, and 
accordingly provides no powerful means of comparing fabric patterns. 

The method has the advantage that a deviation from randomness 


can be detected before the kind of preferred orientation (the fabric 
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pattern) is known. The method is well suited to the detection of 

very weak preferred orientations (Flinn, 1958, p.531). According 

to Flinn (1958, p.532, 533) it takes about 12 minutes to contour a 
100-point diagram using his method, and about 100 minutes per 100 
points plotted to obtain both the areas and the frequencies of the 
concentrations. Before the two-cell zonal test can be used, the fabric 


must be rotated, but does not need to be contoured. 


Reproducibility as a Test of Significance 


Flinn (1958, p.526) listed three basic kinds of test of signifi- 
cance of preferred orientation: (1) Similarity of successive samples. 
(2) Comparison with preferred orientation. (3) Comparison with random- 
ness. Method 3 is that employed by Flinn, and also by Watson (195€b) 
in his test of the vector resultant, R. A few words will be said here 
about method 1. 

Flinn (1958, p.526) argued that comparison of successive samples 
as a test of significance of a fabric pattern requires the knowledge 
or assumption that the fabric in question is homogeneous. If the homo- 
geneity of the fabric were not first established, differences between 
successive fabric samples could be interpreted as indicating a lack of 
significance of the fabric characteristics involved, or as indicating 
spacial variation of the fabric. 

However, Kauranne (1960), in till fabric studies in Finland, 
tested the reproducibility of fabric samples by the following pro- 
cedure. In a set of 200 observations of pebble axes, all the odd- 


numbered observations were assigned to sample A, while all the even- 
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numbered observations were assigned to sample B, thus producing two 
samples of 100 measurements each. The fabric diagrams of samples A and 
B were then compared to test the reproducibility of the fabric (Kauranne, 
1960, p.91 and Table 1, p.90). As a result of this procedure, any 
Spatial variation in the fabric will appear in both diagrams A and B. 
The two samples A and B are independent duplicate samples of the same 
fabric. If similarities are found, they can be interpreted as indicat- 
ing non-randomness. Any differences must be the result of random 
sampling, since they are not due to spatial variation of the fabric. 
This technique could be made yet more diagnostic by taking larger 
samples and assigning the measurements to three or four groups instead 
of two. 

The above method could be used to establish the minimum sample 
size required to define the fabric with any desired degree of reliability. 
The method would be particularly valuable in the case of, for example, 
a complex fabric not amenable to direct mathematical treatment. Kamb 
(1959, p.1900), in discussing ice-crystal fabrics, emphasized the im- 
portance of reproducibility as a test of significance of fabric pro- 
perties. 

Reproducibility in duplicate samples taken during the present 
study to establish the reliability of samples of 50 long axes is dis- 
cussed in Chapter 4. In evaluating the reproducibility of sample 
characteristics in order to determine whether or not the sampled pop- 
ulation is non-uniformly distributed, it is important that examples of 
random samples from uniform populations aiso be studied. These 
illustrate the degree of reproducibility, or non-reproducibility, to 


be expected if the population is uniform, and thus provide a basis for 
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recognizing a significant degree of reproducibility. 


Direct Methods 


Vector Summation Method 


Fisher (1953) suggested a probability distribution as a basis 
for the statistical treatment of vectors scattered about a mean 
direction. Watson and Irving (1957) employed Fisher's distribution 
in a detailed development of a suitable statistical procedure for the 
treatment of observations of the direction of magnetization of rocks. 
Fisher (1953) suggested that those directions will, when regarded as 
points on a unit sphere, be distributed over the sphere in accordance 
with the function 

C.exp(Kcos ¢) 
where K is a precision parameter, g is the angle between the polar 
vector (the centre of the distribution) and the observation vector, 
and C is a constant. The density thus has axial symmetry about the 
polar axis, and attains a maximum at the pole and a minimum at the 
anti-pole. When K=0 the density is uniform over the sphere. When i 
is large the density is confined to the region about the pole and has 
approximately a Spherical normal destribution in this region (Watson, 
1956a, p.153). 

If (hs oms5.n3.) are the direction cosines of the i-th observation 
of N observations, and (1,m,n) are the direction cosines of the maximum 
likelihood estimate of the polar vector, then Fisher (1953) has shown 


that 
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1 =2 1; ee Oogam i anaeeY, 
R R R 
where R* = (215)? + (x ms)? + (Ting) 4. Thus R is the length of the 


vector resultant of all the vector observations, and (1,m,n) are the 
direction cosines of this vector resultant (Watson, 1956a, p.154). 
Fisher also showed that the best estimate k of K, the precision para- 


meter of the population, is given by, for K greater than 3, 


where ll! is the size of the sample (Watson and Irving, 1957, p.292). 

Watson (1956a) gave approximate tests for the equality of the 
K's and polar directions for any number of populations. These tests 
of significance were elaborated by Watson and Williams (1956). Watson 
(1956b) gave a table of significance points for a test of K=0 using R, 
the length of the vector resultant of the observation vectors. 

Watson and Irving (1957) described the application of the above 
methods to paleomagnetic data. Since the method they described is 
the chief method used in the present study, it will be delineated in 
some detail below. After calculating the best estimates of the polar 
direction and precision parameter, they applied Watson's (1956b) test 
of randomness. This test is based on the argument that, ''Given a 
sample of size N, the length R of the vector resultant will be large 
if the sample shows a preferred direction and small if it does not. 
Assuming that there is truly no preferred direction (that is, K=0), a 
value Ro, say, may be calculated which will be exceeded by R with any 
Stated probability.'' (Watson and Irving, 1957, p.292). Watson (1956b) 


has tabulated Ro for N=5(1)20 and probabilities 0.05 and 0.01. He has 
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also shown that for N>20, Ro is successfully approximated by the 
expression 
a2 
Ns 
3 


where Me stands for a chi-square variate with three degrees of freedoin. 

It should be carefully noted that the above test of randomness 
was designed to test a set of unit vectors whose sample space corre- 
sponds to a sphere. The importance of this fact will be discussed in 
subsequent sections. 

The next step in the analysis of Watson and Irving (1957) is the 
comparison of the precision parameters of different populations. For 
two samples consisting of [!; and ilg observations that give precision 
parameters ky and ko, the ratio k;/kg has the F-distribution with 
2(tlo - 1)/2(ily - 1) degrees of freedom, and may be used to test the 
hypothesis that Kj=K2. Values of kj/ko far from unity suggest that 
K)#K5. For several populations, the ratio of the largest to the smal - 
lest k may be referred to tables of the maximum F-ratio (Hartley, 1950; 
David, 1952), in order to test the hypothesis that K is constant over 
the populations. Next, the accuracy of the calculated mean direction 
is estimated (Watson and Irving, 1957, p.292) as follows. If c is the 
cosine of the angle between the calculated mean direction and the true 


mean, Fisher has shown that this cosine will be less than c with proba- 
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Zp Noe, 
may be referred to the F-ratio tables with 2(p - 1) and 2( Ny - p) 
degrees of freedom. It is supposed that the sample from the i-th 
population contains tl; specimens and has a resultant of length Rj and 
that R is the length of the vector sum of resultants of the separate 
samples. Again, large values of the statistic suggest that the assump- 
tion of identical true mean directions is false because the algebraic 
- sum of the sample resultants R; will then be much greater than the 
length of their vector sum, R.'' (Watson and Irving, 1957, p.293). 

Watson and Irving then described a method of testing the fit of 

the observations to Fisher's distribution, using chi-square. Finally 
they presented a technique that may be used when observations are made 
at several sites, and that permits calculation of within-site and 
between-site precision parameters. This technique can also be used 
to determine, on the basis of initial test samples, the optimum sampling 
pattern to define an overall mean direction with a given accuracy 


(Watson and Irving, 1957, p.298). 


The above method of treating sets of three-dimensional vectors 
is compared with the required operations previously listed. The R 
test for randomness serves the purpose of establishing the presence 
or absence of a preferred orientation. The method is designed to be 
applied when the basic form of the distribution is known to be unimodal 
and approximately spherical normal. Fit of the data to a spherical 


normal distribution can be tested. If the fit is very poor, or if the 
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distribution is bi- or multi-modal, the method cannot be applied as 
described. The method provides statistics describing the location of 
the distribution and its shape; it also allows comparisons to be made’ 
of different distributions. [lence the method takes care of all required 
operations except determination of the kind of preferred orientation, 
that is, the basic forr of the distribution. It is applicable only if 
the distribution is known to be unimodal and approximately spherical 
normal. 

Watson (1960) extended the above techniques to construct a test 
for coplanarity of population means, a property sometimes of importance 
in geophysical studies. (le further adapted the method for the study 
of folded strata (\’atson, 1965). In the latter application, unit 
vectors normal to various portions of a surface having nearly cylindri- 
cal folds are distributed on the sphere in the form of a partial girdle. 
This differs from a girdle in lacking rotational symmetry (Watson, 1965, 
p.193). The problem of girdles has been dealt with by Selby (1964). 

In all the situations just mentioned, the form of the distribu- 
tion must be known or determined before the appropriate method can be 
used, and if the pattern is complex (for example, a point concentration 
and a girdle may both be present) the methods are not directly appli- 
cable. Furthermore, the methods were designed to deal with observations 
of directed quantities, and great caution must be used in applying them 


to a set of non-directed axes. 
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Eigenvector Method 


A technique designed specifically to treat a set of non-directed 
axes was presented by Scheidegger (1965). As in the previous method, 
the procedure is valid only for unimodal distributions. Scheidegger 
postulates a distribution of the form 

p= const .exp(k.cos“¢) : 
where P is the probability at an angular distance ¢ from the central 
axis of the distribution, and k is a precision parameter. For k=0 the 
distribution is uniform. The distribution differs from that of Fisher 
(1953) only in having cos? in place of cosg. This makes Scheidegger's 
distribution a slightly different shape from Fisher's. More important, 
it means that the former is periodic in ¢ with a period of 180 degrees 
(Scheidegger, 1965, p.C1€5), whereas the latter has a period of 360 
degrees. Stated differently, this means that whereas Fisher's distri- 
bution has as its sample space the whole sphere, Scheidegger's distri- 
bution is restricted to a hemisphere. The reason for this difference 
is as follows. Fisher sought to describe the distribution of the 
directions of directed lines, whose sample space is the whole sphere. 
Scheidegger's purpose, on the other hand, was to describe the distri- 
bution of the orientations of non-directed lines, whose sample space 
is a hemisphere since any orientation may be represented in a given 
hemi sphere. 

The problem of finding the mean and variance of a given sample 
may be regarded as the problem of fitting an appropriate hypothetical 
distribution to the sample. The distribution is located in such a way 


as to maximize the probability of the sample. Scheidegger's method 
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operates as follows. Let Pj; be the probability of the i-th axis in a 
sample of N axes, and let the smaller angle (that is, less than 90 
degrees) between the i-th axis and the central direction of the distri- 
bution be g;. The probability of the sample is given by the product 
of the probabilities of the individual axes that make up the sample. 
That is 

Pare Prod (P, ) 
where P_ is the probability of the sample and Prod (P; ) means the product 
of all the P;. Writing this in terms of Scheidegger's density function, 


P 


S Prod (const .exp(k.cos“¢;)) 


const .exp(=ik.cos¢;) 


const .exp(k. =cos“¢, ) 

Clearly, Pg will be a maximum when = cos*6; is a maximum. If 
(1,m,n) are the direction cosines of the central direction of the dis- 
tribution, and (1;, mj, n;) the direction cosines of the i-th axis in 


the sample, then 


cos p+ = ir pte mom. + nen; 
= ety by 2 2 2 
and cos g: = a11! + agom + 233n + 2a1oIm 
+ 2a; 3I1n + 2a mn 
where ree Shh 
= 2 
agg = 2m, 
3 2 
a33 = 2; 


a12 = 291 =P Im, 
813 = 83,.> 2 1.n, 
823 = 832 => min; 
Thus the problem is to maximize the above quadratic in 1, m and n 


under the condition that 17 + m2 + n@ = 1 (since 1, m and n are direction 
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cosines). The solution to this classical problem is that (1l,m,n) are 
the direction cosines of the eigenvector of the matrix a; corresponding 
to the largest eigenvalue v;. This eigenvalue is the corresponding 
value of Lcos*¢;. Hence 

Vpiut dicos2¢. 


N ON 


which is the mean square of the cosines of the deviation angles. This 
mean square cosine may be converted to an angle, which may be regarded 
aS a Standard scattering angle. The direction of the eigenvector cor- 
responding to the smallest eigenvalue is the direction in which Zcos*¢; 


iS a minimum. A corresponding scattering angle may be calculated. 


The method provides an estimate of the mean of a sample, and an 
estimate of scattering. The distribution must be known to be unimodal 
and to conform approximately to the distribution const .exp(k.cos“¢) . It 
would probably be possible to construct a test for randomness using the 
maximum value of Dicos@g., and tests for comparing distributions using 
the statistics given by this method, but it appears that such tests have 


not yet been devised. 


Further Discussion of Fisher's and Scheidegger's Distributions 


It is interesting to compare Scheidegger's (1965) method of 
locating the mean of a set of directions to that of Fisher (1953). If 
the same approach of maximizing the sample probability is taken using 
Fisher's distribution, the problem of locating the mean of a sample 


reduces to that of finding the direction that maximizes 2cos é:. But 
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cos g; is simply the component in the reference direction of a unit 
vector in the direction of the i-th observation, and the sum of the 
components in a given direction of a set of vectors iS a maximum when 
that direction coincides with the direction of the resultant of the 
vectors. Hence Fisher's solution is obtained: that the best estimate 
of the mean of a set of directions is the direction of the vector sum 
of a set of unit vectors in the given directions. 
Mathematically the above can be stated as follows. If (1,m,n) 

are the direction cosines of the required mean direction, and (1;,m:,n;) 
the direction cosines of the i-th observed direction, then 

cos g;= 1.15 + mom; + nen; 
and Di cosepjo= Vey etom.limy otene inj teoF 
and the problem is to maximize F under the condition 12 + m@ + n2 = 1, 
This condition may be used to eliminate one variable from F, making it 
a quadratic in two variables. The solution for the maximum is then ob- 
tained by equating each of the two partial first derivatives to zero. 
This gives 

oul. = _2 mj zn, 


ae Bion R R 


where R* = (2a\e + (Zim, )2 + (21n,)-. Thus (1,m,n) are the direction 
cosines of the resultant of unit vectors in the given directions, which 
is the result obtained by Fisher (1953). 

If one attempts to follow a similar procedure, again using 
Fisher's distribution, for a set of non-directed axes, it is necessary 
to use }cos d; | (the absolute value of cosine ¢;), since a non-directed 
axis always makes an angle of less than 90 degrees with any given 


direction. ow the problem is to maximize eee é;| = Mt + mom; +nn|. 
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This, however, cannot be expanded as was S:cos ¢;, and therein lies the 
real advantage of Scheidegger's distribution. Since the latter distri- 
bution uses cos*¢; instead of cos g;, it makes no distinction between 
cos g; and -cos ¢;, and hence is suitable for non-directéd axes. 
Furthermore, 2'cos“¢. may be expanded as previously indicated, allowing 
the straightforward eigenvector solution. In many computer systems 
subroutines are available for calculation of eigenvalues and eigen- 
vectors of a real symmetric matrix, such as the one employed in this 
procedure. 

One of the computer programs written during the course of this 


study would, with a simple modification, make possible a graphical 


solution to the problem of maximizing =| cos g:| , thereby permitting 
the treatment of non-directed orientations using a model distribution 
having the same shape as Fisher's distribution. Such a procedure 


would still, however, be limited to unimodal distributions. 
The Treatment of Non-Directed Orientations 


Parts of this report are concerned with evaluating a published 
procedure for applying the vector summation method described in a 
foregoing section to the orientations of the a-axes of till pebbles, 
and with devising an improved procedure. It seems logical, therefore, 
to consider at this stage the more general problem of applying the 
vector summation method to any set of non-directed orientations, as a 
basis for later discussions. 

Fisher's vector summation method was designed for estimating 


the mean of a set of directed orientations grouped about a single 
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direction on the Sphere. Any orientations to which the method is to 
be applied must be in the form of vectors or directed orientations. 
To make a non-directed orientation into a directed orientation it is 
necessary to assign a sense, that is, to choose one end of the axis 
or the other. This fact was pointed out by Steinmetz with reference 
to c-axes of quartz grains in metamorphic rocks: ''The c-axes could 
become vectors only if a sense were assigned by some external criterion.'! 
(Steinmetz, 1962, p.801). It was also recognized by Scheidegger 
(1965, p.C164): "A vector is a directed quantity, an axis is not. 
Hence, in order to assign a vector to each axis, it iS necessary to 
make a convention with regard to the side of the axes in which the 
corresponding vector is Supposed to point.!'! 

The only method of assigning sense to a number of axes while 
preserving the characteristics of the original distribution of the 
axes is to cut the distribution in half using a plane passing through 
the origin of the reference axes. In other words, instead of con- 
sidering the distribution as Sccupying a sphere, one considers only 
that part of the distribution that lies in a given hemisphere. It 
should be observed that any plane will cut a set of non-directed axes 
into two identical hemispheres. It is important to remember that, 
while such a hemispherical distribution completely defines the original 
spherical distribution and therefore adequately represents it, it is 
not identical to it, since the original distribution consisted of two 
such hemispheres. 

Thus if a set of non-directed orientations is unimodally dis- 
tributed, that is, the orientations are grouped about a single orienta- 


tion, it is usually possible to obtain satisfactory results by consider- 
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ing only that part of the distribution that lies in a selected hemisphere, 
assigning unit vectors in the observed directions in the hemisphere, and 
then utilizing Fisher's vector summation technique to locate the best 
estimate of the mean. It is only necessary that sufficient care be taken 
in selecting the hemisphere to be used. 

Distributions of non-directed orientations that are not unimodal, 
that is, in which the orientations are grouped about two or more centres, 
cannot be handled in this way. The vector summation method can be ap- 
plied only if the observed orientations are divided into a number of 


groups each of which has a unimodal distribution. 


The Question of Reference Planes 


Much of the confusion that has characterised the literature con- 
cerning geological applications of the various methods of treating 
orientation data has two chief causes: (1) The failure of workers 
applying the methods to recognize the differences between different 
kinds of orientation data, particularly the distinction between directed 
and non-directed orientations and the limitations imposed on the methods 
by this fundamental distinction. (2) The failure of workers applying 
the methods to recognize the distinction between the several different 
functions of reference planes used in the treatment of orientation data. 

The different kinds of orientation data commonly used in geology 
were described previously; this section will describe three separate 
functions performed by reference planes. 

1 When orientations are measured in the field they are referred to 


the horizontal plane; this is the first reference plane. Directions 
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in this plane are usually referred to true or magnetic north. 
These are conventions, adopted for convenience and to ensure 
reproducibility and compatibility of observations. 

Measurements of orientations in two dimensions can be plotted 
directly on a two-dimensional diagram. tlowever, before orienta- 
tions measured in three dimensions can be represented on a 
two-dimensional diagram they must be projected in some way onto 
a selected plane. This projection plane is the second reference 
plane. It need not be the same plane as that used as a refer- 
ence for the field notes; the data may be projected onto any 
plane. However, use of the horizontal plane as a projection 
plane offers certain advantages: (a) the data may be related 
directly to compass directions, (b) directions of dip are readily 
apparent, and (c) orientations measured using a horizontal 
reference plane are most easily plotted when the projection 
plane is also horizontal. 

As described previously, it is possible to use certain treatment 
methods designed for directed quantities in the treatment of 
non-directed quantities provided appropriate care is taken to 
select a suitable origin. In the case of two-dimensional data, 
this origin is a line in the data plane; when the observed 
quantities are three-dimensional, this origin takes the oe of 
a plane. This is the third reference plane. It has nothing 
whatever to do with either the measurement reference plane or 
the projection plane used for display purposes. All three planes 
are distinct and completely independent. There is no reason, 


for example, why three-dimensional data cannot be projected onto 


bestelq od més enotenemtb owt nt enctsesneive Ye etnemeiuessh) So % 
-sinsite e1oted .vavewoll .mergsib lenotdaamtb-owt 6 ao yfsze1tb |! 
6 no beinegerges ed nao edotenemtb seidt nt bewessm 2nott — 


ono yew omoa' nt betzeLorq 8d teum yots meapeth Teno tenemtb-owd tah 


sonsetey baesse oad 2F enetq nottostorg, ett sane! q: betelez il ol 
-19799 6 2B ‘eeu sentt | 26 enstq smee sit od jon bson tI sonst 
xns osno ala ta th at siteb a bist? eds ai sick e 
nottosjorg 6 28 snelq Isfnosimion st AB ocu , 19vSwoH aie 
bstelen od yem eisb ot (6) veapsdnkVbe nistres e1stto enetq’ 7 
yl tbsen ov qib.3o enottoe th (d) .edotdostib 2zeqmoo ot yl toortb 
fetnosred,6 pateu bewsesm enotternetyo (9) brs ,Inotsqqe 
notsoetorg et ner beatotg I Fess teom sis snsiq sonatsteay = = 
| Jletnosirod cals 27 snstq 6 70 
Insmise nt nisies eeu od widtezeq et ti Nievotverg bediroesb 2A” Cas 
to Inemiee ts sit ni zettisneup bedostib yo? benpigeb eborijeam © 
oF raves er ows. eteingorqgs bsbivoyg estttineup betogxtb-nen © 
etsb lsnotedemtbsow? Vo e260 edz nl intetio stdestue 6 testes |) 
. bevieedo sit new yonetq eteb ad? ni snit & et nigtso efits © "% 
20 ano? att 2edet atte eine iene dtanasialshal ssa 
ba rc er S 


é | W | . 
een ey" ay a my id _ 


ia oe | sion “a 


‘age ip wae 


AT 


a horizontal plane for display purposes, while being mathemati - 
cally treated using a vertical plane or any other plane as a 


reference for assigning sense to the orientations. 


Before the advent of electronic computers, a worker wishing to 
adjust orientation measurements on a set of axes so that the Specified 
ends of the axes would fall in some hemisphere other than the conven- 
tional one below the horizontal would have used a graphical procedure. 
The data would have been plotted on a horizontal projection, and, with 
the use of a stereonet, adjusted graphically so that the horizontal 
projection would, in effect, be transformed into a projection of the 
same data on the desired plane. The data would then be ''projected 
back'' jnto three dimensions by reading off the new trends and plunges. 
This is probably the reason why the analysis of data using other than 
the conventional lower hemisphere is inseparably bound, in the minds 
of some workers, to the procedure of "'rotating'' the data to forma 
new projection. However, the new projection is not required for the 
mathematical treatment: it is only necessary that the three-dimensional 
data be adjusted to the desired hemisphere. This can be accomplished 
directly without the use of a new projection and without rotating the 


data. 


The orientation of a pebble a-axis is normally measured as trend 
and plunge with the plunge measured downward from the horizontal. Thus 
the resulting values of trend and plunge actually indicate the direction 
of the downward-pointing end of the axis. It should be remembered, 


however, that the use of the horizontal plane as a reference, and the 
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measurement of plunge downward rather than upward from the reference 
plane are merely conventions, established to permit unambiguous 
reporting of three-dimensional orientations. There is no implication 
that the axis is a downward-pointing directed quantity. Treatment 

of the measurements directly as downward-pointing directed quantities 
constitutes use of the horizontal plane as a reference plane for 
assigning sense to the axes, and this may or may not be a suitable 


way of treating the particular set of data. 


Summary 


Flinn's (1958) procedure is designed to detect preferred orienta- 
tion and test the relationship of preferred orientations to particular 
directions. It is specifically intended for use in marginal cases in 
which the preferred orientation is weak. The method was not designed 
to determine the preferred orientation or any accurate estimate of 
the strength of preferred orientation. 

Scheidegger's (1965) procedure was designed to determine the 
preferred orientation and amount of scattering in a set of non-directed 
axes. Its disadvantages are: (1) It does not provide precise tests 
for comparing samples based on the statistics produced. (2) It treats 
only unimodal distributions having an axially symmetrical grouping. 

Watson and Irving's (1957) procedure was designed to treat vectors, 
and does provide both descriptive statistics and precise tests using 
these statistics for comparing samples. It also requires unimodal 
distributions with axial symmetry. 


Hence the most complete procedure is provided by Watson and 
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Irving's method. However, just as Scheidegger's method was not designed 
to treat vectors, that of Watson and Irving was not designed to treat 
non-directed orientations, and the application of it to, for example, 
till fabrics must not be done without careful consideration of the 
applicability of each step of the procedure. 

It is important to recognize that each of the above methods was 
designed specifically for either directed or non-directed orientations, 
and that a given method may not be directly applicable to different 
kinds of orientation data. Important too is the recognition of the 
three functionally distinct reference planes used in the treatment of 
orientation data: (1) The horizontal reference for measurements. 

(2) The projection plane used for graphical display of the data. (3) The 
plane used for assigning sense to non-directed axes if they are to be 


treated by a method requiring directed lines. 


Methods Examined for Applicability 


to Till Fabric Data 


Previous Application to Till Fabrics 


Only one of the above-described methods has previously been 
used in the treatment of till fabrics: the method described by Watson 
and Irving (1957) for the analysis of paleomagnetic vectors was experi- 
mentally adopted by Andrews and Shimizu (1966) for the analysis of till 
fabrics. They used a FORTRAN IV program which assigned unit vectors 
to the downward-pointing ends of the till stone long axes as specified 


by their azimuth and plunge, and then computed the direction and 
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magnitude of the vector resultant. In addition, it calculated an 
estimate, k, of Fisher's precision parameter K using the formula 


given by Watson (1956a) : 


where N is the number of observations and R the length of the resultant 
vector. The program also provided the 95 per cent confidence radius 

for the calculated mean direction, and the magnitude, R, of the resultant 
vector. The analytical prodedure was described by Andrews and Shimizu 
(1966, p.156) as follows: 

"The approach to the analysis of the computer results is best 
seen through the series of progressive steps listed below: 

"(1) Obtain till fabric data and feed it into the computer 
program. 

'(2) Test the hypothesis that the three-dimensional distribution 
is random by evaluating R against an acceptable confidence level, e.g., 
90, 95 or 99 per cent. 

''(3) Accept or reject the null hypothesis; if ''accept'' the analysis 
must stop, if "reject'' the analysis proceeds to (4). 

"(L) Plot the values for orientation and dip as well as the limits 
of the confidence area which has been determined prior to the analysis. 
This involves rerotation if rotation has already been applied to the 
observations. 

'(5) Ascertain whether the sample approximates a spherical- 
normal distribution, the lower limit of acceptance of this hypothesis 


being k23. If it does not approximate a spherical-norinal distribution, 


then the analysis must cease or an attempt be made to improve the value 


of k by rotation of the data. Implicit in this lower limit of accep- 
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tance is the fact that, with k equal to 3 for 95 per cent of the time, 
observations will lie within a radius of + 90 degrees about the mean 
vector. As k increases this limit becomes smaller. 

'(6) Examine other samples in the same manner and setuc the 
results by use of the maximum F test (see above). If variance is 
homogeneous proceed to (7). 

'(7) Use the F test to ascertain if sample orientations belong 
to the same parent population. 

'(8) Use the results in an analysis of variance design to 
determine precision parameters ® and & (Watson and Irving, 1957).'! 

The three till fabrics analyzed by Andrews and Shimizu (1966) 
were processed using the data as measured in the field, processed 
again after the data had been rotated through 30 degrees about a 
horizontal axis and again after rotation through 90 degrees. It is 
not €lear what axis of rotation was used, but it seems to have been 
a horizontal east-west line. Andrews and Shimizu observed that an 
increase in the calculated R and k values resulted from both of these 
rotations, the maximum values being obtained with the 90 degree rota- 
tion. They concluded that such a 90 degree rotation should be used 
to analyze till fabrics. 

In order to take care of pebbles lying in the B-lineation (that 
is, transverse to the ice-movement direction), Andrews and Shimizu 
(1966, p.160) suggested that, in order to keep sample numbers constant, 
'l0 per cent of the observations lying + 90 degrees from the preferred 
orientation after a 90-degree rotation has been completed, be omitted 
from the calculations and the results rerun,'' 


This method was again used by Andrews and King (1968), with the 
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addition of some calculations based on Watson and Irving's (1957, p.297) 
analysis of dispersion table for the calculation of within- and between- 
site precision parameters and the determination of optimum sampling 
patterns. Andrews and King (1968) regarded the scatter of pebble axes 
in a fabric sample and the variation in the mean directions of different 
fabric samples as equivalent to Watson and Irving's within-site and 
between-site dispersion respectively. 

Andrews and Shimizu's (1966) procedure was used by Cowan (1968) 
and Caine (1968). These works added nothing to the procedure itself 


and will be discussed in a subsequent section. 


First Method Examined 


Since the method employed by Andrews and Shimizu (1966) appeared 
to offer certain advantages over conventional, strictly visual methods 
of evaluating and comparing till fabric diagrams, it was decided to 
investigate its applicability to the data used in the present study. 
The descriptive statistics obtained were examined in relation to 
graphical displays of the same fabric data (the point density diagrams 
obtained previously) in order to determine how meaningful and how useful 
the calculated statistics were. This was regarded as a necessary pre- 
caution before using the statistics in various tests concerning the 
similarity of population means and variances. 

An attempt was made to calculate meaningful statistics to re- 
present the fabric data using the approach described by Andrews and 
Shimizu (1966) with two modifications. Firstly, an improved computer 


program was written that avoided the necessity of rotating the fabric 
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data and, secondly, another computer program was written in an attempt 

to eliminate the arbitrary selection of a reference plane for use in 

the calculation of the statistics. These two modifications are described 
in detail below. 

As discussed in the section on directed and non-directed orienta- 
tions, before any set of nonedirected axes can be analyzed by a method 
designed to treat vectors, the axes must be assigned a sense by choosing 
a reference hemisphere. Frequently, a set of non-directed axes grouped 
about a single direction is transformed into a set of vectors in order 
that the resulting set of vectors may be treated as having a Fisher 
distribution. The closest approximation to a Fisher distribution is 
obtained when sense is assigned to the axes using a dividing plane 
normal to the central tendency of the grouping. 

Andrews and Shimizu used a computer program that assigned vectors 
to the axis orientations as specified. Since the axis orientations 
were Specified by their trend and plunge (measured downward from the 
horizontal), the resulting set of vectors necessarily occupied the 
lower hemisphere. Hence the resulting set of vectors had approximately 
a Fisher distribution only if the axis distribution had its central 
tendency vertical. For this reason, it was necesSary for Andrews and 
Shimizu to rotate each set of axes until its central tendency was nearly 
vertical. 

For the purpose of the present study, a computer program was 
written that calculated the same statistics as those calculated by the 
program used by Andrews and Shimizu, with the addition of the vector 
magnitude, L, equal to R/N. However, instead of assigning vectors to 


the axes directly as specified by their trend and plunge so that the 
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resulting set of vectors occupies the lower hemisphere, the program 
is capable of assigning vectors so that they lie in any desired hemi- 
sphere. The desired hemisphere is Specified by the trend and plunge 
of its axis, that is, the normal to its equatorial plane. Thus, the 
axis of the hemisphere may be made to coincide with the central 
tendency of the distribution. This achieves the desired relationship 
between the data and the reference hemisphere without the tedious and 
time-consuming process of rotating the fabric data. This program is 
referred to as the ''mean vector program''. Details of it may be found 


in Appendix D. 


The second modification to the method described by Andrews and 
Shimizu (1966) concerns the selection of the axis about which the data 
was rotated, and the angle of rotation. In terms of the new program 
described above, this means the selection of the position of the 
reference hemiSphere. Andrews and Shimizu recommended a rotation of 
90 degrees, and in their examples used an eaSt-west axis of ‘rotation. 
No reason was advanced for choosing this axis of rotation, other than 
that it increased the tightness of grouping when the data were pro- 
jected onto a horizontal plane. Only two rotations were tried: a 
rotation of 30 degrees about an east-west axis, and a rotation of 90 
degrees about and east-west axis. The 90 degree rotation was considered 
to be the better, since it produced a horizontal projection with a 
greater degree of cluster. No other rotations were considered. 

In the present study, selection of the position of the reference 
hemisphere corresponds to the Selection of a suitable rotation in the 


procedure of Andrews and Shimizu. One approach to the problem of 
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selecting a good position for the reference hemisphere would have 
been to try two or three possible positions and then choose the one 
that produced a set of vectors with the highest degree of cluster. 
However, this seemed like something of a hit-and-miss approach, and it 
was decided to attempt to devise a method of locating the optimum 
position for the reference hemisphere. The basic aim was to produce 
a set of vectors having a distribution most closely approximating a 
Fisher distribution; in practice this was done by choosing the set 
of vectors that had the highest degree of cluster. Vector magnitude 
(Curray, 1956, p.119-120) was chosen as a numerical indication of the 
degree of cluster of a given set of vectors. A computer program was 
then written that would, in effect, place the reference hemisphere 
successively in 333 different positions, and for each position would 
calculate the vector magnitude of the resulting set of unit vectors. 
The 333 values of vector magnitude were then printed on a stereographic 
projection, each value being located at the point corresponding to 
the pole of the reference hemisphere used in its calculation. The re- 
sulting plot thus provided a graphical display of the variation in the 
tightness of grouping of the vectors as the reference hemisphere was 
placed in different positions. This plot was roughly contoured and 
the peak value located. The position of this peak value indicated the 
position of the reference hemisphere that produced a set of vectors 
having the highest vector magnitude, that is, the highest degree of 
cluster. This position was then specified when calculating the statis- 
tics using the program described above. 

This approach of finding the optimum position of the reference 


plane based on the highest degree of cluster (or minimum variance) was 
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also suggested by Cowan (1968, p.1149) and Andrews and King (1968, 
p.443) but was not employed by them. 

Chayes (1954) put forward the principle of using the 'minimum 
variance origin'' in the treatment of non-directed orientations in two 
dimensions. The method described above is a direct extension of that 
principle to three dimensions. 

This program for producing a plot of vector magnitude variation 
is referred to as the ''vector magnitude program'', Details of it may 


be found in Appendix D. 


Thus the procedure initially used to obtain descriptive statistics 
may be summarized as follows. The data for each sample were processed 
using the vector magnitude program. The diagram produced was roughly 
contoured and the trend and plunge of the peak position determined 
using a Schmidt net. This trend and plunge were then punched onto a 
card and subjected, along with the fabric data, to treatment by the mean 
vector program. The statistics obtained using this method are here 
called ''whole-sample'' statistics and are listed in Table 2. The mean 
direction calculated by this method will be referred to as the ''whole~ 
sample mean direction'', or simply, the 'whole-sample mean''. The 


significance of these results will be discussed in the next section. 
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TABLE 2. Results of first method: characteristics of vector magnitude 
diagrams and whole-sample statistics calculated using reference 
plane giving maximum vector magnitude. ''Plane'' = plane giving 
minimum vector magnitude, specified by dip direction and dip; 

N = sample size; k = estimate of K, precision parameter of popula- 


tion; R = length of resultant vector; L = vector magnitude, R/N. 
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TABLE 2 (continued) 
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Evaluation of First Method 


The Method 


After the statistics for the writer's samples had been calculated 
and compared with the eual-area projections of the data, it became ap- 
parent that the method, even with the improvements described in the 
previous section, had several flaws. Briefly, these are: 

1 The method is incapable of satisfactorily handling multimodal 
distributions, yet these are the rule rather than the exception 
for the a-axes of till pebbles, with which the method is intended 
to deal. 

2 The method misinterprets and misuses the precision parameter, k. 

3 Watson's (1956b) significance points for R, the length of the 


resultant vector, are not valid for non-directed orientations. 


Comparison of the results obtained using the approach of Andrews 
and Shimizu (1966) with the point density diagrams for the same data 
revealed that in many instances the calculated statistics were neither 
useful nor meaningful in terms of the actual preferred orientations of 
the till pebbles. This led to a closer examination of the suitability, 
for the treatment of the a-axis orientations of till pebbles, of the 
method devised by Watson and Irving (1957) for the treatment of paleo- 
magnetic vectors, and to the clarification of certain important differ- 
ences between the two applications. 

Briefly, the relevant aspects of the differences between magnetic 


vectors and pebble axes may be stated as follows. (1) While magnetic 
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vectors are directed, pebble axes are non-directed. This necessitates, 
as described previously, the selection of a suitable hemisphere before 
the orientations can be treated as vectors. (2) While paleomagnetic 
vectors are tightly grouped about a single direction, the a-axes of till 
pebbles are loosely grouped about (in most cases) two or more directions. 
The way in which these differences affect the procedure recommend- 
ed by Andrews and Shimizu (1966) for the treatment of a-axis orientations 
will be discussed below, considering in sequence each of the eight steps 
in their procedure. 
(1) Obtain till fabric data and feed it into the computer program. 
'(2) Test the hypothesis that the three-dimensional distribution 
is random by evaluating R against an acceptable confidence 
level, e.g., 90, 95 or 99 per cent. 
''(3) Accept or reject the null hypothesis; if '‘accept' the analysis 
must stop, if 'reject' the analysis proceeds to (4),!! 
Watson's (1956b) test for randomness of directions was designed 
for vectors distributed on the sphere according to Fisher's probability 
function (Watson, 1956b, p.160). In this situation, non-randomness 
would appear as a tendency for the observed vectors to cluster toward 
one side of the Sphere, rather than being distributed over the whole 
Sphere. 
Pebble axes are non-directed, and, before they can be treated 
using a vector method, must be "'cut in half"! to lie in a selected hemi- 
Sphere. Thus, what is actually dealt with vectorially is a set of vectors 
distributed over a hemisphere. But a set of vectors lying in a hemi- 
Sphere constitutes a clustering or preferred orientation when compared 


to randomness over the whole sphere. Hence any set of vectors ina 
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hemisphere should yield a significant value of R, indicating non- 
randomness on the sphere. That is, a set of randomly distributed non- 
directed orientations would, when sense was assigned by selecting a 
hemisphere, give rise to a Set of vectors distributed randomly over 
the hemisphere. However, this set of vectors would yield a significant 
value of R, since it is not randomly distributed over the sphere. 

If R is the magnitude of the resultant of a set of N unit vectors 
distributed randomly on the sphere, R iS approximately distributed as 


NXi3 

fer 
where x stands for a chi-square variate with three degrees of freedom 
(Watson, 1956b, p.161). For N=50, the 95 per cent, 99 per cent, and 
99.5 per cent points of this distribution are respectively 11.4, 13.7 
and 14.6, This means that if 50 unit vectors are selected randomly 
from a population of unit vectors that is uniformly distributed over 
the sphere, the probability is only 0.5 per cent or one in 200 that the 
magnitude of the resultant will exceed 14.6. That this is not true in 
the case of a distribution in a hemisphere is demonstrated by samples 
505, 506 and 507 (fig. 5), each consisting of 50 orientations selected 
randomly from a population which was uniformly distributed over a hemi- 
sphere. The R values for samples 505, 506 and 507 are respectively 28.6, 
28.8 and 28.1. 

Clearly, the use of this test in the procedure of Andrews and 
Shimizu (1966) is not valid. In any case, the correct approximation for 
nis 

2 
NX3 
3 


(Watson, 1956b, p.161) 
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fitig (Andrews and Shimizu, 1966, p.156). 
Ss 


(4) Plot the values for orientation and dip as well as the limits 


not 


of the confidence area which has been determined prior to the 
analysis. This involves rerotation if rotation has already 
been applied to the observations. 

'(5) Ascertain whether the sample approximates a spherical. normal 
distribution, the lower limit of acceptance of this hypothesis 
being k23. If it does not approximate a spherical normal 
distribution, then the analysis must cease or an attempt be 
made to improve the value of k by rotation of the data. Im- 
plicit in this lower limit of acceptance is the fact that, 
with k equal to 3 for 95 per cent of the time, observations 
will lie within a radius of + 90 degrees about the mean 
vector. As k increases this limit becomes smaller.'' (Andrews 
and Shimizu, 1966, p.156) 

This step involves an incorrect interpretation of the precision 
parameter, k. When the distribution is known to be spherical normal 
(this implies unimodality) k gives information about the shape of the 
spherical normal distribution, that is, the degree of scattering of the 
data. k gives no information about the kind of distribution; in fact, 
if the distribution is not spherical normal, the parameter k, as used by 
Fisher (1953) and Watson and Irving (1957), is meaningless. The value 
of k cannot be used to distinguish spherical normal distributions from 


distributions of other kinds, since a Spherical normal distribution may 
have any value of k from zero to plus infinity. 


Hence, k < 3 does not imply that the distribution is not spherical 
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normal, nor can k>3 be used to justify the assumption that the distri- 
bution is spherical normal. The assumption of a spherical normal distri- 
bution must be justified in terms of the purpose of the entire procedure, 
and whether the assumption aids or defeats that purpose. This question 
is discussed fully later. 

What Andrews and Shimizu intended in step (5) is revealed by the 
second two sentences of the step. The significance of the value k=3 is 
that a Fisher distribution with k >3 has more than 95 per cent of the 
total probability within 90 degrees of the mean direction, that is, more 
than 95 per cent of the probability lies in one hemisphere. As mentioned 
previously, a set of non-directed orientations gives rise to a set of 
vectors that is confined to a hemisphere. If this is to be approximated 
by a Fisher distribution, it is clearly necessary that the approximating 
distribution be essentially confined to a hemisphere, that is, it must 
be essentially zero over the hemisphere farthest from its mean direction. 
If this condition is deemed to be satisfied by a Fisher distribution 
having at least 95 per cent of the total probability within one hemispher, 
then the minimum permissible value of k may be readily computed. Watson 
and Irving (1957, p.290) have shown that the probability that a direction 
will be observed which makes an angle of ¢9 or more with the true mean 
direction is, when k>3, given with good accuracy by the formula 

Pr(g>¢,) = exp } -k(1-cos $o)} 
If $,=90 degrees, 

Pr(g > 90 degrees) = exp { -k | 
Hence, if Pr(¢>90 degrees) =0.05, then 

0.05 = exp { -k | 


and k=3. This result must be regarded as only approximate, since the 
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formula used is not accurate for k< 3. 

Thus, Andrews and Shimizu were correct in saying that the pro- 
cedure must stop if k<3. However, the true reason for this is the one 
just described, and has nothing to do with the validity of the assump- 
tion of a Spherical-normal distribution. 

''(6) Examine other samples in the same manner and compare the 
results by use of the maximum F test (see above). If 
variance is homogeneous proceed to (7).'' (Andrews and 
Shimizu, 1966, p.156) 

In principle this step is valid, assuming the statistics used 
really represent the samples. However, Andrews and Shimizu incorrectly 
stated that the precision parameter k ''measures homogeneity of variance 
and can be tested by the maximum F ratio (Hartley, 1950) where the 
degrees of freedom are 2(N-1) and all samples have an equal number of 
observations.'' The parameter k represents the variance of the sample; 
the ratio of the largest k to the smallest k in a set of samples is 
the 'maximum F ratio'', and it is this that measures homogeneity of 
variance and can be referred to F-tables, not k itself. Furthermore, 
Hartley (1950, p.312) stated that when estimates of variance are based 
on different degrees of freedom (that is, the samples have different 
numbers of observations) the maximum F ratio can still be used as a 
"rough indicator of heterogeneity!’ (Hartley, 1950, p.312) if the mean 
number of degrees of freedom is used and provided the degrees of freedom 
"do not differ considerably". 

'(7) Use the F test to ascertain if sample orientations belong 

to the same parent population.'' (Andrews and Shimizu, 1966, 


p.156) 
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This step is also valid in the treatment of a-axis orientations, 
provided the correct statistic 


(N-q) (ZR; -R) 
(q-1) (N- ZR; ) 


where N is the total number of observations, q the number of samples, 
Rs the length of the resultant of the i-th sample, and R the length of 
the resultant of all the observations combined (Watson and Williams, 
1956, p.348) is used. 
'(8) Use the results in an analysis of variance design to 
determine precision parameters © and & (Watson and Irving, 
1957)." 

This is a valid and useful procedure when applied to sets of 
orientations that have Fisher distributions about their means, the 
purpose of which is the estimation of the overall mean direction and 
the prediction of the optimum sampling pattern. However, it must be 
remembered that ''This treatment is approximate, being valied for smal] 
dispersions only, that is“ and G targe" (Watson and Irving, 1957, 
p.297), unless the modified procedure for situations involving variable 
© is used (Watson and Irving, 1957, p.299). It is not a valid pro- 
cedure when applied to sets of orientations that do not have Fisher 
distributions about their means. Andrews and Shimizu did not test the 
fit of their data to a Fisher distribution. It will be shown below 
that a model of a single Fisher distribution is invalid for all but 


a very small proportion of till fabric samples. 


The biggest single flaw in Andrews and Shimizu's (1966) approach 


of subjecting till fabric data to Watson and Irving's (1957) procedure 
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for the statistical treatment of paleomagnetic vectors is that the 
observed orientations of pebble a-axes are treated as though they 

had a Fisher distribution in a selected hemisphere. However, Andrews 
and Shimizu (1966, p.156, 159), Andrews and King (1968, p.446), Kauranne 
(1960, p.94), Holmes (1941, p.1312), and Harrison (1957, p.295), using 
data from such widely scattered areas as northern England, central 
Finland, and northeastern U.S.A., have commented on the prevalence of 
distributions with both parallel and transverse maxima. Holmes (1941) 
further suggested that pebbles of certain shapes tend to align them- 
selves obliquely to the ice movement direction, thus implying the 
possibility of a-axis distributions having more than two maxima. 

The effect of assuming a Fisher distribution when the actual 
distribution is bimodal can be illustrated by the resuits obtained 
using Andrews and Shimizu's (1966) procedure initially applied to the 
data. For example, sample 1 (fig. 6) shows two maxima at nearly 90 
degrees. The most reasonable interpretation of this, in the light of 
both previous work and the results of the present study, is that the 
stronger maximum at 352 09 represents pebbles aligned parallel to ice 
movement and the weaker maximum at 064 25 pebbles transverse to the 
movement direction. The whole-fabric mean calculated on the assumption 
of a Fisher distribution lies at 017 09 and clearly does not give a 
good indication of the ice movement direction. 

Sample 5 (fig. 6) is another example of a good parallel/transverse 
type of distribution. The two maxima lie at 034 00 and 139 08, the 
first of these being the stronger. The whole-sample mean based on the 
assumption of a Fisher distribution is 184 03, again failing to indicate 


the movement direction. 
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A third good example of this type of distribution is sample 10 
(fig. 7). Here the parallel maximum is at 009 09 and the transverse 
maximum at 104 00. The whole-fabric mean is 041 07, about halfway 
between the two maxima. 

In some samples of this type, for example, sample 7 (fig. 6), the 
whole-fabric mean does lie close to the parallel maximum. However, this 
occurs only when the proportion of transverse pebbles is quite small and 
the distribution has almost perfect symmetry about the plane passing 
through the parallel maximum and normal to the transverse maximum. 

One of the aims of treating till fabric data by the vector summa- 
tion method is to determine the ice movement direction as indicated by 
the strongest preferred orientation of the long axes of the till pebbles. 
The above results show that the effect of transverse pebbles, that is, 
those lying in a weaker preferred orientation at 90 degrees to the ice 
movement direction, is to distort the calculated mean away from the 
direction of maximum density of observations so that it lies somewhere 
between the two preferred directions. Thus, if it is required that the 
calculated mean coincide with the stronger preferred direction, (that is, 
the parallel maximum), then pebbles lying in the weaker preferred 
orientation (that is, the transverse maximum) must be deleted before 


the mean is calculated. 


Andrews and Shimizu (1966, p.159-160) mentioned the problem of 
transverse pebbles, and proposed the solution ''that 10 per cent of the 
observations lying + 90 degrees from the preferred orientation after a 
90-degree rotation has been completed, be omitted from the calculations 


and the results rerun.'' They reran their samples using this procedure, 
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obtaining increased values of k and R/N, and concluded that ''Clearly, 
all statistics increase in significance." 

As nearly as can be determined from the statement of the procedure 
quoted above, and considering the computer program used by Andrews and 
Shimizu, it appears likely that those pebbles were deleted whose azimuths 
lay closest to the directions 90 degrees from the azimuth of the calcu- 
lated mean vector. This procedure has two minor drawbacks: (1) the 
actual proportion of transverse pebbles in tills varies greatly from 
the assumed 10 per cent, and (2) if the mean vector is not horizontal, 
deleting the pebbles on the basis of a 90-degree difference of azimuth 
does not delete those pebbles that are actually 90 degrees from the mean 
vector. This error becomes more serious as the plunge of the mean vector 
increases. The procedure also has one major drawback: if there are 
transverse pebbles, as is usually the case, the whole-sample mean direc- 
tion will frequently not coincide with either maximum of the distribu- 
tion, as described and illustrated above, and hence the transverse 
pebbles do not lie at 90 degrees to the calculated mean direction. In 
short, deletion of transverse pebbles cannot be carried out until the 
transverse direction has been located, but the directions of the fabric 
maxima cannot be located by Andrews and Shimizu's procedure until trans- 
verse pebbles have been deleted. The statistics k and L must increase 
in value if pebbles lying close to 90 degrees from the calculated mean 
are deleted and the data rerun; this does not imply that the new mean 
is any closer to the stronger maximum of the distribution. 

Andrews and Shimizu suggested that k< 3 indicates the presence of 
enough transverse pebbles to be a nuisance and warrant deletion. This 


is revealed by their step (5) in which the value of k is taken to indicate 
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how well the data fit a Fisher distribution, and by the following 
statement in their "Summary and Conclusion": ''in some samples there 
is a large number of transverse pebbles which decrease the precision 
parameter to <3.!'' (Andrews and Shimizu, 1966,p.161). 

Several of the samples used in the present study yielded a value 
of k<3 on initial processing using Andrews and Shimizu's procedure. 
These samples were rerun with all observations lying more than 60 de- 
grees from the calculated mean deleted. Although an arbitrary angular 
cutoff was used instead of an arbitrary percentage cutoff, the pro- 
cedure is basically the same as that used by Andrews and Shimizu, 
except that the true deviation angles were used instead of azimuth 
differences, thus making the procedure more precise. The results ob- 
tained before and after deletion are listed in Table 3. 

In every case, v] 95 (the 95 per cent confidence radius on the 
position of the mean) decreases and both k and L(=R/N) increase. 
However, it is interesting to examine the effect of the deletion on 
the positions of the means. In all cases, the means move by only a 
few degrees, so that if the mean fell near a maximum before deletion, 
it remained near the maximum after deletion, and if the mean was not 
close to a maximum before deletion it was not close to one after dele- 
tion. In other words, the only significant effect of the deletion was 
to increase the statistics k and L and decrease 0 95° 

Andrews and Shimizu appear to have been satisfied by this in- 
crease in numerical value of the parameters k and R/N. However, the 
reason for deleting transverse pebbles is to eliminate their effect 
of distorting the calculated mean away from the region of maximum density 


of observations. Andrews and Shimizu failed to ascertain whether or not 
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this objective had been achieved; the results described above and 
listed in Table 3 indicate that for the data used in the present study 
(apparently not abnormal for till fabric data) Andrews and Shimizu's 
method of deleting ''transverse'' pebbles is unsuccessful. 

Thus, if Andrews and Shimizu's recommended procedure is fol lowed 
for the nearly-ideal parallel/transverse patterns of samples 1, 5, 7 
and 10 (figs. 6, 7), the unrealistic mean of 017 09 for sample 1 and 
the unrealistic mean of 041 07 for sample 10 will both be accepted 
since in both cases k>3. The unrealistic mean of 184 03 will be re- 
jected for sample 10 since k< 3, but after ''transverse'' pebbles have 
been deleted, k >3 and the still-unrealistic mean of 181 02 will be 
accepted. The realistic mean of 046 06 for sample 7 will be accepted 
since k >3. 

These samples represent nearly ideal simple fabrics, yet in three 
cases out of four the mean calculated by Andrews and Shimizu's method 
fails to indicate the ice movement direction or have any other useful 
meaning. 

Table 4 shows a number of hypothetical a-axis distributions of 
till pebbles, their possible appearance in a sample, the useful informa- 
tion that can be obtained from the sample, and whether or not Andrews 
and Shimizu's method of representing the sample can be expected to 
reveal this information. The only cases for which the method is useful 
are those in which all maxima of the distribution are parallel maxima 
produced by ice movement and only the mean direction of movement is 
required. In all other cases some other method of representing the 


distribution must be used. 
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Vie: 
Applications of the Method 


Andrews and King (1968) employed the method of Andrews and 
Shimizu (1966) in studying till fabric variability in northern England. 
They carried the method through to the stage of calculating the within- 
and between-site variability parameters \ and & to determine the optimum 
sampling pattern. 

The fallacious tests based on R and k are still used; however, the 
method of selecting an axis for rotating the data was an improvement 
over the fixed east-west axis used by Andrews and Shimizu (1966). A 
90-degree rotation was used, ''with the axis of rotation passing through 
the area of minimum observations'' (Andrews and King, 1968, p.444). This 
is a better attempt to approximate the ''minimum variance'' position of the 
reference hemisphere. 

Because of the problems already discussed, Andrews and King's 
calculation of w and é and the estimation of the optimum sampling pattern 
cannot be regarded as valid. Andrews and King's claim that their calcu- 
lated optimum sampling pattern is ''surprising but is statistically 
sound" (Andrews and King, 1968, p.457) is itself surprising in view of 
their own statement that the calculated means and k values do not re- 
present Fisher distributions or even unimodal distributions (Andrews and 


King, 1968, p.446). 


In addition to using the method of Andrews and Shimizu (1966) 
(with a fixed axis of rotation), Cowan (1968) determined the directions 
of fabric maxima using a x* test ona plot of a-axis trends, as well as 


plotting point diagrams on polar projections. No fabric diagrams were 
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presented in the paper, however; only the statistics calculated using 

"ea on the trend distribution and Andrews and Shimizu's (1966) method on 
the three-dimensional distribution. Although Cowan was thus able to 
check the mean directions given by Andrews and Shimizu's (1966) procedure 
against graphical representations of the data, he was misled into paying 
unwarranted attention to the mean directions given by Andrews and 
Shimizu's method by the supposedly high degrees of significance of the 
parameter R: ''Orientation and vector strengths indicate that none of 

the fabrics were random and therefore some explanation should be attrib- 
uted to the mean orientations derived from these.'' (Cowan, 1968, p.1154). 

Cowan was misled also by Andrews and Shimizu's (1966, p.156) in- 
correct statement concerning k: ''This statistic k, measures homogeneity 
of variance and can be tested by the maximum F ratio.'' Led by Andrews 
and Shimizu's use of the condition k 23 (Andrews and Shimizu, 1966, p.156) 
into believing that values of k > 3 had some kind of statistical signifi- 
cance, Cowan concluded "'k values (where k = (N - 1)/(N - R)) for the 
homogeneity of variance are significant for both proximal and distal 
fabrics in this example. This means that they may be tested further, 
if necessary, to determine whether or not the mean directions of popula- 
tions are identical’ (Cowan, 1968, p.1150). It is the ratio of k values 
that tests homogeneity of variance, not individual k values. 

Since the numerical procedure for representing till fabric data 
used by Cowan (1968) is known to be of little value in many cases, the 
reader is left without any reliable indication of what the observed 
distributions were actually like. This situation emphasizes the neces- 
sity of presenting till fabric data graphically, at least until a numeri- 


cal procedure has been devised and widely accepted that is valid in all 
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cases and adequately describes the distribution. 


Caine (1968) must be given credit for recognizing the inability 
of Andrews and Shimizu's (1966) method to locate the preferred orienta- 
tions when significant numbers of transverse pebbles are present, and 
for apparently attaching no significance to the calculated value of R. 
In applying the method to the orientation of periglacial blockfield 
material of Tasmania, Caine (1968) used stereographic projections and 
the results of two-dimensional analysis to determine the preferred 
orientations of the blockfield particles. Hethen applied Andrews and 
Shimizu's (1966) procedure separately to the observations lying in the 
two preferred orientations (parallel and transverse to the local slope 
direction), separating the two by a constant 30 degree/60 degree cutoff, 
after rotating the data 90 degrees about a line normal to the slope 
direction. 

Caine cited the results of the three-dimensional analysis as 
tending to corroborate other evidence that the blockfield material was 
preferentially oriented parallel to the local slope direction. However, 
suppose Caine's procedure is followed using a set of random orientations. 
The data are rotated 90 degrees about a horizontal line perpendicular to 
the local slope direction, and the mean of the orientations then having 
trends + 60 degrees from the local slope direction calculated. This 
mean will lie close to the vertical. When it is rotated with the data 
back to the original position, it is very likely that the final position 
of the mean will be close to the local slope direction; the procedure 


followed ensures that this will be so. 
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Second Method Examined 


The first method sets up a single Fisher model for a whole till 
fabric sample. This approach fails because most till fabrics are multi- 
modal. It is suggested that if the Fisher model is to be used to 
describe till fabric samples, the modes must first be located and a 
Fisher model then constructed for each separate mode. An attempt to 
design such a procedure is described below. 

A procedure for numerically representing a distribution of orien- 
tations must be based on (1) the nature of the distribution and (2) the 
information required. For the long axes of till pebbles the relevant 
considerations are: (1) The long axes of till pebbles characteristical ly 
exhibit a high degree of variability in their three-dimensional distri- 
bution pattern, even over relatively short distances within a single 
till unit. This variability is evident not only in the preferred orien- 
tations of the long axes but also in the number of preferred orientations 
and the angles between them. (2) As indicated in Table 4, in most cases 
the information required is given by the positions and the relative 
strengths of the modes of the distribution. 

Thus, a procedure for treating the data should (1) locate the 
important modes of the distribution, (2) determine the relative strengths 
of the modes, (3) describe the modes numerically in a way that permits 
statistical comparisons with other samples. A procedure for statistical 
representation and comparison of samples must be based on a model that 
is at least approximately valid for all samples tested. It has been 
demonstrated in foregoing sections that a model of a single Fisher dis- 


tribution for a sample of long-axis orientations of till pebbles is 
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rarely valid and a procedure based on this model does not give the in- 
formation required. Any unimodal model would be similarly invalid. 

None of the methods presently available for calculating descriptive 
Statistics for orientation data is directly applicable to bimodal or 
multimodal distributions. Even if the methods were so applicable, they 
could not be applied until each sample had been appraised visually to 
determine which model was valid for that sample, since no statistical 
methods are available for testing a three-dimensional distribution of 
orientations against models any more complex than a single Fisher dis- 
tribution. 

Thus none of the numerical procedures available is suitable for 


treating directly a sample of long axis orientations of till pebbles. 


It is suggested that bimodal and multimodal distributions can be 
treated by a method based on a unimodal model if the observed distri- 
bution is first subdivided into a number of separate unimodal distri- 
butions. Such an approach formed the basis of Andrews and Shimizu's 
(1966) suggested technique of separating the modes of the distribution 
by vertical planes using a two-dimensional plot of the trend distribu- 
tion. This would be adequate for the purposes of a two-dimensional 
treatment, but it is felt that when a three-dimensional treatment is 
being used, a method of subdividing the observed distribution based 
on its three-dimensional characteristics is far more suitable. A pro- 
cedure will be described whereby this can be accomplished easily and 
effectively. 

The basis of the suggested approach is the hypothesis that the 


observed bimodal or multimodal distribution can be approximated by two 
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or more superimposed Fisher distributions, each having a different mean 
position, different K, and different relative magnitude. A generalized 
procedure for dealing with such a distribution will first be described. 
The validity of the model for till fabric samples will be discussed in 
the following section, ''Evaluation of Second Method''. 

The first stage of the suggested procedure is the production of a 
contoured point density diagram. This can be accomplished using the 
computer program described in Chapter 2 and given in Appendix D. The 
second stage involves examination of the point density diagram and the 
selection of a model for the sample. This step is discussed fully in 
Chapter 4. Basically the procedure is one of comparison with various 
models, and selection of the simplest model that is consistent with the 
sample. The selected model will have zero, one or two modes if the 
sample size is 50 (50 measurements are not adequate to indicate reliably 
more than two preferred directions). Using a stereonet, the trend and 
plunge of the approximate centre of each mode can be read off, together 
with an angular radius defining the range of the mode on the surface of 
the reference sphere. These values thus define a circle on the sphere 
(or a cone whose apex lies at the origin of the reference axes) that 
contains the observations that define the mode. Such a circular or 
conical region is thought to be the most realistic way of delimiting 
each mode, since the model to be used for the mode is a Fisher distri- 
bution, which has axial symmetry. Further, such a delimitation is very 
suitable for use in the computer program that constitutes the next stage 
in the procedure. 

The mean vector program described earlier with a slight modi fica- 


tion, is now used. As previously described, this program calculated the 
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position of the mean and other statistics on a set of orientations using 
the specified hemisphere. The hemisphere was specified by the direction 
of its central axis. This program was modified by adding an angular 
limit, @(< 90°). The modified program now calculates the same statis- 
tics, but uses only those observations making an angle not greater than 
QO with the specified axis. Thus it is now possible to represent those 
observations lying within a cone with its axis in any desired position 
and with any desired semi-angle, QO (that is, observations lying within a 
circle of any radius <90 degrees in any position on the sphere). The 
number of observations lying in the specified region is also recorded, 
thus providing an indication of the relative strengths of the modes. 
After these statistics have been calculated for each mode, they may be 
used in the relevant tests described by Watson and Irving (1957), pro- 
vided the model can be shown to be valid. Thus all the requirements 
stated above have been satisfied: the modes have been located (graph- 
ically), an indication of the relative strengths of the modes has been 


obtained, and statistics calculated for comparison with other samples. 


Evaluation of Second Method 


The validity of comparisons of samples using the statistics 
calculated by the method just described depends upon the validity of 
the assumption that the observations selected by the specified cone 
actually have a Fisher distribution about their mean. So that this 
could be tested, the mean vector program was further modified to incor- 
porate a test of fit of the data to a Fisher distribution with the 


estimated parameters. Details of this modification and the test of 
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fit may be found in Appendix D. 

Tests were run on 78 selected modes on 42 different diagrams. The 
cones used and the results obtained are listed in Table 5. The first 
point to be observed about these results is that 40 of the 78 cones con- 
tained less than 20 observations, too few to calculate a chi-square 
value for either the polar or axial distribution. A further 20 contained 
more than 20 observations (permitting calculation of a chi-square value 
for the axial distribution) but still too few for calculation of a value 
for the polar distribution. Hence only 18 of the 78 cones contained 
sufficient observations to permit calculation of chi-square parameters 
for both polar and axial distributions. 

If the hypothesis being tested is rejected when a chi-square value 
is obtained that exceeds the value at the 95 per cent point, every one 
of the 18 sets of observations that could be tested would be rejected on 
at least one count (that is, either the polar test or the axial test or 
both). If the 99 per cent point is used as the criterion of rejection, 
13 of the 18 would still be rejected on at least one count. 

Two facts become evident from these results: (1) For this method 
to be useful, much larger samples would be necessary than were used in 
the present study. (2) Few of the modes tested even approximate a 
Fisher distribution. This implies that statistical tests using the 
sample means and k values calculated are not valid. in view of this, 


no such tests were attempted. 
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TABLE 5. Statistics calculated using second method and results 


of test of fit to a Fisher distribution. Column headed 
'PrZ" gives the probability of obtaining the associated 
value of X* under the hypothesis being tested. Four points 
were used: 5%, 2'3%, 1%, and 3%. Each probability is stated 
as 'less than'' the smallest of these four probabilities that 
exceeds it. Probabilities greater than 5% are so indicated. 
All x2 values have one degree of freedom. A = semi-angle 

of cone. n = number of observations within cone. CS = 


chi-square. 


8 | 


cee * stat an 


xa we nix, a a 


iy ia "teas ent hie 


batele 2 vai idee ned om pag: ee et iS aks 
pa at : re we 
todd cot etree “wot: seers Ae ioe ta one ae a ! 


i 
‘oh : a 
-boiso roar: 02. ore Ke nae Sateong weit Led 
desde tye 92 ree yao) to nenpab ono aves desis 
Vey +2009 4 ateiatw ‘ariottavisado. 4 Setlea pleas 


me f = fa 8 if ‘ - sys 49 


oa oh A. aN " . 
ioe nel 7 
| 


STA 


ae 
oO 


O10 DO ONN WWW FPwwhrd Wd — I 


028 


———=> ome coe 


RESULTANT 


CONF. 


82 


022 


° ° ° ° ° ° ° ° ° 
“SO Unk iE t— ENO EG) 


° ° ° ° 
W Ww CW 


NNi CO RE NMW FENWO WW 


° 
— 


— 


— 


COO 0 — OO W1o NO © —~S 
° ° ° ° ° ° ° ° ° ° c) ° 


OOo —- HnNOWN —- WN O 


CO CO © 


Oxn~s 


— (oe OO COO 10 CO (oe) 


\O CO ON 


\O NI \O 
° fo} e 


MwW™NW 


\o 


ow fe 


CO: BS G6. SO. €6y oO 40m gO GO 0° ge 


° 
— 


° 
CO 


° ° 
NO 


oom 6 
Ww NW 


14. 


he 


° © 6 «6 SO. go. 40 “oO ec ° 
—) (Oe) NNO CO EN CO = 


Nwu— Oo 


\o © \O OVO NUT WOO ON — W WOV’O 


° ° 


° 
WW 


PESPLOF: FIT 
POLAR AXIAL 
CSE PRC Sas PR, 
203 pou 553 O25 

B28 <2'5 
3201875 
5.8 <2 
45 65 
8.4 <4 11.5 <b 
heh & 45.9 ls 
eo US Ala e < ls 
A-8 die 9.6 6% 
4.2 <5 
Lay a5 WV eS <5 
Ach oe - 6.9 <1 
10.2 <5: 26.5 < 
2.8 »5 4.9 <5 
Bu3 “2's 46.6 (2 
Poa ss 320.220 
14,8).<% 
9.0 <b 
2-025 aS tis 
5.6 <23 lO ce 
5.5 <2 
6L8e 


pred ae "aaa aa 
i> ee & - se a 8 Se on 


28> B. 3 } a : aan ae 8 qe 9.8 Bia! OLS cay 80 ’ 

aN Se ee 8 A 

2S) 8.8" ; ” “Pe TS, ee 7 tS 80. 00° 4 
2) aH baad PLO LTT % e SS 8! €ep ie 
| ene M8.81'@.07 8.8 (dl so ES 4 

y.#e 


4 Yede QBN 2 Ge Step Bel Bi oe 
o> TT HB A108 20.8 esi ie Naot eek to ae eas 
a> Oat Po 8.88 eel Pe hey Be. 30 eo an oe 
o> S.af @ Be f.8y Tee at Buh VASE Asal Xa (08 @ ESE 
: if i vara le pas re = He ‘ RFR Sate hap egy: 
> 3.8 % 8. 3.08 O.0f Ae 2.) ae eo so. a ce 
a> S.A 1.82 EX. @b 8.8 oe it 80 ond- 28 o8 @0 Ol 
2.20 88.08 a.0f £8 ROR) fo Sol dt RE 00 | 
S.2e ae.e! Sif Se a ef ass af Of ar’ ss 
J “i eae | ee. ote 
“OT. ah: at Ho OVS St a So SVs 
TN od ae a.) see | { ae Af ee0 
a> €.8 2@ Tit @.@8 T6.8 O.0f 058 EVE )S0 0S TE G2 £0 O80 a 
Bs) Bde KY.TF AST VR EET 00 VSO #f 2 oo Le 
fy 9 Re \. S Cle Sr eB «OSE gol ae! ae a Bo Ter & 
| 18 B@eS Bill Le sft Vl OVS St es ot eas 1 i 
| SNE 20.TE EN T.8  f.@, BOSS! OF ef 30 0S! 
ft? ris as. Mia ae of $.€8 $0.0. oft ER . PRE as SSE . ve ot TE 08 = 
a ed Re 8. $ €.68 40.8 a€1°8,01 Aves 20 ade ef 08 £0 SOE. OS 
RRO OOL2F Bet ay f'E. H 00 280° St 26 ant 7 os 
aS a8 3? C2 BG SSF Hor S.8S OF EO OF i ig 
: Bae AROS BATT OR #,€1 @0 Geo A yO 
(.a¢ fe.vS Bot 6.8 8.01 a hee we: 
nla wee 2 O08 0 
| 0 Ss AP 8 
e< €.€ ar ie 
eh Bal 
gle £8 
H#> [.00 8 


83 


TABLE 5 (continued) 
CONE TEST OF FIT 
SHA» PAXTS RESULTANT  CONF.RAD. POLAR AXIAL 
NOe wi Re Pleas in “TRo PE PR Ob 99h Koh) CSe PRE CS: (PRA 
Smee Ose On toc) weno na 1956 1 69. Salat 1210.92 31 Beto 
SU, 105 lo 25 9 
Bi 05000125 10, 229201 9.6. 10.0 13.0. 24.53.96, 3 
31 350 18 30 12 346 18 11.6 8.9 11.5 24.69 96.3 
Bee leu20 630) 16009 23.1025 10.1. 13.0.21,.53. 95.8 
BZ, a5 12 40 22 "toc Is 2Oel “oro sl2cd roroo ore 2S) 5 
B50 209 29720. 9 
Bor 559 TS Hey) 4 Ty oe Gro =la.o COTO? SOs 5.1 <2& 
34 024 10 26 29 024 09 28.1 Ge “Go0 32 3197-0.) be. 45° 2 3932 5 
3425 O24 EH 271 BS 02507 +17, FE 47 IPM O20 P24 63°96 S2 
BueHeO25 10 26sl2)025 11) li.7 ~6.9 8.9 40.93 97.8 
36 00220 "30 923 000 2)? 22.0 6.8 8.5 21.01 95.4 9.4 <4 
Sees lo Hos) eeslOn lH 3/25" fel £829 10,95 91> 108) xe 8.5) <1 
37) 20871225 916 * 206 Ck TS COM ONs ey: 79342 163.'197 58 
a/eeoGEl 9320821. -066. 1602070". -7,3. 9.2 20.12 9523 18.8) <3 
3723 320 TOMO =2T@ 3hpo20 1029 ' PPPP Pe, Perorre S96 636 <2 
57235202. 10 26.7 
37-3 074°15 30 9 
Biees0C, OOO 1225 307,09 2002 9.4971.97911.78 91.9 AERO Ue 
S720) 15 25°lO 207 To Stoney Oo roel sero OTe 
37-4 056 20 23 8 
38 076"24745" 27" "000 26°24. 49 "Oss" tS? 9.88 90.3 ° 3.3595 ° 5.8 <2 
Bem so M1 15 20) e163 10.2) 10.2 13227 10.60 691.2 6.0 <2 
39° 357 00 35 19 “l50rTOl t7e9 8.7 11.0 15.92 94.0 
BO mEOH IO 35220 8.040 09 4Ntd 245 6927 12.3 12.21) 92.2 Side SL 
ho O45 11 35 30 O42 12 28.1 7 OPO SS 21551619356 ONG Zaks 68 3NK1L 
POO 100s 45 Ou O50 1529 @12.0516.4 16.16 86.4 
Ei 019" 60" 45°27" "024°01 “2450 ""1090 12.6°*8476"89!0 P86 “cls 
Dees ede 3Oe ws 350r 155469" eel e> Oe Liesl Gliese 166-95 976.9. <1 
BEM 29607 30820 9295808 20.3" 7.1 950: 15.66 93.6 6.6 <2 


Ait ae) 
4 ae A 


8 
| rae o/h ae 5 i 
ae th ee Lea 
fi - we ; a Ah ba) a a py 


- "Ant 
ae, .! 


TTT IO eat 
JAIXA a wen) 
eee seh neste ‘- 


ee +h P* 
eT 23 $4 “23 (aa I Kee 
aia ti baat i ey ee st ees 
FZ aan ba a Of: st i 12. * tes 
he £88 saat pet ® i 7 se “fo ess 00 O¢0 
ee oo) ie ee ate St OE BF age 
«Bae ERS uit 2,01 aed: | J 
( 8.8 e.1e 2e.0F #8) 8.e fos et Ser Ss 
| nv wes) nw om @ ose 
28> 1.2 1.98 8.8 esre.e a8 tI te’ ae 
aC EE 28> ef «OMe TESE D8 BY 7.63 Co MSo OS Oe 
$,08 €8.8S O18 Ff.) é. t} fo 280 Bt EH 
Bye €e.08 2.8 8.5. V.TF Tt eso St oe 
> de  § ae 10.78 28 8.8 0.88 1S OO ES 0 a4 
> a8 o> HOt f.42 CGP e.8- 1.0 3 are af Of€ if OF BT: 
..% 9.x Ea. SH EN N.2., (Oger of dos ot eS st 
> S.tt ee ed ek a ‘BY dd0 TS OF e 330. 
a> 3.3 3.08 HF.@ Tet PTT. os TIE 
j ih F Re . ; * 3: ot soe en 
gis> id e,fe BX.ff e.1i Ae 8.08 20 Soe ss ob 8 Boe aa 
T.¥e B.C 1.01 BY Be et TOS OF bis 
een Avi. be 2 en eee ee Th iy ee .- 2 “eS 
O> 8.2 88 GE. ECR Re Wert t.e Ais “as 80 «TS Ld 
gS? 0.9 | $46 OBO! Sef A.OT S$.81 df SAE OS ea 8e 
cake gfe 80.27 Hf N.8 8 @.01 fo Get Gl BE 00.7! 
a> 0.8 S90 tS.S1 BST Te 4.61 20 odo OS 2€ O 
i> €.8 af Gud S.8eer.27°S8 Oy | 1.82 Sk SHO Of zt | 
PBR STB AST et eat fo gor Bi | 
ee CF §.88 a7.8 aeSf O.Of 0.AS JO #g0 YS & 
> e.d a at Efe Sil ee OX eS El ee 
WS) O05 Bvge de,2f O.2 1.5 L& as: ms a | 
5 Y ' 7 | ; at i, ; 
; ry if : 
i 


84 
Summary and Conclusions 


Andrews and Shimizu's (1966) suggestion that the treatment of till 
fabric data could benefit from use of a three-dimensional statistical 
procedure based on work by Fisher (1953) and others concerning probabil- 
ity distributions on a sphere was a good suggestion. They recognized 
the problem posed by the bimodal nature of most till pebble a-axis dis- 
tributions, and proposed the solution of deleting transverse pebbles 
before treating the data. 

The procedure they recommended, however, is not suitable for treat- 
ing a-axis orientations of till pebbles for the following reasons: 

1 The condition k 23 is used as the sole criterion of applicability 
of the method, that is, as indicating that the sample is essen- 
tially unimodal. This is not true. 

2 The recommended procedure for deleting transverse pebbles fails. 

3 Watson's (1956b) test for randomness of directions based on the 
length of the resultant vector is not valid for hemispherical 
distributions. 

At least one subsequent publication depended heavily on the statis- 
tics obtained by Andrews and Shimizu's recommended procedure for repre- 
sentation of till fabric data. Since it has been shown that such 
statistics are quite unreliable and in many cases physically meaningless, 
the actual nature of the data used in that study was not communicated. 
Furthermore, since the conclusions of the study were in large part based 
on the statistics, those conclusions must be regarded as questionable. 

It is important to recognize that the method fails not because it 


is a bad or unsound method but because it is used in a situation for 
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which it was not designed. Any statistical procedure is based on 
certain assumptions about the data being treated; if these assumptions 
are not valid then the method is not sound. In this case the under- 
lying assumption is that the data conform to a Fisher distribution. 
Hence Andrews and King's application of the procedure to bimodal dis- 
tributions (Andrews and King, 1968, p.446) is not valid, contrary to 
their assertion that "This result may be surprising but is statisti- 
cally sound.!'' (Andrews and King, 1968, p.457). The importance of 
constructing valid models for statistical studies of geological orienta- 
tion data was emphasized by Pincus as early as 1953. 

An alternative approach is suggested based on the assumption that 
a multimodal axis distribution can be adequately described by a model 
consisting of several Fisher distributions, one corresponding to each 
mode. A graphical procedure has been described for setting up such a 
model for a sample of till stone long axes. Testing of till fabric data, 
however, has shown that the assumption that a Fisher model can adequately 
describe even a single mode of a till fabric sample is not justified. 
This may be due to the elongate nature of most till fabric modes, result- 
ing from the preference for a near-horizontal plane observed in the 
majority of samples. Thus, statistical comparison of samples using the 
parameters calculated by this method cannot be regarded as rigorous. 

The calculation of descriptive statistics for till fabric samples 
based on valid three-dimensional models must await further work on 
either (1) models for elongate groupings or (2) procedures permitting 
a grouping to be treated separately from a girdle occurring in the same 
sample. 


Until adequate techniques for the numerical representation of till 
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fabrics are devised and widely accepted, till fabric data should be 
presented graphically. This ensures that a publication will not become 
entirely useless should the numerical procedures used in it be shown 

to be invalid. Computer programs are now available that reduce drasti- 
cally the man-hours involved in the production of stereographic diagrams, 


as well as eliminating inaccuracy and subjectivity. 


Possibilities for Future Work 


Work that needs to be done on the numerical treatment of distri- 
butions of long axes of till stones may be divided into two kinds: (1) 
new procedures based on presently available models, and (2) establish- 


ment of new models. 


New Procedures with Present Models 


Watson (1956b) derived the distribution of R, the length of the 
resultant vector, for a set of unit vectors distributed randomly on the 
sphere, and determined some of the percentage points of this distribution. 
For the treatment of non-directed axes it would be very useful to know 
the distribution of R for a set of unit vectors distributed randomly on 
a hemisphere, and some upper and lower percentage points of this distri- 
bution. Observed values of R (or vector magnitude, R/N) could then be 
compared with the known percentage points and be assessed as significant 
or non-significant at a suitable confidence level. Such an ''R-test'! 
for a hemispherical distribution would actually be more versatile than 


the same test applied to a spherical distribution. On the hemisphere, 
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high R values would indicate a tendency for the observations to group 
about one direction, and low R values would indicate a tendency for the 
observations to avoid a certain direction or be concentrated into the 
the plane normal to that direction. Random distributions would have 
intermediate R values. In the case of a spherical distribution, high 

R values indicate a tendency for the observations to prefer one side of 
the sphere; low R values, however, do not indicate deviation from random- 
ness, since low R values are expected for random distributions on the 
sphere. Since low R values are also expected from distributions having 
a preferred plane, the ''R-test'' on the sphere is not suitable for de- 
tecting this form of non-randomness. 

In the treatment of paleomagnetic vectors a preferred direction 
is the only form of non-randomness that is of interest; hence R provides 
a convenient test. In the case of till stone long axes, however, a 
preferred plane is common and must be considered as a relevant alterna- 
tive to randomness. Thus the distribution of R for random sets of 
vectors on a hemisphere would be particularly valuable for till fabric 
work since (1) pebble axes, being non-directed, must be treated in a 
hemisphere if treated vectorially, and (2) the value of R in a hemi- 
spherical distribution is useful for detecting non-randomness of the 
forms most frequently found in till fabrics. 

Since the hemispherical distribution of a set of non-directed 
axes depends on the position of the hemisphere, a suitable procedure 
for using the value of R as an indicator of non-randomness would be as 
follows: (1) Using many different hemispheres, calculate many values 
of R. (2) Determine the maximum and minimum values of R. (3) Compare 


these with known significance points of R. 
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The vector magnitude program already used performs step (1) 
(printing R as a percentage of N, the sample size), and permits a 
rapid graphical solution to step (2). Thus if significace points 
of R (and therefore of R/N) were known, the vector magnitude diagrams 
already produced and described elsewhere in this report could become 
a rapid and simple means of evaluating the probability of non-randomness 
of the two forms most common in till fabrics -- a preferred direction 
and a preferred plane. 

A complete procedure for evaluating a till fabric sample might be 
based on this method as illustrated in Figure 8. The existence of a 
preferred plane and/or a preferred direction would be established by 
comparing the maximum or minimum value of R determined from the diagram 
with a suitable percentage point of R. It may be that a strongly pre- 
ferred plane containing a random distribution of azimuths could contain 
random groupings that would produce a significantly high value of R. 
(In other words, a significantly high R value would indicate non- 
randomness, but not the kind of non-randomness.) For this reason an 
indication of a preferred plane (that is, a low R value) should be 
investigated before a significantly high R value is interpreted as 
indicating a preferred direction. Only if there is no indication of 
a preferred plane can a high R value be taken directly as indicating a 
preferred direction. 

If no significantly low R value is found, thena significantly high 
value is looked for (fig. 8). If none is found, the distribution is re- 
garded as random. If one is found, it is concluded that there is a 
preferred direction and the distribution may be represented by a single 


Fisher distribution. The appropriate parameters are calculated. (Only 
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R diagrams with a single peak are considered at this point in the 
procedure. If the R diagram has two peaks a significantly low value 

of R will probably have been found, since there is therefore a pre- 
ferred plane. Three or more preferred directions that are far from being 
coplanar are not regarded as an alternative that needs to be considered 
when dealing with till fabrics.) 

If a preferred plane is indicated, the next step is to examine the 
distribution of azimuths in that plane and look for any significant de- 
viation from randomness. A procedure should be used that overcomes the 
difficulty presented by the possibility of the observed distribution 
being non-random but having two-fold or four-fold symmetry in the pre- 
ferred plane. Such a procedure might be similar to that employed in the 
mean vector program (Appendix D). 

If no significant deviation from randomness can be found, the 
distribution is regarded as consisting of a uniform girdle, and, if 
desired, parameters of the distribution may be calculated using one of 
Selby's (1964) models. If the azimuth distribution is found to differ 
Significantly from randomness, the next step is to determine the pre- 
ferred direction or directions. Ultimately it may be possible at this 
stage to separate the pattern into its linear and planar components 
mathematically and calculate parameters for each component (see below). 
At present, however, one must be content with subjective assessment 
based on visual examination of either the three-dimensional density 


diagram or a diagram of the azimuth distribution in the preferred plane. 


Both Fisher's (1953) model of a grouping about a line and Selby's 


(1964) model of a girdle distribution are useful models when they are 
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applicable. However, one of the main problems of dealing with till 
fabrics is that the pattern is not only complex but highly variable. 

If a fabric pattern is considered as being composed of a combination of 
planar and linear components (girdles and preferred directions), till 
fabric patterns commonly have a girdle and two preferred directions, a 
girdle and one preferred direction, or a girdle with no preferred direc- 
tion. Samples are considered in this study that have no girdle but only 
a preferred direction (for example, samples 22 and 26 (fig. 9) and 466-2 
and 466-4 (fig. 10)). Fisher's model is thus applicable in these latter 
cases, but for the majority of till fabric samples neither Fisher's nor 
Selby's distribution alone is adequate to describe the sample. Since 
Fisher's and Selby's distributions describe the two basic components of 
till fabric patterns -- the girdle and the grouping about a line -- one 
is led to hope that it will one day be possible, probably using a com- 
puter, to fit a complex model consisting of a combination of the two 
basic forms to an observed set of orientations. In the meantime, the 
first step would be to devise procedures (again taking advantage of the 
computer) to simulate random samples of limited size from hypothetical 
populations that fit the basic models. This would perhaps lead to pro- 
cedures for evaluating bimodal and multimodal samples to establish 
whether two adjacent modes in a sample should be regarded as represent- 
ing one or two modes in the population. At present this is a matter of 
guesswork. Simulation could then be carried out using combinations of 
the basic models, so that a reasonable complex model could be constructed 
to fit a complex observed distribution, using repeated approximations 
and visual comparison. 


One advantage of such a composite model would be that a single 
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component of the fabric could be described separately, and compared 
separately with single components from other samples, without its 
properties being masked or distrorted by the other components of the 
fabric. Relationships between different components of fabric patterns 


could then be studied, since each component could be treated independent - 


ly. 


Improved Models 


Many of the fabric patterns considered in this study consist of 
one or two near-horizontal groupings that are more or less elongated 
horizontally. This elongation is sufficient in most cases to constitute 
a significant deviation from a Fisher distribution, thus invalidating 
statistical procedures based on the assumption of a Fisher distribution. 
It was suggested above that these elongate groupings might be represent- 
ed by a model consisting of a Fisher distribution superimposed on a 
uniform girdle. The compound nature of such a model might, however, make 
the problem of estimating appropriate parameters very awkward. An alter- 
native approach might be to use as a model a spherical normal distribu- 
tion without axial symmetry (Pincus, 1953, p.504). If statistical 
comparisons could be set up on the basis of such a model similar to 
those presently in use based on a Fisher distribution, they could 
probably be applied to many till fabric samples that cannot be repre- 


sented by a Fisher distribution or a combination of Fisher distributions. 
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Test of Fit Taking Measurement Errors into Account 


Since the recorded orientation of an axis is not necessarily 
the true orientation of the axis, showing that recorded orientations 
do not have a Fisher distribution does not imply that the true orienta- 
tions do not have a Fisher distribution. Ideally, a test of fit should 
be devised that allows for the probable deviations of the observed 


orientations from the true orientations. 
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THE PROBLEM OF INTERPRETING 


PEBBLE FABRIC PATTERNS 


Reliability of Samples 


Fabric at Sample Site 


Introduction 


The problem in interpreting a fabric diagram is: What can be 
inferred about the fabric represented by the sample? This involves 
separating the random from the non-random elements in the sample, 
that is, characteristics peculiar to the particular sample from those 
of the population, in order to establish the presence or absence of 
a preferred orientation in the population and to determine the preferred 
orientation if there is one. There are two possible ways of tackling 
this problem: (1) Comparison of the sample with various models. This 
is done by setting up, either physically or mathematically, artificial 
populations with known distributions, then producing random samples from 
these artificial populations and comparing them with the real sample 
being evaluated. (2) Comparison of successive samples from the popula- 


tion being studied. 
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Comparison with Models 


One valuable aid is the simulation of random samples from model 
distributions of various kinds. This gives an idea of the degree to 
which the form of the distribution is distorted by the random element 
of the sample for various sample sizes, and thus provides a basis for 
interpreting a given sample. The following models are useful in dealing 
with a-axis orientations of till pebbles: 

1 Axes uniformly distributed. 

2 Axes concentrated in or near a plane, with all parts of the plane 
occupied with equal frequency (uniform girdle). 

3 Axes concentrated in or near a plane, and in that plane near one 
particular orientation. 

4 Axes concentrated in or near a plane, and in that plane near each 
of two orientations at right angles. 

5 Axes concentrated around one orientation, with no plane preferred. 

A type-1 distribution and a crude model of a type-2 distribution 
may be easily simulated. Distributions of types 3, 4 and 5, however, 
are difficult to simulate, and are best tested for by comparison of 
successive samples. Since distributions of types 1 and 2 are the only 
ones in the above list that have no preferred orientation, comparison 
with them can be used to establish that an unknown population does have 
a preferred orientation. Since one of the problems in the present study 
is the detection of a preferred orientation for the pebble axis popula- 
tions sampled, comparison with type-1 and type-2 distributions is a 
useful technique. Simulated samples from such distributions were pro- 


duced for comparison with the real fabric samples, and are illustrated 
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in Figure 5. 

In comparison with uniformity over the hemisphere, observed 
values on a point density diagram may be assigned a definite probability 
using the binomial distribution or the Poisson approximation to it. It 
is felt that such a technique is of limited usefulness, however, and 
that a qualitative visual comparison of the point density diagrams is 
just as effective and more versatile. The procedure described by Flinn 
(1958) was designed for testing samples against the model of uniformity 
over the hemisphere, but it requires the use of a particular method of 
contouring not used here. 

Random samples from a uniform distribution over the hemisphere 
were obtained in the following way. A 10-inch rectangular grid was 
superimposed on a 10-inch diameter equal-area projection. Points were 
plotted on the projection using two-digit rectangular coordinates ob- 
tained using a table of random numbers. Points falling outside the 
circumference of the projection were ignored, and the procedure was 
continued until 50 points lay within the circle of the projection. The 
trends and plunges rauidedhwad by the points were then read off using 
an equal-area or Schmidt net, and a point density diagram prepared in 
the same way as for the fabric samples. 

The population sampled consists of the points defined by all 
possible pairs of two-digit coordinates, that is, all possible sequences 
of four digits. Since these are the non-negative integers from 0 to 
9999 inclusive, each has the same probability of being obtained from the 
random number table. Hence the population consists of 10,000 points 
distributed uniformly over the grid, and each has the same probability 


of being selected for the sample. Since the projection is equi-areal 
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the points within the projection represent points distributed uniformly 
over the reference hemisphere. In this way, samples of 50 random orienta- 
tions were obtained. They are illustrated in diagrams 505, 506 and 507 
(fig. 5). 

A simple model of a distribution of type 2 consists of a uniform 
distribution over that part of the sphere less than a specified angular 
distance from a given plane. Ina real girdle distribution, the proba- 
bility density would be a maximum in the plane, decreasing away from the 
plane. Since this is not easily simulated, the uniform distribution 
was chosen, using a limit of 30 degrees on either side of the plane. 

This limit was drawn on the projection (using the plane of the projection 
as the plane of the girdle) and the above-described procedure for obtain- 
ing randomly distributed points was followed until 50 points lay within 
the annular area representing the girdle. Five such samples were obtain- 
ed; they are illustrated in diagrams 508-512 (fig. 5). 

A x2 test was used to check the randomness of the above samples, 
from both distributions. For the random distributions over the entire 
hemisphere both the azimuthal distribution and the polar distribution 
(Frequency in concentric equal areas) were tested. For the girdle distri- 
butions, only the azimuthal distribution was tested. Two of the "i 
values were Significant at the five per cent level: the polar distri- 
bution in sample 506 has a x? value with a probability of less than 
one per cent, and the azimuthal distribution of sample 512 has a ye 
value with a probability of slightly greater than one per cent. 

Samples of 100 points from the girdle were obtained by combining 


sample 508 with sample 509, and sample 510 with sample 511. These are 


labelled 513 and 514 respectively (fig. 5). Clearly, field samples of 
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100 axes must be compared with test samples of 100 orientations for the 
comparison to be valid. All 250 points from the uniform girdle were then 


run as one sample; this diagram is labelled 515 in Figure 5. 


Comparison of Successive Samples 


An alternative method of evaluating the random and non-random 
elements is the selection of several independent samples from the un- 
known distribution and comparison of these with one another. The non- 
ramdom element, that is, the density distribution of the population, is 
common to all the samples, whereas the superimposed random ''noise!'' will 
be different for each sample. Comparison of several independent samples 
enables the common features to be singled out. This may be called a test 
of "'reproducibility'' of the samples. It is of particular value when the 
population under study departs only slightly from uniformity, so that 
Comparison with artificial samples from uniform populations reveals no 
great deviation. In such a case, the reproducibility of the slight non- 
randomness may be regarded as significant. 

This technique is of particular value also when no assumptions can 
be made about the probable distribution in the population, or when the 
models that would provide the most useful comparisons cannot be easily 
simulated. For example, comparison with random samples from a model of 
a uniform girdle (that is, one having a uniform azimuth distribution in 
the plane of the girdle) may indicate that the population under study 
has a non-uniform girdle. Further tests of this sort would require use 
of models of different kinds of non-uniform girdle; such distributions, 


however, are not easily simulated. In this case the form of the distri- 
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bution must be established by comparison of successive samples. 

Some valuable information about the adequacy of the 50-axis 
samples used in this study is provided by the duplicate samples taken 
at sites 34, 35 and 37 at location F (fig. 11). Each of these samples 
was initially plotted as two groups of 50, using the first 50 measure- 
ments and the second 50 measurements, to test the reproducibility of a 
50-axis sample. These plots are shown in Figure 11 as subgroups 1 and 
2. As pointed out by Flinn (1958, p.526) this test of the reproducibil- 
ity of a sample requires the knowledge or the assumption that the fabric, 
over the distance from which the sample was taken, is homogeneous; only 
then can differences between successive groups of observations be inter- 
preted as due to sampling rather than spatial variation of the fabric. 
This difficulty can, however, be overcome by a procedure employed by 
Kauranne (1960). This involves placing all the even-numbered observa- 
tions in one group and all the odd-numbered observations in a second 
group, then comparing these two groups. Any spatial change of the 
fabric will thus equally affect both groups, and the two groups are 
truly duplicate samples of the same fabric. Hence any differences can 
be interpreted as being due solely to sampling. The groups obtained in 
this way are shown in Figure 11 as subgroups 3 and 4. 

The differences between subgroups 1 and 2 are not noticeably 
greater than the differences between subgroups 3 and 4 at either of 
sites 34 or 35, suggesting that the fabric at these sites was actually 
homogeneous throughout the volume of till from which each sample was 
taken. The difference between subgroups 1 and 2 at site 37 does appear 
to be somewhat greater than the difference between subgroups 3 and 4, 


suggesting that the fabric was not in fact homogeneous throughout the 
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volume of the diamicton from which the measured pebbles were taken. This 
is not altogether surprising, since it supports the stratigraphic assign- 
ment of site 37 to the basal diamictons of the Lake Edmonton sediments 
instead of to the upper till. 

The reproducibility of the samples will be evaluated by comparing 
subgroups 3 and 4 for each site. At site 34 the maximum that appears at 
024 10 in subgroup 3 and at 023 10 in subgroup 4 is clearly reproducible 
with fairly high accuracy. The two smaller maxima lying on either side 
of the principal maximum and having somewhat larger plunges are also re- 
producible but with lower accuracy. 

Samples 35-3 and 35-4 both indicate a preference for north-north- 
easterly and south-southwesterly trends. The grouping here is evidently 
not as tight as that at site 34, there being no prominent reproducible 
maximum in the samples of 50 axes. Only the general preference for 
north-northeasterly and south-southwesterly trends is reproducible. 

At site 37, the reproducible features are the three groupings 
trending northwest, east-northeast and south-southwest. The positions 
of the density peaks are not constant. 

The above experiment shows that a sample of 50 axes is usually 
sufficient only to indicate approximately the preferred orientation 
or orientations. Only very high densities in the sample can be re- 
garded as accurate indications of a preferred orientation in the pop- 
ulation. It is interesting to note that while the peak at 023 10 in 
sample 34-4 and the peak at 207 12 in sample 35-3 are comparable in 
both extent and magnitude, the former is accurately reproduced in the 
duplicate sample whereas the latter is not. In both cases the maximum 


density exceeds 15 per cent in a three per cent area. Densities of this 
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magnitude appear in three of the five samples taken randomly from a 
uniform girdle. 

Thus a density of this magnitude may be closely reproducible or 
not reproducible at all; only additional observations can determine this. 
Hence, assuming the criteria of pebble selection used in this study, 50 


measurements are barely adequate to define the pebble fabric of a till. 


Summary 


In view of the preceding discussions, interpretation of the fabric 
patterns represented by the 50-axis samples discussed below will be done 
using the following procedure. 

First the three per cent density diagram of the sample is compared 
with the test samples from a uniform distribution and a uniform girdle. 
If it appears that the fabric is not uniform and does not consist simply 
of a uniform girdle, then the observations are assumed to be grouped 
about one or more axes. The problem then is to determine the number of 
preferred axes; comparison with the duplicate samples from sites 34, 


35 and 37 is found to be useful in this step. 


Fabric at Exposure 


Representation of the fabric of a till unit at an exposure by a 
single one-point sample is valid only if the fabric is homogeneous 
throughout the unit at the exposure. Homogeneity of the fabric of a 
till unit at any exposure implies that the probability of finding a 


pebble with any given orientation is the same at every point in that 
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unit at the exposure. An attempt was made to test the homogeneity of the 
fabric of the upper till at location F (fig. 12). The till is exposed 

on a slope with a northward dip of about 18 degrees that was made in pre- 
paration for mining of the underlying Saskatchewan Gravels. There was 

no visible evidence of any mass downslope movement. 

Fabric samples were taken from the upper till at the nine sites 
indicated in Figure 12. The samples were grouped both horizontally and 
vertically in order to detect any systematic variation in the fabric 
either vertically or laterally. Clearly the patterns shown by the sam- 
ples differ from place to place. The problem is to determine to what 
extent these differences are due to ''random noise'' and to what extent 
they are due to differences in the fabric. 

The fabric patterns are interpreted using the principles estab- 
lished in the above discussion. Comparison of sample 40 with the test 
samples (fig. 5) shows two things: (1) the observed axes are concen- 
trated close to a near-horizontal plane, and (2) in this plane, the 
direction 045 14 is definitely preferred; the grouping about this 
direction is so tight as to suggest that the two smaller groups at 
100 11 and 150 00 are not merely part of the ''tail'' of this grouping 
but in fact indicate one or more other elements in the fabric pattern. 
Whether they are simply accidental groupings in the preferred near- 
horizontal plane, or whether they indicate two more preferred orienta- 
tions Cannot be established. The maximum dip of the preferred plane, 
as determined from the vector magnitude diagram, is 14 degrees toward 
037 degrees (Table 2). 

The axes in sample 38 are concentrated close to a plane dipping 


29 degrees towards 039 degrees (Table 2). Comparison with test samples 
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508-512 indicates no great deviation from the patterns expected from a 
uniform girdle. Sample 42 is similar to sample 38 in being concentrated 
close to a plane but having no reliable indication of a preferred 
direction. The preferred plane dips 29 degrees towards 011 degrees. 

The axes of sample 6 are concentrated close to a near-horizontal 
plane (the position of the minimum on the vector magnitude diagram in- 
dicates that this plane dips 5 degrees towards 125 degrees). In this 
plane there is clearly a preference for the position indicated by the 
density maximum at 066 00. 

Comparison of sample 34 with test samples 513 and 514 (fig. 5) 
indicates a preference for a near-horizontal plane (dipping 07 degrees 
towards 055 degrees) and a clear preference for northeasterly and south- 
westerly trends. As discussed previously, the constancy of the position 
of the density maximum at 024 10 in independent samples of 50 axes from 
this site strongly suggests that this is a preferred orientation of 
the till pebbles at this site. The sampt#e suggests that there is a 
second preferred orientation near the density maximum at 099 05, but 
its position is indicated only approximately. 

The axes of sample 35 are clearly, though loosely, grouped about 
a near-horizontal line trending about 026 degrees. The vector magnitude 
diagram indicates that the preferred plane dips 53 degrees towards 295 
degrees, but this plane is only weakly preferred. 

Sample 36 shows a preferred plane that dips 29 degrees towards 
038 degrees. The double-peaked grouping in the northern part of the 
diagram probably represents a preferred direction. 

Sample 39 shows a preferred plane dipping 14 degrees towards 079 


degrees. The gaps in the girdle at 080 and 118 degrees and the low at 
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006 degrees suggest a sparsity of pebble axes in the population with 
north-south or east-west trends, that is, a preference for northeast- 
southwest and northwest-southeast trends. 

Sample 41 also has a preferred plane, dipping 05 degrees due 
north. In this case the two gaps in the girdle at 060 and 155 degrees 
suggest preferences for the east-southeasterly and north-northeasterly 
parts of the girdle. 

Thus the nine samples indicate preferred orientations ranging 
from north to east-northeast. If the degree of reproducibility demon- 
strated for sites 34 and 35 (fig. 11) is assumed for all nine sites, 
then a substantial variability of the fabric must be inferred. The 
implication. of this conclusion is that a one-point sample from the 
upper till at Location F is likely to give misleading information if 
it is used to infer ice movement direction. Sampling procedures designed 
to overcome this difficulty are discussed below. It is interesting to 
note also that there appears to be no greater consistency of the fabric 
horizontally than vertically. 

The approach of using a single one-point sample to represent the 
fabric of a till at an exposure may be further evaluated by considering 
other locations at which several samples were taken from the same unit. 
Reasons for the differences or similarities of the samples are not 
discussed; it is only the fact of the variability that is of interest 
here. Discussion of fabric origins will be found in a later chapter. 

Eleven samples of 50 axes are available from the till exposed at 
Location B (fig. 10). Considerable variability is evident in both units. 
Even the three samples from the top of the till include both northwest- 


southeast and northeast-southwest preferred orientations. Thus the same 
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high variability is present laterally as vertically. 

Samples 4 and 5 (fig. 6) were both taken from the lower till at 
Location E, about 200 feet apart. The lower till is about 12 feet thick; 
sample 4 is from the top two feet, and sample 5 from three to four feet 
above the base. The two samples bear no resemblance whatever. 

Samples 7 and 21 (fig. 6) from the base of the lower till at 
Location F were taken about 300 feet apart and are quite similar, both 
showing a preference for a northeast-southwest trend. 

At Location G (fig. 7) the writer originally obtained samples 
31 and 32, respectively from the lowermost two feet and the upper- 
most two feet of the 14 feet of upper till, sample 32 being vertically 
above sample 31. No similarity is evident. Just over a year later, 
sample 466-5 was obtained by another worker from the same gravel pit at 
about 10 feet below surface. The exact site is almost certainly at 
least 100 feet from sites 31 and 32. Assuming seven feet of Lake Edmon- 
ton sediments as mapped by the writer, site 466-5 would be about three 
feet below the top of the upper till, and thus stratigraphically close 
to site 32. The similarity of samples 32 and 466-5 is striking. 

Samples 8 and 9 (fig. 7) from Location H were taken from two 
different ''phases'' of the lower till that are easily distinguishable 
in the field (Plate 2-G). Phase A (sample 9) is dark grey and has a 
fractured, blocky structure. Phase B (sample 8) is brown and has a 
characteristic granular structure. Samples 8 and 9 are clearly very 
different. However, no information is available concerning the vari- 
ability of the fabric within each ''phase'' of the lower till at this 
exposure. 


At Location I (fig. 9) the lower ''till'' consists of three members 
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separated by bedded sands, silts and clays. The entire lower till 
complex is 23 feet thick. Sample 18 is from the middie member (4 

feet thick); sample 26 is from the lower member (6 feet thick) directly 
below sample 18; sample 44 is from the middle member about 10 feet 
horizontally from sample 18. The similarity of all three of the above 
Samples is remarkable. 

Samples 22 and 23 (fig. 9) from the lower till at Location K 
at first sight appear dissimilar, but actually are alike in having 
a fairly high proportion of observations trending east-northeast or 
west-southwest. 

Although differing in detail, samples 466-3, 466-8, and 466-11 
(fig. 13) from the upper till at Location L all indicate a general 
preference for a northeast-southwest trend. Sample 466-1 differs in 
having prominent maxima in both the northeast and the north-northwest. 
For the lower till at the same location, samples 45 and 466-10 (fig. 13) 
both indicate a preferred west-northwest trend. Sample 43, however, 
differs greatly from 45 and 466-10. 

Thus it can be seen that, although different samples from a till 
unit at a given exposure are in some cases very similar, this can by 
no means be counted upon, since frequently great variation is evident. 
It must be concluded that a single one-point sample should never be 
depended upon to indicate the preferred trend of elongate pebbles in 
a till unit at an exposure . Indeed, the assumption that the pebbles 
in a till unit at an exposure have a single preferred trend must be 
questioned. They may have many preferred trends, or they may have 
different preferred trends in different parts of the exposure, both 


laterally and vertically. 
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Effect of Pebble Size 


and Shape 


Possible effects of pebble size and shape on the orientation of 
the a-axis were looked for by measuring the lengths of the a-, b- and 
c-axes of 226 pebbles from sites 34, 35 and 36 in the upper till at 
Location F. These dimensions, together with the a:b and b:c ratios 
and the attitude of the a-axis are given in Appendix C. In addition, 
the dimensions and attitudes of a number of flat pebbles with low 
elongation were measured. Appendix C also gives the shape classifi- 
cation of the pebble according to the system of Holmes (1941). 

Two variables were examined for a relationship with the a-axis 
orientation: the a:b ratio (here called ''elongation'') and the a-axis 
length. Pebbles were grouped according to their azb ratio and point 
density diagrams prepared for the a-axes of each group (fig. 14). The 
top row of diagrams in Figure 14 shows the effect of raising the lower 
limit of the azb ratio. As the lower limit is raised from 1.15 to 1.25 
the maximum at about 024 14 becomes stronger, but no other changes in 
the a-axis distribution are evident. As the lower limit of a:b ratio 
is raised from 1.25 to 1.35 the maximum at about 024 14 becomes weaker, 
and the grouping centred at about 351 25 becomes stronger. These two 
trends continue as the lower limit of asb ratio is increased to 1.45 
and then to 1.55 at which point the maximum at 349 19 is stronger than 
that at 029 15. In addition, southwest-trending axes appear as a sepa- 
rate group centred at about 208 25, and the grouping that was evident 
at about 055 27 has disappeared completely. Thus there appears to. be 


a tendency for the north-northwesterly and south-southwesterly trending 
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maxima to be more strongly preferred by pebbles with higher a:b ratios. 
It is worth noting that the weak maximum at about 104 09 on diagram 106 
(a/b 2 1.15) survives right through to diagram 124 (a/b 2 1.55) and 
still trends 104 degrees. 

The bottom row of diagrams on Figure 14 illustrates groups with a:b 
ratios restricted to the ranges 1.15 to 1.25, 1.25 to 1.35, 1.35 to 1.45, 
and 1.45 to 1.55. If diagram 121 is compared with test samples 505-507, 
it appears to be random. On the other hand diagram 122 shows a preferred 
orientation of 017 17. This suggests that 1.3 would be a better lower 
limit of a:b ratio than 1.2 for use in future studies. However, the 
question of a:b ratio clearly requires further study, since it appears 
that the a:b ratio can affect the orientation of the pebble. A possible 
explanation for the relationship between a:b ratio and a-axis orientation 
observed here is discussed below in the section dealing with the move- 
ment of the glacier that deposited the upper till. Johansson (1968) 
also classified according to a:b ratio pebbles from till deposited by a 
glacier whose movement direction had changed. However, his data indicate 
no apparent correlation between a:b ratio and a-axis orientation. 

Figure 15A illustrates the effect of a-axis length on a-axis orien- 
tation. Three lower limits of a-axis length, 5mm, 25mm, and 40mm, were 
used, and these subgroups studied using two different lower limits of 
a:b ratio, 1.15 and 1.45. The top row of diagrams in Figure 15A illus- 
trates the effect of an increasing lower limit of a-axis length using 
all pebbles with a:b 2.1.15. As the lower limit of a increases, the 
maxima become stronger (except the small one at 104 09), but no change 
in the relative magnitudes of the maxima is evident. The same effect is 


seen in the bottom row where a:b 27 1.45; here the strengthening of the 
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maxima iS even more pronounced. Again, no substantial change in rela- 
tive strengths of the maxima is seen. Thus, while there is no. apparent 
relationship between the size of a pebble and the orientation preferred 
by that pebble, larger pebbles do display a far smaller amount of 
scatter than do smaller ones. Clearly, the sizes of pebbles measured 
for a fabric sample will therefore partially determine the minimum 
number of observations required, since the more scatter there is in the 
observed orientations the more observations are necessary. Johansson 
(1968) classified till pebbles according to size. His data indicate 
that larger pebbles have more scatter in their orientations than smaller 
pebbles. This difference may be due to the fact that the pebbles 
measured by Johansson were large compared to those used in the present 
study and in higher concentration, so that considerable interaction of 
pebbles probably took place during glacier movement. 

Figure 15B shows the distribution of the c-axes (or poles of the 
a:b planes) of 84 '"'flat'' pebbles from the same three sites -- 34, 35, 
and 36 -- in the upper till at Location F. The highest density is very 
close to the centre of the projection, indicating that the most strongly 
preferred attitude of pebbles with high b:c ratios (thateis;, 'f tax" 
pebbles) is with the asb plane nearly horizontal. A similar observation 
was made by Harrison (1957) and Johansson (1968). Diagram 128 in Figure 
15B shows a second grouping of c-axes with plunges from 18 to 34 degrees 
and trends from about 090 to about 125 degrees. These pebbles have 
their a:b planes dipping steeply to the northwest and striking from 
northeast to north. A possible explanation for this grouping of c-axes 
will be discussed below in relation to the movement directions of the 


glacier that deposited the upper till. 
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Recommended Sampling Procedures 


The above data illustrating the variability of fabric pattern over 
short distances within a single stratigraphic unit clearly indicate that 
the approach of using a one-point sample to characterize the fabric 
pattern of a unit at an exposure is not valid. The data support the 
similar conclusions of Kauranne (1960), Andrews and Smith (1966), and 
Andrews and King (1968). 

The question thus arises of what valid approach can be taken to the 
problem of measuring till fabrics. Clearly, if single-point samples are 
inadequate, an approach involving multiple-point sampling must be adopted; 
that is, many small samples from different parts of an exposure must 
replace a single large sample from one point in the exposure. Andrews 
and King (1968) suggested such multiple-point sampling and the use of 
Watson and Irving's (1957) procedure involving calculation of within- 
and between-site precision parameters and to determine the optimum 
sampling pattern. As was pointed out earlier, this method is not appli- 
cable to till fabric data since the latter cannot be described adequately 
by Fisher-type models. The problem of how many points should be sampled 
and how many pebbles measured at each site must be solved by some other 
means. 

To this end, some experiments in sampling were performed using the 
data available from the upper till at Location F. A total of 550 long 
axeS were measured at nine sites (fig. 12). Although substantial vari- 
ability of the fabric pattern was evident from site to site, the density 


diagram for the whole group (number 118, fig. 12) showed a single maximum 
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at 029 10, coinciding with the trend of groove molds (Plate 1-A, 2-E) 
from 028 to 033 degrees seen on the base of the upper till at several 
localities within one mile of Location F. This coincidence suggests 
that the whole group does indicate the dominant ice-movement direction 
in spite of the inhomogeneity of the fabric apparent from the individual 
samples. Although it cannot be assumed that this relationship holds at 
other localities and for other till units, it was used as a basis for 
experimenting with some possible sampling patterns. 

First, samples of 50 axes were made up by taking 10 from each of 
the five sites 39, 40, 35, 41, and 42 (fig. 12). The first 10 observa- 
tions from each site were used to make the first sample, the second 10 
observations from each site to make the second sample, and so on. The 
resulting samples are numbered 101-105 and appear in Figure 16. Sample 
101 shows no clear indication of non-randomness in its near-horizontal 
girdle when compared with test samples 508-512 (fig. 5). Samples 102- 
105 show a more or less definite preference for northeasterly and south- 
westerly trends, but the trend of the strongest maximum ranges from 028 
to 048 degrees. Thus these samples do not provide a reliable estimate 
of the maximum of the whole group at 029 10. 

For the next experiment, the sample size was increased from 50 to 
80 and all the available sites included. To make each experimental 
sample, 10 pebbles were taken from each site except 6 and 34, each of 
which was represented by only five observations since they lie close 
together. Again, five samples were constructed; they are shown as 
samples 131-135 on Figure 16. Only three of the five, 132, 133, and 135, 
indicate a preference for a northeasterly trending orientation. Number 


131 shows a girdle dipping slightly to the northeast, but could not be 
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interpreted as indicating any preferred orientation. Number 134, if 

seen alone, would be interpreted as indicating two preferred orientations, 
the stronger dipping slightly toward 350 degrees, the other dipping 
slightly toward 080 degrees. It would be reasonable to infer from this 
sample ice movement from 350 degrees, regarding the maximum at 080 
degrees as representing transversely oriented pebbles. Thus a single 
sample of 80 long axes taken from the nine sites used here cannot be 
depended upon to indicate a preference for a northeasterly trending 
orientation. 

Finally, the sample size was increased to 200, 25 observations 
being taken from each site except 6 and 34. 25 observations were taken 
from sites 6 and 34 together, the split being first 12 to 13, then 13 to 
12. The resulting samples are numbered 136 and 137 on Figure 16. Com- 
parison with test sample 515 (consisting of 250 axes; fig. 5) shows that 
both samples 136 and 137 indicate (1) a girdle dipping slightly to the 
northeast, (2) a preference for a northeasterly trending orientation, and 
(3) a high degree of scatter, the scatter being greatest in the plane of 
the girdle. The strongest maximum of sample 136 is at 042 22; that of 
sample 137 is at 040 12. These two samples are thus very similar; they 
are also similar to the total group (sample 118, fig. 12) except for the 
difference (11 to 13 degrees) in the trend of the maximum. 

It is concluded that a sample of less than about 200 long axes 
from the same nine sites could not be relied upon to indicate the north- 
easterly trending preferred orientation. Fabric and other evidence lead 
to the conclusion that the fabric of the upper till at Location F is 
complex, that is, it has been produced by ice movement in more than one 


direction (in this case ranging from northeast-southwest to north-north- 


sty 

TE HET todmul <nakyedneino bea 
.2nofisineire beristeq ows gniseatbni 26 | 19 | 

pniggib radio eft ,2enypeb 02€ brswot Qladgtte peserquerersi 

2tdt mor? ietat of sidénease;.od biuow 71 .eseypeb: dtonnl aii 


080 t6 mmixem edt pnibispey .2eexpsb OE mort inemevom 997 olqmez 


sipnie 6 eudT ,2elddaq beinsiao yisevevens12 pniinezeiqey 26 2earpab 
ed fJonns2 sven, beau. 2stte antn eit mor? nese? Cone wl: 08 to olan 
pnibnert ylseteseditvon 6 10? aonarsteig 6 sispibat bed moqu bebriageb 
alt not teine tre 

2 olqmez oft .yllontt 
nedet o1ew enottevieedo 2S Ag bos r sean oti2 does: mort neadst pated 
ot Ef not ,€f oF ST t2exF% pated tiga, ae: alwaritngod #8 bre & 293fe. mort 
.3t owt, ao Yer bos ath tha-iadleis 916 2elqmea patiiues) oAT.. St 
sent ewonz (2 .prt yeaxs O@S to pritetenes) Ot alqmee deat Aviw noe treg 


2nottevisedo 28,008 oF icant, 


ert ot yl tdeite eniqgih sibyip 6 {n) steatbhat VET bine det ‘apiqnae dtod 


brs ynofistneine pnibrions \inedaeendvon 6 ” ORO IATENG 6 «s) <tasartion 
to ensiq ont nt teotserp pnisd detisoe wie, 7994692 to saxpeb dgid.e (€) 
to tedt ;S8.S00 46 2i O&f stqmee to. muniieem. tepnnate ent .albiig ods 
yent tel imie yey aaa sis 2etqmee ows ees Sf O00 46, at V&l elqmse 
odd net, tqoone (St «p77 BIT elqmee) quere: aca eeiieed wslimte oats 916 
.momixem odd Yo bneys on! at (zaorpsb €1 of i) eonens ib 

 23%68 enol 00S swods nari enh ie ii biennial se 
-ltron elt ateotbni ot nogy baile. ed ton.b 
_ bee! sansbive redto, ne iis’ 


=f 


113 


west-south-southeast). This may be the cause of the high variability 

of the fabric over short distances that in turn results in the high 
scatter observed in the total group. This high degree of scatter is 
the. reason why large samples are necessary to define the preferred 
orientation when the sample is spread over the whole exposure. Thus 

the complexity of the fabric at this exposure might be cited as the 
underlying reason for the necessity of such large samples, and the 
exposure regarded as a special case. However, it can never be assumed 
that an unknown fabric is not such a complex fabric. Few till fabrics 
have been studied in sufficient detail to establish either their complex- 
ity or their variability. Therefore, the worst must be expected at each 
new exposure, and the complexity and variability of the fabric tested. 

A decision concerning the minimum adequate size of the total sample 
should then be made on the basis of this information. 

Rigorous methods for determining reliability of sample means and 
adequate sample sizes have been devised (Dennison, 1961) but are appli- 
cable only to unimodal distributions and are probably of little value 
for till fabrics. The method of Harris (1969) for determining ''!minimum 
Significant orientation count'' may be a useful technique. It involves 
the use of sequential sampling (Pincus, 1953, p.487) to estimate minimum 
satisfactory sample size. Pincus (1953) recommended use of a system of 
Sequential sampling in situations involving more than one concentration 
of orientations. 

One possible procedure involving sequential sampling at a new till 
exposure might be as follows. A sample of 100 axes is taken, measuring 
20 at each of five sites located, perhaps, at the horizontal and vertical 


extremes of the exposure with one in the centre. The trends of the first 
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50 axes (that is, the first 10 from each site) are plotted on prepared 
polar coordinate paper. Similar plots are made of the first 75 axes 
(the first 15 from each site) and of the whole group of 100 axes. By 
comparing these plots, the investigator can decide whether or not the 
features of the fabric pattern in which he is interested have been esta- 
blished with the desired degree of reliability. If not, additional sets 
of five axes from each site can be added to the whole sample until the 
investigator is satisfied that the relevant aspects of the fabric have 
been determined with sufficient accuracy. 

A procedure such as this ensures that sufficient measurements wil] 
be made in the field to determine the required information. In addition, 
valuable data are obtained concerning the degree of variability of the 
fabric. If sample sites are judiciously chosen, independent information 
can be obtained about horizontal and vertical variation. The number of 
sites should ideally be related to the variability of the fabric; if an 
operator began with five sites and found great variability, he might be 
wise to add four more sites to form a grid of nine. 

It is emphasized that the significance of such a multiple-point 
sample when the fabric is variable is by no means established. The 
coincidence of the trend of the density maximum of the total group and 
the trend of basal groove molds for the upper till at Location F cannot 


necessarily be extrapolated to other locations or other units. 


The data examined here indicate that pebbles with a:b ratios of 
less than about 1.3 are not good indicators of ice-movement direction; 
more elongate pebbles should therefore be used. In addition, the a:b 


ratio of each pebble should be recorded, since pebbles with different 
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a:b ratios may have different preferred orientations. 

The sizes of pebbles used for a fabric study should also be 
recorded, so that any correlation between size and preferred orienta- 
tion or degree of preferred orientation can be discovered. This 
information could save much unnecessary work and improve results in 


later stages of the investigation. 


Work Needed 


Extensive investigation of small-scale pebble fabric variability 
in tills is clearly necessary. Suitable sampling procedures for deter- 
mining the overall fabric of a till exposure must be established, and 
the meaning of such an overall fabric when the fabric is variable must 
be investigated. 

The a:b ratio of a pebble apparently has some effect on its 
orientation in a till. It was suggested above that more elongate 
pebbles achieve their ''stable'’ orientation more rapidly than less 
elongate pebbles, so that when the direction of movement of a glacier 
changes the more elongate pebbles are the first to be re-aligned. 
Further investigation of this possibility is required, since it could 
prove to be a powerful tool in the determination of the movement history 


of glaciers. 
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SUMMARY OF RECOMMENDATIONS 


1 Fabeiciveriabi lity 

Till fabrics are locally variable, and a spatially limited 
sample (e.g. 2'x2") may be grossly misleading if taken as indicating 
the dominant movement direction of the glacier at that locality. Pebbles 
should be taken from several points spaced as widely as possible across 
the exposure to estimate fabric variability. Combining them into a 


single sample may indicate the dominant movement direction. 


2 Size of pebbles 

Present study indicates larger pebbles have less scatter in 
orientation than smaller pebbles. Hence the largest available stones 
should be used for determining ice movement direction. The sizes of all 


stones used should be recorded. 


3 Shape of pebbles 

Present study imei catee a:b ratio may affect orientation of 
pebble. A pilot study should be carried out in which pebbles are grouped 
according to size and a:b ratio and possibly other axial ratios, so that 
any relationship between shape and orientation can be detected. The 
most suitable combination of axial ratios and size for the purpo - of 
the study can then be selected, and the efficiency of subsequen rk 
greatly increased. The present study indicates that pebbles with a:b 
ratios less than about 1.3 are not good indicators of ice movement direc— 


tion, but this may not be true in other areas. The axial lengths of all 
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pebbles used should be recorded. 


4 Facetting 
Many glacial stones have a "pointed end" and a “blunt end!. 
The direction in which the pointed end of such a stone points should 


be recorded, as this will help to determine the sense of ice movement. 


5 Pebble density 

The density of pebbles in the till may. affect the reliability 
of the pebble orientations as indicators of ice movement direction, and 
an estimate of the pebble density should be recorded at each sample 


site. 


6 Sample size 

Minimum adequate sample size will vary according to fabric 
variability, amount of scatter in pebble orientations, and other factors. 
|t should be determined by a procedure of sequential sampling at each 


sample site, 
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ICE MOVEMENT DIRECTIONS 
AND FABRIC PATTERNS OF GLACIAL SEDIMENTS 


IN THE EDMONTON AREA 


Local Stratigraphy 


Bedrock 


The bedrock of the area consists of poorly consolidated Upper 
Cretaceous bentonitic shales and sandstones with some coal seams and 
bentonite beds. The beds dip at about 20 feet per mile to the south- 
west (Westgate, 1969). Thalwegs of the major valleys in the bedrock 


surface are shown in Figure 17. 


Quaternary Deposits 


The succession of Quaternary deposits in the area is shown in 
Figure 18. These deposits have been described in detail by Westgate 
(1969). The Quaternary landforms and surficial deposits of the district, 
shown in Figure 19, have been described by Bayrock and Hughes (1962). 

The oldest Quaternary deposits are fluvial gravels and sands known 
as the Saskatchewan Gravels. Vertebrate fossils suggest that the young- 
est beds of this formation are of late Pleistocene age (Reimchen, 1968). 
A greyish-brown, dense, Clay-loam till with some inclusions of stratified 
sand overlies the Saskatchewan Gravels or lies directly on bedrock. It 
varies in thickness from zero to more than 20 feet, is highly jointed, 


and possesses some folded joint surfaces. The exact age of this 
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Thalwegs of preglacial valleys in the Edmonton area 


Figure 17 
and locations of measured sections. 
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lower till is not known but it is probably of early Wisconsin age. 
Stratified sediments, 40 feet thick in places, and referred to 
as the Tofield Sand, commonly separate this till from an overlying 
till. This upper till, which is loamy and yellowish brown, has a 
pronounced columnar structure; the greatest observed thickness is 25 
feet. This till is believed to be of late Wisconsin age. 
Up to 50 feet of well-bedded Lake Edmonton sands, silts and 
clays rest on the compact upper till. In some areas, however, deformed 
stratified beds with lenses of till-like material, interpreted as 


ablation deposits, immediately underlie the lacustrine sequence (fig. 18). 
Stratigraphic Setting of Fabric Samples 


Detailed location and stratigraphic information for all sections 
will be found in Appendix B, and the fabric diagrams appear along with 
simplified stratigraphic sections in Figures 6, 7, 9, 10, 12, and 13 
(in pocket). In this section, all other relevant information, such as 
structural observations, will be listed, and consideration given to 
stratigraphic interpretation where necessary. 

Fabric samples were taken mostly by the writer from 41 sites at 
11 locations (A to L except C, fig. 20). These samples (and corresponding 
sites) are numbered 1 to 45 with the numbers 14, 15, 16 and 19 not used. 
12 more samples, including one from an additional location, M (fig. 20), 
were obtained as part of a research project by students in the Geology 
466 class at the University of Alberta in the fall of 1969, and are 
numbered 466-1 to 466-12. 


Each sample will be assessed to establish what reliable inferences 
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Figure 20 Locations of measured Sections. 
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can be made about the distribution of pebble axis orientations at that 
site. The degree of reliability of such inferences will also be con- 
sidered. The techniques used are: (1) Comparison with test samples 
from a uniform distribution and a uniform girdle. (2) If the distri- 
bution is not similar to either of the two models considered in 1, 
comparison with samples 34, 35 and 37, each of 100 observations, or 
with subgroups of these, is found helpful in establishing the form of 


the distribution. 
Location A 


Both upper and lower tills are present at this section, being 
separated by 12 feet of Tofield sand, but only the fabric of the upper 
till was determined (sample 12, fig. 6). If sample 12 is compared with 
the test samples from hypothetical models shown in Figure 5, it is 
readily seen that it differs from samples 505, 506 and 507 from popula- 
tions distributed uniformly over the hemisphere in lacking orientations 
with steep dips. In this respect it resembles samples 508-512 from a 
uniform girdle. However, sample 12 has a density of more than 8e in a 
three per cent area Where e is three per cent, the expected density for 
a uniform distribution), while the highest density in samples 508-512 
(all using a three per cent area) is less than 7e and this appears only 
once. Apart from this high density, sample 12 resembles samples 508-512. 
It is inferred that the a-axes of pebbles at site 12 are concentrated 
into a near-horizontal plane, and in this plane have a preference for 
north-northeasterly and south-southwesterly trends. It is difficult to 


make any inferences about the strength of the preferred orientation. 
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If a density of more than 6e in a three per cent area is to be expected 
in a sample of 50 orientations from a uniform girdle, then probably only 
a slight preference for one orientation in that girdle could result in 
a density of 8e to Ye. In other words, the degree of preference may be 


weak or strong; 50 measurements are insufficient to determine this. 


Location B 


At this location approximately 50 to 55 feet of till is found 
between the underlying Saskatchewan Gravels and the overlying Lake 
Edmonton sediments. Although both upper till and lower till appear to 
be present here, the unconformity is not clearly visible and its place- 
ment at 10 feet below the base of the Lake Edmonton beds is tentative. 
The fabric diagrams and simplified stratigraphic section are shown in 
Figure 10. A total of 11 fabric samples are available from the till; 
their stratigraphic positions are shown in Figure 10. 

Sample 2 shows a preferred orientation indicated by the density 
maximum at 033 24. There is probably a second, weaker preferred orienta- 
tion at 330 18. Sample 33 indicates a preference for trends in the 
northwesterly quadrant, all parts of the quadrant apparently being 
equally preferred. Sample 466-4 shows a strong preference for the 
orientation 036 00. Sample 466-6 has a near-horizontal girdle, but no 
preferred orientation. Sample 466-12 shows a preferred orientation at 
008 14, and a weaker preferred orientation at 276 10. Sample 466-7 shows 
a clear but quite disperse grouping about approximately 068 30. Sampie 
30, when compared with test samples 505-507 (fig. 5), appears to be non- 


random on the hemisphere by virtue of the almost complete absence of axes 
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in a wide band from the northeasterly quadrant to the southwesterly 
quadrant. The centre line of this band would be the trace of a plane 
dipping steeply to the northwest; it is therefore inferred that there is 
a preferred orientation of pebble long-axes dipping slightly to the south- 
east. This inference is supported by the appearance of a density value 
of more than Se at 155 41; a density of this magnitude did not appear on 
any of the three test samples 505-507. Sample 29 has a near-horizontal 
girdle; the double break in this girdle trending northwest-southeast 

and the near break in the northeasterly quadrant suggest a preference 
for north-south and east-west trends, but this is not well indicated. 
Sample 28 shows a preference for the position of approximately 001 04. 
Sample 27 indicates a preferred orientation at 020 00, and sample 466-2 


a relatively strong preference for the position 025 12. 


Location D 


At this location, only the top few feet of the lower till was 
exposed (fig. 6). It was overlain by 10 feet of sand and silt, which was 
in turn overlain by a complex, partially deformed sequence of brown dia- 
mictons interbedded with silts and fine- to medium-grained sands. Sample 
1 was taken from the lowermost diamicton in this sequence. The apparent 
lateral equivalents of these diamictons a few hundred feet east have a 
granular texture with granules up to about one half inch, so that they 
resemble 'mud pebble conglomerates". This fact, and the interbedding of 
the diamictons with silts and fine sands, suggest that the diamictons 
are lacustrine turbidites and not tills. For this reason sample 1 is 


not considered in the discussion of ice-movement directions. 
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The fabric of sample 1 shows a preferred orientation at 351 07 and 
possibly a weaker one at 057 00. Sample 13 is from the top of the lower 
till; it indicates two preferred orientations, the stronger at 134 16, 


the weaker at 021 00. 


Location E 


Both tills are present at this exposure, with no intervening beds 
(fig. 6). A large mass of Cretaceous Edmonton Formation is contained 
within the upper till at one end of the exposure. 

The upper till fabric, sample 3, shows a preference for north- 
easterly trends, possibly for the two positions 008 19 and 066 19. 
Interpretation of sample 4 is uncertain; it may indicate a preferred 
orientation at 090 20, may indicate only a preferred plane dipping 
28 degrees towards 029 degrees, or may be random. The lower till fabric, 


sample 5, shows two preferred orientations at 032 05 and 138 05. 


Location F 


At this locality the upper till is overlain by Lake Edmonton sedi- 
ments and rests directly on lower till. One fabric sample (37) was 
taken from the base of the Lake Edmonton sediments (fig. 6), nine from 
the upper till (fig. 12) and two from the lower till (fig. 6). The 
samples from the upper till were described in connection with fabric 
variability. Sample 37 from the Lake Edmonton deposits shows a pre- 
ference for a near-horizontal plane, and comparison with test samples 


513 and 514 (fig. 5) does not justify inference of any preferred direcs 
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tion in this plane. 

Sample 7 from the lower till shows preference for a near-horizontal 
plane and a well-defined preferred orientation at 044 04. A weaker pre- 
ferred orientation at 308 03 is suggested. Sample 21, also from the 
lower till, indicates only a preference for trends in the northeasterly 
quadrant. 

The base of the lower till at this locality consists of a clayey 
layer that is highly sheared (Plate 2-A, B, C). The top of this highly 
sheared layer is marked by a prominent, flat, shiny shear plane. Slick- 
ensides on this shear plane form two major sets with trends of 033+6 and 
000+2 degrees. A less prominent set of slickensides has a trend of 
about 294 degrees; curved slickensides were also observed. 

The base of the lower till is undulated in places. Two adjacent 
undulations with a wavelength of about two feet had trends of 335 and 
337 degrees. A third undulation about 500 feet from the other two had 


a trend of 310 degrees. 


Location G 


Samples 32, 466-5, and 31 are from the upper till at this locality 
(fig. 7). Up to two feet of lower till separate the upper till and 
Saskatchewan Gravels, but no samples are available from it. Both samples 
32 and 466-5 indicate a preference for southeasterly trends with a high 
degree of horizontal scatter. Sample 31 reveals a preference for low 
plunges, but no preferred trend. 

At one point where the upper till rests directly on sand of the 


Saskatchewan Gravels, its base forms well-defined groove molds that 
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trend from 028 to 033 degrees. 


Location H 


Here, about three feet of upper till rests directly on about 
three feet of lwer till, and is overlain by about 10 feet of Lake 
Edmonton deposits (fig. 7 and Plate 2-D). The exposed face is about 
140 feet long and trends east-west. Fabric samples 10 and 11 are from 
the Lake Edmonton sediments, while samples 8 and 9 are from the lower 
tat. 

Sample 10 shows preferred orientations at 016 06 and 282 05. 
Sample 11 has a near-horizontal girdle, with two gaps trending 
northwest-southeast, a maximum at about 273 05, and a weaker maximum 
at about. 224.05; 

Groove molds on the base of the upper till (Plate 2-E) trend 
035 degrees. Most of the lower till is dark grey, has a fractured, 
blocky structure and contains numerous irregular masses of clean, fine- 
to medium-grained sand (Plate 2-F). Part of the lower till at one end 
of the exposure, however, is brown and has a fine granular texture that 
contrasts with the blockiness of the rest of the till (Plate 2-G). 

Sample 9 (fig. 7) comes from the blocky lower till at the middle 
of the exposure. It shows a near-horizontal girdle with a preferred 
orientation of 036 05. Sample 8 was taken from the granular-textured 
lower till; it also has a girdle with a preferred orientation, but in 
this case the girdle dips 60 degrees toward 021 degrees, and the pre- 
ferred orientation is at 322 48, Figure 21 shows the distribution of 


the poles of 47 joint planes in the blocky lower till close to the site 
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Figure 2]. Poles of joints in fractured lower till at Location H. 
Diagram prepared using a 3% counting circle. Contour interval 
=o = 57, The expected density \for a untform distribution. N= 


number of poles plotted. North at the fop. 
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of sample 9. It has two maxima at 325 23 and 250 12. There is no 
apparent relationship between these maxima and the fabric of the lower 


tall. 


Location I 


At this locality about 11 feet of upper till rests directly ona 
15-foot thick sequence of three dark grey diamicton members separated 
by two sand, silt and clay members (fig. 9). Sample 17 is from the upper 
till, samples 18 and 44 are from the middle dark grey diamicton, and 
sample 26 is from the lower dark grey diamicton. The horizontal distance 
between samples 18 and 44 is 10 feet. 

Sample 17 has a girdle dipping 30 degrees toward 205 degrees. 
There appears to be a preference for west-northwesterly and east-south- 
easterly trends, though not pronounced. Sample 18 has a relatively 
strong preferred orientation at 327 29 and a preferred plane dipping 
31 degees toward 356 degrees. Sample 44 shows a preference for the 
orientation 350 11 and for a plane dipping 12 degrees toward 004 degrees. 
Sample 26 indicates only a preferred orientation of approximately 342 15. 

The interbedding of the dark grey diamictons with sand, silt and 
clay suggests that this part of the sequence is of lacustrine origin. 
The lower diamicton has a granular structure similar to that seen in the 
brown diamictons at Location D; this structure is believed to be indica- 
tive of a subaqueous gravity mass flow deposit (turbidite). The entire 
sequence of diamictons, sands, silts and clays between the Saskatchewan 


Gravels and the upper till is regarded as probably lacustrine in origin. 
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Location J 


Here about two feet of upper till is separated from about 15 feet 
of lower till by one foot of sand. Sample 20 is from the upper till 
and sample 25 from the top of the lower till (fig. 9). 

Sample 20 has a girdle dipping 10 degrees toward 344 degrees. 
Northwesterly and southeasterly trends are preferred, with considerable 
horizontal scatter. Sample 25 has a peak at 346 41, and a smaller peak 
at 206 11. The "best fitting'' plane dips 52 degrees toward 288 degrees. 
The lack of continuity of the density high suggests that there is not 
a uniform girdle, but that the density peaks indicate two preferred 


positions. 


Location K 


A sketch of this exposure showing its structural features appears 
in Figure 22. The bedrock surface slopes quite steeply (10-20 degrees) 
downward to the south. Isolated remnants of massive lower till and of 
rhythmically bedded brown to buff clay rest directly on the bedrock sur= 
face. These are overlain by up to 15 feet of upper till which, in its 
upper part, contains thin sand beds that are folded as indicated in 
Figure 22. It is overlain by clean fine- to medium-grained sand. This 
sand and the folded upper till are truncated by a nearly level uncon- 
Formity on which rests a diamicton about three feet thick that thins 
rapidly southward. 

Fabric samples 22, 23 and 24 were obtained from the points in- 


dicated in Figure 22. Sample 22 (fig. 9) shows a preference for a plane 
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dipping 50 degrees toward 359 degrees and a preferred orientation of 
068 23. Sample 23 is from a few feet above the base of the upper till. 
Comparison with test samples 505-507 (fig. 5) reveals no reason for 
inferring other than a random distribution of the pebble axes. Sample 
24 is from the uppermost diamicton above sample 22. It clearly shows 
a relatively strong preference for northerly to northeasterly trends. 
In view of the high concentrations of observations in this region of 
the projection, it is thought probable that the bimodal nature of the 
sample reflects a bimodality in the population. 

The near parallelism of the fabric maximum from the upper part 
of the upper till to the lineations on the bedding surfaces suggests 
that both the fabric and the lineations were produced by the same 
movement. This movement could have been either a series of subaqueous 
gravity mass flows (density currents) in a lacustrine environment, or 
shearing due to ice movement in a sub-glacial environment. In the latter 
case, the sand beds could perhaps be explained by the "'undermelt theory'' 
of Carruthers (1953). 

The origin of the uppermost diamicton is not clear. Its strati- 
fication, granular structure and rapid thinning southward suggest a 
local gravity flow. Such a flow could have produced the grooves on 
the base of this unit. However, the truncation of the folded upper till 
with its overlying sand (although possibly due to erosion by density 
currents) suggests overriding by a late glacial readvance. Evidence of 
such a late readvance is seen in the form of moraines in the Redwater 
area northeast of Edmonton. It is not known at this time whether this 


readvance extended as far as Location K. 
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Location L is a large gravel pit where numerous exposures of both 
tills have been made. No sediments intervene between the tills. Four 
fabric samples are available from the upper till and three from the 
lower till as indicated roughly by the sketch in Figure 13. 

The lower till contains three structurally distinct vertical zones. 
At its base is an intensely sheared, continuous, black clayey till layer 
one to two inches thick. Plate 1-B shows the basal portion of a till 
block collected from the base of the lower till at site 45. The contact 
with the underlying Saskatchewan Gravels can be seen, and the basal 
sheared zone is clearly discernible. At one point within this layer, 
near-horizontal slickensides were observed, having a trend of 015 de- 
grees. 

The middle zone is dense and massive, and has a columnar structure 
(Plate 1-D, E). Its thickness varies from zero to about four feet. 
Samples 45 and 466-10 (fig. 13) are from this zone. Both show a near- 
horizontal girdle and a preference for northwest trends with fairly 
high horizontal scatter. 

The upper zone of the lower till is highly fractured and contains 
many shear surfaces (Plate 1-D, E). Lineations (Plate 1-C) and grooves 
(Plate 1-D) on near-horizontal shear surfaces have trends between 020 
and 025 degrees. Sample 43 is from this zone. Preference for north- 
northwesterly and northeasterly trends is apparent. The two groupings 
in the sample probably represent two preferred orientations in the popu- 
lation, since the high concentration of observations in this region 


affords a better than usual degree of reliability. The two groupings 
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are centred at approximately 346 35 and 016 10. 

Sample 466-1 from the upper till shows two well-defined preferred 
orientations: one at 340 16 and one at 031 10. Sample 466-3 shows a 
preference for trends in the northerly and southerly quadrants between 
northwest-southeast and northeast-southwest. The maximum at 213 06 may 
represent a preferred orientation. Sample 466-8 indicates a preferred 
orientation of 209 11 and a preference for a near-horizontal plane. 
Sample 466-11 apparently indicates two preferred orientations at 010 06 


and 075 38. 


Location M 


At this gravel pit only one till, the upper till, was present in 
the exposure mapped. Sample 466-9 (fig. 7) shows a preference for a plane 
dipping 06 degrees toward 310 degrees, but no reliable indication of a 


preferred trend. 


Upper Till Fabrics and Glacier Movement 


Well developed groove molds on the base of the upper till in the 
Clover Bar area just east of Edmonton (Plates 1-A, 2-E) trend from 028 
to 035 degrees and indicate that the sense of movement was from northeast 
to southwest (Westgate, 1968). Such a line of movement for the glacier 
that deposited the upper till is also indicated by flutings (Bayrock and 
Hughes, 1962). 

The preferred trends of elongate pebbles in the upper till, in- 


ferred from the present study as discussed in a previous section, are 
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indicated in Figure 23. To these have been added sample AFL (Lorberg, 
1967) and sample RBR-1b (Rains, 1969) measured by other investigators. 

(Rains (1969, p.18) tested the aximuth distribution of sample 
RBR-1b, using 18 ten-degree classes, against uniformity using chi-square 
and found no significant deviation. However, while a high value of 
chi-Square can justify rejection of the hypothesis model, a low value 
of chi-square does not necessarily justify acceptance of the model. 

The low chi-square value calculated by Rains does not justify acceptance 
of the model of uniformity. The writer tested the same sample using a 
different set of classes and obtained a chi-square value that did 
justify rejection of the model of uniformity.) 

Of the seven samples from the upper till at Location F, five show 
well-defined preferred trends. These five preferred trends range from 
000 degrees to 066 degrees (fig. 12). The four samples from the upper 
till at Location L show six preferred trends varying from 340 degrees 
to 075 degrees (fig. 13). One sample is available from the upper till 
at each of the other locations shown in Figure 23; except for those at 
Locations G, I, and J, these show preferred trends ranging from 008 
degrees to 066 degrees. It cannot be assumed that the degree of fabric 
variability seen at Locations F and L does not exist at these other 
locations. In the absence of information about the variability of the 
fabric at these locations, the preferred trends observed in single 
samples cannot be regarded as accurate indications of ice-movement 
directions. 

All the preferred trends shown in Figure 23 lie in the range 340 
degrees to 075 degrees, except for samples 33 (Location bs tide 10), 32 


and 466-5 (Location G, fig. 7), 17 (Location I, fig. 9) and 20 (Loca- 
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Figure 23 Map showing preferred trends of elongate pebbles in the 
upper till. 
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tion J, fig. 9). The preferred trends shown by these five samples 
are all west-northwest-east-southeast to northwest-southeast, that is, 
approximately normal to the direction of movement of the glacier as 
indicated by sole markings and surface features. Thus they may represent 
pebbles oriented transverse to the movement direction. All five 
samples show a high degree of horizontal scatter; it is possible that 
this is characteristic of transverse fabrics. 

Striations on a shear surface within the clayey layer at the 

base of the lower till at Location F indicate two principal lines of 
movement: 033+6 degrees and 000+2 degrees. It will be shown in the 
later discussion of the lower till that these striations are probably 
related to movement of the glacier that deposited the upper till. The 
inference of two principal directions of movement for this glacier is 
supported by the orientations of the more elongate pebbles from the 
upper till at Location F (fig. 14 and 15), which show preference for’ 
the positions 022 07 and 349 22. Moreover, as pointed out previously, 
as the elongation of the pebbles increases (fig. 14, top row) the group- 
ing at 349 19 increases in strength while that at 029 15 becomes weaker 
and finally disappears. If it is supposed that the more elongate pebbles 
will be the first to become re-aligned to a new movement direction, then 
the above observation indicates that movement from about 349 degrees 
occurred later than movement from about 022 degrees to 029 degrees. This 
1S consistent with the striations described above. The weakening of the 
grouping at about 055 27 with increase in elongation of the pebbles may 
be regarded as suggesting that the movement from 022 degrees to 029 
degrees was preceded by still earlier movement from about 055 degrees. 


Figure 15B shows the distribution of c-axes of flat pebbles, again 
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from the upper till at Location F. The majority of these c-axes are 
nearly vertical, confirming the observations of Harrison (1957a) and 
Johansson (1968). A second group of c-axes is seen to trend south- 
easterly to easterly with plunges from 34 degrees to 18 degrees, 
indicating a group of flat pebbles with their a:b planes striking 

From northeast to north and dipping steeply to the northwest or west. 
If it is supposed that the steeply dipping azb planes of such flat 
pebbles will tend to assume an orientation such that they strike nearly 
parallel to the direction of ice-movement, then the observed north- 
easterly to northerly strikes are consistent with ice movement ranging 
From northeast-southwest to north-south. Presumably the azb planes 
initially assumed a near-vertical position with their strikes parallel 
to the early movement direction from northeast to southwest. It is 
suggested that as the movement direction shifted toward a north-south 
line, some of these a:b planes became re-aligned to strike parallel to 
the new movement direction; those that did not become re-aligned were 
tilted away from their vertical attitude as a result of shearing or 
Flow of the till. The longer the period of movement oblique or trans- 
verse to the original strike of an azb plane, the more it would tend to 
be tilted if it did not become re-aligned to strike parallel to the new 
movement direction. Thus, asb planes with the "oldest" strikes would 
be expected to have the greatest deviation from the vertical, or the 
smallest dips. Such a relationship is seen in Figure 15B if the suggest- 
ed history of ice movement is correct: the a:b planes stiking north- 
easterly (parallel to earlier movement) dip less steeply than the asb 
planes striking northerly (parallel to later movement ). 


One possible explanation of the apparent change in ice movement 
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direction from northeast-southwest to approximately north-south at 
Location F is the increasing influence of local topography on the 
glacier as it thinned. The bedrock valley trending north-south about 
one mile west of Location F (fig. 17) may have been responsible for 
channeling the thin waning glacier into a north-south flow pattern. 
Such late-stage changes in flow direction due to increasing influence 
of local topography would clearly be different in different parts of 
the study area. Johansson (1968, p.210) believed that control of ice 
movement by local topography was related to the existence of a relative- 
ly thin ice cover. 

The wide spread of preferred trends seen in the upper till in the 
Clover Bar area may thus be due to the wide range of movement directions 
to which this area may have been subjected as a result of the local 


topography, and may not be found at other localities. 


Lower T1711 Fabrics and Glacier Movement 


The preferred trends shown by samples believed to be from basal 
lower till are shown in Figure 24. Fabrics of units thought to be 
turbidites were excluded. Sample 8 was also excluded; its steeply 
dipping girdle, being uncharacteristic of till fabrics, suggests that 
the granular-structured ''phase'' of the lower till in which it was 
measured is an inclusion of previously deposited sediment retaining 
its original fabric but having been rotated. Before attempting to 
infer any ice movement directions, the first fact that must be con- 
sidered is that the pebble fabric of the lower till is thought to 


have been altered by the overriding later glacier at several of the 
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Figure 24 Map showing preferred trends of elongate pebbles in the 
lower till. 
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locations studied. The evidence for such alteration at Location L 

has been described in detail elsewhere (Ramsden and Westgate, in press) 
and is briefly as follows. The fabric (fig. 13, samples 45, 466-10) of 
the middle zone of the lower till (Plate 1-B,D,E) is thought to be the 
original depositional fabric of the lower till for two reasons: Cp)isite 
consists of a near-horizontal girdie with a preferred direction in this 
girdle, a pattern described by several workers from till in areas of 
ground moraine (Holmes, 1941, p.1313; Harrison, 1957a, p.283; Kauranne, 
1960, p.87), and (2) the middle zone is much less sheared and fractured 
than the upper zone, it is comparatively massive and dense, and its 
columnar structure resembles that of the undisturbed upper till. 

Both structural and fabric observations are important in inter- 
preting the fabric of the upper zone of the lower till. The near- 
parallelism of the slickensides in the basal sheared zone of the lower 
till (with trends of between 015 and 035 degrees) to the grooves and 
lineations in the upper zone (Plate 1-C,D) (with trends of 020 to 025 
degrees) suggests that all these similarly oriented features were 
produced at about the same time by the same agent of deformation. 
Furthermore, the close parallelism of the trends of these features 
to the known direction of movement of the later glacier (fig. 13, 
top row; Plate 1-A; Plate 2-E) strongly suggests that the movement of 
this glacier was the cause of the deformation in both zones. 

The grooves and lineations in the upper zone of the lower till 
are also nearly parallel to the trend of about 016 degrees of the fabric 
maximum of this zone (fig. 13, sample 43). This suggests that shearing 
altered the original fabric of the till, producing the present con- 


figuration. The other maximum in the fabric of this upper zone may 
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represent a remnant of the northwesterly trending maximum of the original 
fabric. 

Other fabric samples from the lower till are thought to represent 
fabrics partially altered by the overriding later glacier. Sample 13 
(fig. 6), from the top of the lower till at Location D, has a secondary 
maximum at 021 00 that may be due to overriding by the later glacier. 
Since the diamictons that overlie the lower till at this location are 
thought to be lacustrine turbidites, it cannot be verified that the later 
glacier actually moved from 021 degrees at this place; however, it is 
not unlikely in view of the upper till fabrics in the area (fig. 23). 
Sample 25 (fig. 9) from the top of the lower till at Location J has a 
maximum at 346 41 and a secondary maximum at 206 11, and may represent 
another partially altered fabric. If the preferred trend of about 306 
degrees in sample 20 from the upper till at this location represents 
pebbles transverse to ice movement as it is thought, then it is consis- 
tent with the postulated partial alteration of the lower till fabric. 

In addition to the above-mentioned fabrics that are thought to be 
partially altered, two of the samples from the lower till are thought to 
represent fabrics completely altered by the overriding later glacier. 
This inference is made on the basis of parallelism or near-parallelism 
of the fabric maximum to the known direction of advance of the later 
glacier, combined with intense fracturing and shearing of the lower 
till parallel to this direction of advance. Such fabrics are samples 
7 and 21 from the lower till at Location F (fig. 6; Plate 2-C). Both 
of these fabrics have a near-horizontal girdle and a preference for 
trends of about 045 degrees. Sample 7 shows a transverse peak. Thus 


these altered fabrics show the same characteristics as primary till 
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fabrics.e Other samples taken from highly fractured lower till show 
maxima with trends within the known range of movement directions of the 
later glacier, but where shearing of the lower till parallel to this 
direction could not be demonstrated beyond doubt the origin of the lower 
till fabric cannot be determined. Indeed, perhaps the biggest single 
problem in Ae roetina lower till fabrics in the Edmonton area is that 
of differentiating between primary fabrics and those produced by the 
overriding later glacier. 

The preferred trends of lower till fabrics were re-plotted in 
Figure 25, omitting those described above as wholly or partially alter- 
ed by the later glacier. The remaining preferred trends at Locations 
H, L, E, and D, those on samples RBR-3, RBR-4, and RBR-5, and the more 
easterly maxima on samples RBR-1la and RBR-2a all lie in the two groups 
118 degrees to 138 degrees and 032 degrees to 050 degrees. These two 
intervals are centred at approximately 130 degrees and 040 degrees, 
that is, are close to 90 degrees apart. It is possible that these two 
directions are parallel and transverse to the dominant movement direction 
of the glacier that deposited the lower till. The maxima in samples RBR- 
la and RBR-2a that do not fit the northeasterly trending group could 
have been produced by the later glacier, although no structural informa- 
tion is available from these sites at this time. 

The wide range of preferred trends at Location B may be due to 
changes in movement direction resulting from varying degrees of influence 
of the local topography as the thickness of the ice-sheet changed. Sit- 
uated on a major bedrock valley close to its confluence with other major 
valleys, the location would be susceptible to such changes. 


Only two of the samples whose preferred trends lie in the two 
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Figure 25 Map showing preferred trends of elongate pebbles in the 
lower till, with wholly or partially altered fabrics omitted. 
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groups described above are accompanied by any positive evidence that 
they represent primary lower till fabrics. These are samples 45 and 
466-10, both from Location L (fig. 13), and the evidence is mainly the 
contrast between the structural characteristics of the middle zone of 
the lower till from which these samples where taken and the sheared 
upper zone. It may be worth noting that both these samples show a well- 
defined near-horizontal girdle but a high degree of horizontal scatter 
about the preferred trend. As mentioned previously, the five upper till 
fabrics that have a similar pattern are thought to be transverse fabrics, 
suggesting that the glacier that produced fabrics 45 and 466-10 in the 


lower till was moving from northeast to southwest. 


Value of Fabrics as a 


Stratigraphic Tool 


It can be seen from Figures 23 and 24 that no preferred trend is 
found exclusively in one till unit. Since both tills show a very wide 
range of preferred trends, determination of the preferred trend of the 
pebbles in an unidentified till cannot serve to identify the till. 
Three samples were taken from the Lake Edmonton diamictons -- two at 
Location H, several feet above the lower limit of observable stratifi- 
cation (fig. 7, samples 10, 11), and one at Location F about one foot 
above the lowermost visible stratification (fig. 6, sample 37). These 
three fabric patterns show no apparent characteristics that would allow 
them to be singled out as being different from the rest of the samples. 
However, groups 37-1 and 37-2 of sample 37 (fig. 11) indicate a sub- 
stantial change in the fabric over a distance of less than a foot; the 


sequential subgroups of samples 34 and 35 from the upper till at the 
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Same exposure show no such change. Thus this extreme variability may 
be characteristic of the lake sediments and perhaps of lacustrine dia- 
mictons generally, but no other data are available to illuminate this 
problem further. 

Two samples were taken from units regarded as possible turbidites 
on the basis of structural and (or) stratigraphic evidence. These samples 
are 1 (Loc. D, fig. 6) and 8 (Loc. H, fig. 7). Five more samples from 
two more locations are from units regarded as probable turbidites; these 
are 18, 26 and 44 from Location I and 22 and 24 from Location K (fig. 9). 
These samples all show a relatively strong preferred orientation, lack a 
continuous girdle (except sample 8), and lack transverse maxima (except 
sample 1). To test the possibility that these features are characterist- 
ic of turbidite fabrics, the same features were looked for on other 
samples believed to be from till. Other samples with these features are 
466-4 and) 46622) (Loc. B, fig. 10), 6 (Loc. Fy fig. 12), and 466-1" and 
L66=8r (loct.s Lan figied3)t 

Sample 466-4 is from the top of the upper till immediately below 
Lake Edmonton sediments. It could conceivably represent a lacustrine 
turbidite. Sample 466-2 is from the basal part of the lower till. This 
basal portion of the lower till is partially stratified, contains thin 
sand beds and is quite possibly lacustrine in origin. Sample 6 is from 
the middle of the upper till which is not likely to be a turbidite. 
Sample 466-8 was taken from well below the top of the upper till and is 
also almost certainly a till fabric. The stratigraphic position within 
the upper till of sample 466-1 is uncertain, and it is possible that it 
came from basal Lake Edmonton sediments. Clearly the available informa- 


tion allows no conclusions to be drawn about differentiation of tills and 
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turbidites by means of their pebble fabric. 

Lindsay (1968) simulated 'mudflow'' fabrics using a computer as 
well as studying real ''mudflow'' and tillite fabrics. He concluded that 
whereas the tillite fabrics had only a single mode, mudflow fabrics were 
characterized by a dipping girdle, sometimes with a maximum onthis girdle. 
This contrasts with the results of other studies, including the present 
one, that indicate a near-horizontal or moderately dipping girdle, fre- 
quently with a maximum on it, is the characteristic pattern of till fab- 
rics in areas of ground moraine. Lindsay, et al (1970) noted that all 
"mudflow'' fabrics examined by them had single modes paralleling the flow 
direction, some with and some without a girdle. They thus lacked trans- 
verse modes. All the fabrics of probable turbidites listed above lack 
feabar suse modes. Thus the transverse mode may be peculiar to till 
fabrics and normally lacking in mudflow or turbidite fabrics. 

Glen, Donner and West (1957, p.194) pointed out that ''Laboratory 
experiments and theoretical studies by various authors have shown that 
an initially random collection of elongated objects immersed in a flow- 
ing liquid will develop a long axis distribution with a peak parallel to 
the direction of flow in a very short time, but that if flow is continued 
for a long time a transverse peak will also develop.'' If this is true, 
it would not be surprising to find that turbidite and mudflow fabrics 
lacked the transverse mode, since their period of formation is very much 
shorter than for till fabrics, and the distance of transportation prob- 
ably much smaller, so that they would be much less well developed than 
cid fabrics. 

The fabrics of stratified Lake Edmonton diamictons (samples 10 
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the fabrics of the probable turbidites. All three have a near-horizontal 
girdle; sample 10 is very similar to sample 5 (Loc. E, fig. 6), a prob- 
able till fabric. Sample 11 shows a preference for northeast-southwest 
and east-west trends. Sample 37 probably represents a uniform girdle. 
Microfabric studies might be able to yield far more information 
than pebble fabric studies for the same investment of time, and may 
prove to be a more powerful tool than pebble fabrics in the determina- 


tion of sediment origin. 


Work Needed 


More fabric and other stratigraphic data are needed before the 
direction of movement of the earlier glacier can be determined. The 
effect of local topography on ice-movement directions in the Edmonton 
area could be better assessed if more fabric data were available at 
locations away from the main bedrock valleys. Differentiation of 
primary fabrics and fabrics altered by shearing resulting from over- 
riding by a later glacier could perhaps be achieved through study of 
the micro-structural features of the sediment. Further work is required 
on the fabrics of diamictons of different origins, for example, tills, 


turbidites, and ice-rafted deposits. 
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CONCLUSIONS 


A method of statistically treating till pebble orientations based 
on an assumed probability distribution has been used by several workers 
in recent years. The method was proposed by Andrews and Shimizu (1966) 
and based on the assumption that the distribution of till pebble axes 
may be described by a probability density function due to Fisher (1953). 
The method was examined and found to be unsatisfactory for the following 
reasons: (1) The condition k23, where k is Fisher's precision parameter, 
is used as the sole criterion of applicability of the method, that is, as 
indicating that the distribution is essentially unimodal. This is not 
true. (2) The recommended procedure for deleting transverse pebbles 
fails. (3) Watson's (1956b) test for randomness of directions, based on 
the length of the resultant vector, is not valid for hemispherical dis- 
tributions. 

The failure of the method, adopted from a procedure designed to 
treat paleomagnetic vectors, is attributed largely to the fundamental 
differences between magnetic vectors and pebble axes. Thus a thorough 
investigation of the differences between different kinds of orientation 
data and of the distinct functions of reference planes used in their 
treatment must precede the development of suitable methods of treatment. 

An alternative approach, assuming Fisher's probability distribu- 
tion for each mode of a multi-modal pebble axis distribution, was also 
examined. Testing of till fabric data, however, showed that in the 
majority of cases even the individual modes did not conform to Fisher's 
density function. 


The calculation of descriptive statistics for till fabric samples 
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based on valid three-dimensional models must await further work on either 
(1) models for elongate groupings, or (2) procedures permitting a group- 
ing to be treated separately from a girdle occurring in the same sample. 
Until adequate techniques for the numerical representation of till fabric 
data are devised and widely accepted, till fabrics should be represen- 
ted graphically to avoid possible loss of information. 

Knowledge of the distribution of R, the length of the resultant 
of a set of unit vectors, for a hemispherical distribution would be a 


valuable aid in the interpretation of till fabric samples. 


Fifty pebbies with a-axes of at least one inch and a:b ratios of 
at least 1.2 are barely adequate to define the fabric of a till at a 
single point. Preferred orientations are indicated only approximately. 

Although different samples from a till unit at a given exposure 
are in some cases very similar, this can by no means be counted upon, 
since frequently great variation is evident. It must be concluded that 
a single one-point sample should never be depended upon to indicate the 
preferred trend of elongate pebbles in a till unit at an exposure. 
Indeed, the assumption that the pebbles in a till at an exposure have a 
single preferred trend must be questioned. They may have many preferred 
trends, or different preferred trends in different parts of the exposure, 
both laterally and vertically. 

A sample of less than about 200 long axes from the nine sample 
sites in the upper till at Location F could not be relied upon to indi- 
cate the northeasterly trending preferred orientation. This is attribu- 


ted to the degree of fabric variability of this unit at this location. 
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More elongate pebbles, that is, pebbles with higher a:b ratios, 
appear to be more sensitive indicators of ice-movement direction. They 
are apparently the first to become re-aligned with a new direction of 
ice movement, so that they may be found to prefer a later direction of 
movement while less elongate pebbles still prefer an earlier direction. 

No relationship was found between the size of a pebble and its 
preferred orientation. The size was, however, found to have considerable 
influence on the tightness of grouping of pebbles about their preferred 
orientation. Harrison's (1957a) conclusion that the majority of flat 
pebbles have their c-axes nearly vertical was confirmed. The imbrication 


of these pebbles, if any, could not be determined. 


Five fabric samples from the upper till are thought to represent 
transverse fabrics. They all show little vertical scatter but much hori- 
zontal scatter. This may be characteristic of transverse fabrics. The 
direction of movement of the glacier that deposited the upper till at 
Location F appears to have changed from northeast-southwest to north- 
south during the later part of the glaciation. This change may be due to 
the increasing influence as the glacier thinned of the north-south-trend- 
ing bedrock valley about one mile west of this location (fig. 17). The 
wide spread of preferred trends seen in the upper till in the Clover Bar 
area may thus be due to the wide range of movement directions to which 
this area may have been subjected as a result of the local topography, 


and may not be found at other localities. 


The fabric of the lower till has been wholly or partially altered 


by the overriding later glacier at several of the locations studied. One 
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of the major problems in interpreting lower till fabrics in the Edmonton 
area is that of differentiating between primary fabrics and those pro- 
duced by the later overriding glacier. Insufficient data are available 
to determine the direction of movement of the glacier that deposited the 


lower till. 


Fabric samples from stratified Lake Edmonton diamictons showed no 
features that could serve to distinguish them from till fabrics. Several 
samples from probable lacustrine turbidites, however, suggest that fabric 
data might help to distinguish such deposits from tills. 

More information is required concerning fabric variability, the 
effect of pebble shape on pebble orientation, and the effect of local 
topography on ice movement. 

The potential value of studies of micro-fabrics and micro-struc- 


tural features of tills and other diamictons should be investigated. 
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APPENDIX A 


PLATES 


PLATE 


| 


Groove molds trending-028 to 033 degrees at the base of the upper 


till at Location M. The undercut is about two and a half feet 
deep. 
Ti-lhbleck Collected from thewbase of ihe, lower fill et Location &: 


The basal sheared zone is clearly visible. 
Lineations trending 020 to 025 degrees on a near—horizontal shear 
surface in the upper zone of the lower till (see text). The coin 


iS a dime. 


Undulating shear surface separating the upper zone (b2) from the 


middle zone (bl) of the lower till at Location L. The undulations 
are approximately parallel to the movement direction of the glacier 
that deposited .the upper til) Capproximately 025° fo°055 degrees). 


Photo covers about eight feet vertically at centre. Letters and 
numbers have same meaning as in E. 


The section exposed. af Location L: a. Saskatchewan Gravels. bl. 


Middle zone of lower till. b2. Upper zone of lower till. c. Upper 
till. d. Lake Edmonton sediments. The prominent shear surface 
separating the middle and upper zones of the lower till is indicated 


by dashed Iines. The pick is 17 inches long. 
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A, B 


Stereopair showing the basal. sheared zone of the lower till at 
Location F. The prominent flat shear surface on which striations 
were measured can be seen about one inch above the top of the 
underlying sand of the Saskatchewan Gravels. Vertical separation 
of the two pebbles is about three inches. 


Part of the lower till at Location F showing, its highly fractured 
character. The two pebbles that appear in A and B can be seen to 
the left of the mattock head. The mattock is I7 inches long. 


The section exposed at Location H: a. Saskatchewan Gravels. 
b. Lower (till. ¢. Upper Till. ds Lake Edmonton sediments. The 
rod is four feet long. 


Groove molds frending 035 degrees at the base of the upper fil]. at 
Location: H. . The brush .is about, eight. inches ong: 


The tower till. and:base of the upper 4aail at, Location H,. showing 
inclusions of clean, fine= to medium—grained sand in the lower 
till. The fractured. nature of. the-fower fill (b) contrasts wi th 
the columnar jointing of the upper till (c). The rod is four feet 
long. 


The lower till at Location H, showing normal dark gray, blocky 
lower till below and the brown, granular-—textured "phase", believed 


_ to be an inclusion of previously deposited sediment, above. Diameter 


of large stone is*about one foot. 
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APPENDIX B 


MEASURED SECTIONS 


Location A 
East bank of North Saskatchewan River, SW%, sec.6, Tp.51, R.25, W4. 
Stratigraphy from unpublished report by G. Gabert, Geology Department, 


University of Alberta, 1967. 


Average 
Description Thickness (ft) 
Pleistocene 
Silts and clays, bedded 34 
Till, brown, sandy 2 
Silt and sand, bedded t2 
Till, clayey, dark grey 3 
Saskatchewan Gravels: sand hg 


Cretaceous 


Edmonton Formation + 
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Location B 
North bank of North Saskatchewan River; middle of west side of 


sec. 15, Tp.52, R25, W4, 


Description Average 
Thickness (ft) 


Pleistocene 


Clay and silt, bedded 8 

Till, brown vertical joints 10 

Till, grey-brown to dark grey ke 

Fine sand and silt, current-bedded 2 

Sand, medium-grained, buff, stratified 34 

Saskatchewan Gravels: sand and gravel 11 
Cretaceous 


Edmonton Formation 2+ 
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Location C 


East bank of North Saskatchewan River, middle of west side of 
sec.14, Tp.52, R.25, W4. Stratigraphy from unpublished report by E F 


Lorberg, Geology Department, University of Alberta, 1967. 


Description Average 
Thickness (ft) 


Pleistocene 
Clays and silts, bedded ay | 


Till, dark grey, lighter colours towards 18 
top, columnar jointing 


Saskatchewan Gravels: sand and gravel 66 
Cretaceous 


Edmonton Formation 2 
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Location D 
Excavation for new Biological Sciences building, north end of 
University of Alberta campus; middle of south side of sec.31, Tp.52, 


R.24, Wh, 


Description Average 
Thickness (ft) 


Pleistocene 
Clay, massive, brown, scattered pebbles 10.0 


Clays, silts and fine sands, brown and grey, Bl 
highly contorted 


Silt and fine sand, buff to brown 1.6 
Diamicton, brown, many vertical joints 6.0 
Silt, fine sand, brown to buff, lenses of brown, 4.0 


sandy diamicton 


Diamicton, brown 1.0 
Silt and sand 15 
Diamicton, brown ZO 
Sand, medium-grained, clean, buff 2.0 
Silt and fine sand, brown 4.0 
Sand, fine- to medium-grained 4.0 


Till, clayey, dark grey 2.0+ 
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Location E 


Abandoned gravel pit just east of creek, sec.1, Tp.53, R24, W4, 


Average 


Description 
Thickness (ft) 


Pleistocene 


Silt, clay, fine sand, buff to brown, 8 
scattered pebbles, highly contorted 


Till, sandy, brown, medium blocky structure 

to massive. Inclusions of buff to brown silt 
and sand up to at least 4 ft x 20 ft. 

Includes large mass (larger than 5 ft x 20 ft) 


of Cretaceous bedrock. 26 
Till, clayey, dark grey, fine blocky structure 12 approx. 
6+ 


Saskatchewan Gravels: sand 
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Location F 


Twin Bridges Sand and Gravel Co. Pit No. 2. NE%, sec.6, Tp.53, 


R.23, W4, 


Description Average 
Thickness (ft) 


Pleistocene 
Silt and clay, stratified, brown to buff. 
Scattered pebbles and less stratification 
in lower part 10 


Till, brown, columnar jointing 6 


Till, dark grey, fine to medium blocky 
structure 10 


Saskatchewan Gravels: sand and gravel 30 approx. 


Location G 
Gravel pit owned and operated by G. Kropp in north central part of 


sec.8, Tp,53, R.23, Wh. 


Description Average 
Thickness (ft) 


Pleistocene 


Sand and clay, stratified 7 
Till, sandy, buff to brown. Columnar jointing 14 
Till, clayey, grey. Blocky structure. 1 


Saskatchewan Gravels: sand and gravel 20+ 
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Location H 


Twin Bridges Sand and Gravel Co. Pit No. 1. SW, sec.16, Tp.53, 


R.23, W4. 


Description Average 
Thickness (ft) 


Pleistocene 


Clay and silt, stratified, scattered pebbles. 
Less well stratified and more pebbly in lower 
part. 10 


Till, brown, columnar jointing 3 


Till, dark grey, fine to medium blocky 

structure. Inclusions of clean, medium- 

to coarse-grained sand in upper part. 

Also inclusion of brown diamicton with 

granular structure. 3 


Saskatchewan Gravels: sand 6+ 
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Location I 


| 69 


East bank of North Saskatchewan River, middle of east side of 


seq.34, Tp.53, /R. 


Description 


Pleistocene 


23, W4. 


Silt and clay, stratified 


Till, light brown, columnar jointing 


Diamicton, 


Sand, silt 
4einch bed 


Diamicton, 
Sand, silt 
Diamicton, 


Clay, dark 


dark grey, columnar jointing 


and clay, stratified, includes 
of brown diamicton 


dark grey, columnar jointing 
and clay, stratified 
dark grey, granular structure 


grey, granular structure 


Saskatchewan Gravels: sand and gravel 


Cretaceous 


Edmonton Formation 


Average 
Thickness (ft) 


322 
3.4 


33 approx. 


25+ 
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Location J 


North bank of North Saskatchewan River, middle of north side of 
sec.23, Tp.54, R.23, W4. 


A St ng Sse enn se 


Description Average 
Thickness (ft) 


Pleistocene 


Diamicton, stratified 8.7 

Silt 2.0 

Till, brown, columnar jointing 4.0 

Sand M49 

Till, grey, blocky structure, sand 

and gravel inclusions 15 approx. 
Cretaceous 

Edmonton Formation ot 


Location K 
Gulley formed by surface runoff, on Pointe-aux-Pins Creek in 


extreme SW corner of sec.32, Tp.53, R.22, W4. 


Description Average 
Thickness (ft) 


Pleistocene 


Diamicton, brown, stratified, granular structure 2 

Till, brown, beds of sand in upper part 10 approx. 

Till, dark grey, massive 1 
Cretaceous 


Bearspaw Formation 2+ 
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Location L 


Gravel pit operated by Alberta Concrete Products Limited. NW%, 


sec.8, Tp.53, R.23, W4. 


Average 


Description 
Thickness (ft) 


Pleistocene 


Silt and clay, stratified 9 
Till, brown 10 
Till, dark grey 7 
Saskatchewan Gravels: sand and gravel 10+ 


Location M 


Gravel pit operated by J. Pawluk and Son. NW%, sec.8, Tp.53, 


R.23, W4. 


Description Average 
Thickness (ft) 


Pleistocene 
Silt and clay, stratified Zen 
Till, brown, columnar jointing 12.5 


Saskatchewan Gravels: sand 5+ 
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FABRIC 34 (Dimensions in millimeters) 
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FABRIC 34 (continued) 
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FABRIC 34 (continued) 
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FABRIC 34 (continued) 
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FABRIC 35 (continued) 


Pebble Axial lengths Axial ratios a_axis ab plane 
number 2 AR |e a:b bic trend pit . ste. dip 
O3 | 050". O22? Ol 1.4? > 20 [ihe5* > RF Oye a14e Se 
O32 024: OTF . OOS fae IG O07 O09 O08" +52 SE 
033 O40’ 040° 025 id - le 

034 Ore O15", O09 uae ib OP, Is 

035 OFF -OrK2’ O06 [eee [RS OMe « 26 

036 030’ O24 O14 [Reoe ihe? 252 49 Ive" 62° SW 
O37 ' O45? OSPF - O50 Pee lee O27) : 54 

038 O50. O22 - OF [aes [PG OOS? . 2 213 85 NW 
O39 OZ> Ore’ . O48 fae le Q29* . 25 259 34 NW 
040 055’ . 050° O24 Po? Wh 

O4 | O50" O25’ . O21 2s [RZ eo. 

042 Ory OS’ OOF io ee eS L547. (© 092 |5 SW 
043 O51" O20 O25 ll [PZ 

044 O25’ OFF O44 ce O46 353i O25. Secse 
045 O19. O12 O09 or 1S 220' - IM 

046 O29? O19" O1F fier ee KOZ . 24 

O47 OFT: O14 .OO9 oet.’ RO (S57 40 

048 039. 025 «O17 ee > TS O48 . I2 

O49 O29): 025 O20 Meo. ~-lieae | 

050 Or Olly “O08 ean. Alea So" AZ 

O5\| OSG 032 O21 Woe txt eo OO|l 29 

052 020° Of7 009 2 eo O04] 48 O24* (82-SE 
053 O25 O18 .OT5S lor lee arr. - Re 

O54 O25 Oley.. O14 Hee Aas Ale 269 64 

G).5)5) O34) 032 O24 al iD O52. (OVP 
056 O32. 024) 02 [oy leat 342 06 

a)ey O69 -0356-. 05" io he Og = 35 

058 069 043 O04! Poor. 1 199 30 

059 Opa". O10, O05 Se - 2.0 552° 

O60 025°, O13. 008 oA we Soe) - TO 

06] OFS. OS’. OO 4s TG ee 

062 Ol, O10: OO8 ee mM i, C7 26s 

063 Oise Ol: 005 Pao 28) 60S .42 

064 OTE", O06" 005 [ash LEG Zee : 05 

065 OS) O25 O16 oh? ee 254 57 

O66 O20 OTS - 008 re. IEG ‘> I 

O67 020 018 . 005 hal D6 Ze.) So 

068 OS OR” OTS faite [ed S56) a / 

O69 059 046 .058 Se leer 25575 29 

O70 O25". O20. 7,070 Roo - 2s0 OGG? , 17 


ny yl f 
ie 
fi 7 4 
0 : Pe 


ta) i 


THO a 


tee yeaa 100 


32 $8 ASO 


¥ j 
> 
Tak f 
ie 
aA 
4 
S 


187 


FABRIC 35 (continued) 


Pebble Axial lengths Axial ratios a axis ab plane 
number eee lee | Ce. Sib. BG trend pl str dip 
O7 | O30 .023 OF Leo es Ibo Ot 
O72 OSs ODS OO lefties O25 55 
O73 020 O18 O10 fal [.8 OZae 21 
O74 O28 - OFS - 008 Pio: 250 S46 Wd 
O75 @25 020° . O15 bem [as 238 76 
076 020 O20 Os PRD Wad 216.» @5 
O77 Oe OD OO bea. Tad Gi@ SF 
078 O64 046. 025 aes bee O7@ 2 .D7 
O79 0435 050 030 hea: VO 24 >. 33 
O80 O20 OFS --006 roo ohoo 22 32 
O08 | O25 OFS OFS 252 - £20 O55 52 
082 O45 030 - 020 fae aD 028. 2l 
083 020 O68 Olle in’ 2 lad ONZ - AS 
084 030 025-023 Pe ee 008 37 
085 025 020°. O13 eget Mell We. aie . OS 
086 Zoo. Za 
O87 250} O05 
088 55a) 730 
O89 196 36 
090 Ore? - IZ 
09 | Q245. “(5 
092 185 O09 
093 200) . O6 
094 [735 05 
095 OO} + Zl 
096 2278 28 
O97 OSS) - 44 
098 27 VB 
O99 2a, be 
100 199 22 
10 | 226) - :\o7 
102 2bee is] 
103 2k Q2 
104 O27 28 
105 256,-- 115 


106 285 54 


mn TAC Oe Aw oy Ber Wht an _ | 


vel ' 


tain’ 
Ci aa 


188 


FABRIG 36 (Dimensions in millimeters) 

Pebble Axial lengths AxXtal. rattles a_axis ab plane 
number =e~, Be =e aib bic trend p!| sir dip 
OO | O59, OST O25 ee ae: HFS Of [29° 43 NE 
002 O65 ©5838. 047 eae Lez OFZ 20 

003 O20, CTS OFS pee) he? O22 27 

004 O29... O25" OVS Pe ale ee) OT4 35 

005 ney .: Ors OO Piss Baer: Or6 68 

006 069 O6! 0236 bon! ei s7 OFF $25 (SY .40° NE 
007 (59 « TO6" O85 pea hie 214-01 

008 O29. OTS 2OTZ WO) hes rye 02 [10 ~50: NE 
009 C25 1020 Ol2 Lato? eee OT “he 203 82 NW 
O10 2a Orton Ore ea Gi SAS Tras 

Ol | O27 “925, OTS areal fa4 ri? 53 

Ol2 O57 7025 OFS Perera O87 4} 

O13 045° 050% 025 erik he 2 049 4 

O14 035 ,029 . 024 ave ee [36° 768 

O15 040 024 O18 elt Anas, gyab gt 

O16 O35 20245. 025 ee a ted @) Ol Zero) 

O17 O25 Ol] 7 O45 so aN eS) 048 756 

O18 O55. O21 -.006 (SAG, Bears 0835" 14 (OF 7245Ne 
Ol9 022. DIF ~ O12 aise mat ba: | 54 | ~A5 

020 O64 .050° 2024 oO) a lee wyoyeh 7h 

O21 255° |S7e-096 er eee! PHYO) XO) 

022 O0-56 HOSA ONG leo eee | Peer OS) 192 2| NW 
O23 O25. “OT6" -Ol2 ee Lise LO7: "20 

024 OSO7 FO 27s -O.n6 bail bate LES: «407 NE 
O25 038 O21! . O20 she sa D458) 26 

026 OS 0257 O16 fe ES [O97 10 [15 39 NE 
O27 O22 O18 009 Levee mene <i) O57 + 28 (27) 55> NE 
028 O54 024 020 Aer 2 OTF! 

029 O66 050° 020 202 @25—" 24 217 59 NW 
030 Ore O15" 004 eZee 053 39 267 43 NW 
O3| O20 O19 071 er mapa! 555 427 

O32 O20 O20 008 ae Gay SE OOG +697 'SE 
O55 O29. (O12, 006 eS Vieye a Nees) 349 O08 

034 O538.7.025'- 022 po. + wey S457 27 

O35 048 037.013 POR ee 213 14 205. 253 NW 
036 048 029 022 sy ere VBA 32 | 18 

037 052° 040 027. Peed bao) $56 126 

038 OF Ole Or laD:  WeO 205: 03 

039 055 -O(S O14 [cle te |e 349 O08 

040 O28. .OT 1 —005 2 cae 0094.27 [i729 NE 


aM Ob Cet 


3A. OF | O11 
Wace 0s 


av AS 


IM Ve 
Bi Be 


fanart 
“+ ’ 


I UD 
— — a Cl 


~*~ 


2 


Vor 


+ 


am >- Seem 


Mano 


<< 
a il 


189 


FABRIC 36 (continued) 
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APPENDIX D: 


COMPUTER PROGRAMS 
Unmodified Point Density Program 


This program plots the density of axes or the poles to planes on 
an equal-area projection. The area of the counting circle is specified 
as P, a percentage of the area of the projection. If is the angular 


radius of the counting circle, it can be shown that 


cosa = 1 = P ° 
TOO 


The cosine of the angle ¢g between a grid point with direction cosines 
(14, mj, 1), and an observation with direction cosines (1,m,n) is 
given by 
cos o = 1.1, + mom, + nony - 

Direction cosines (1,m,n) of an axis specified by its trend (TR) and 
hans (PL) are given by 

le= cos(TR) cos( PL) 

m = sin(TR)cos(PL) 

ns San PL). 
During the counting procedure at each grid point, cos g is calculated 
for each observation and compared with cosQA. If cos g2 cosa (that 
is, 6 <Q ) the point counter is incremented. The density values ob- 
tained are plotted as a grid of 333 numbers which can then be contoured 


by hand. 
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Modified Point Density Program 


This modified procedure was described in some detail in Chapter 2. 
The approximation of the probability density surface for a true axis 
employs an expression derived by Watson and Irving (1957, p.293). They 
showed that, for large K, the probability Pr between R, and Ry from the 


centre of the distribution is given by 


Boone e-K(1-cosR1) _e7K(1-cosR2) ; 


Using the specified value of K, Pr is calculated for successive radial 
increments of one degree. The average probability density (that is, 
probability per unit area) within each radial increment is then calcu- 
lated. The areas are determined from the expression 

Aus ez a“ (1-cosR) 
where A is the area of a circle of radius R on the surface of a sphere 
of radius a. Calculations involving spherical probability distributions 
are normally referred to a sphere of unit radius, so that a=l. 

Successive values of probability density are subtracted to obtain 
the thicknesses of the disks used to approximate the probability distri- 
bution (see Chapter 2). 

Using the radius of the counting circle calculated for the 
Specified percent area, P, the volume of the probability distribution 
lying within the counting circle is calculated for centre-to-centre 
distances of zero, one degree, two degrees, and so on until the largest 
disk no longer overlaps the counting circle. These values of volume 
(that is, probability) are stored for use during the counting procedure. 


The counting procedure is similar to that used by the unmodified 
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program, except that the angle g, instead of being compared with the 
radius of the counting circle, is rounded off to the nearest degree and 
referred to the table of probabilities (volumes), and the appropriate 


probability used to increment the counter. 


Common Area of Intersecting Circles on a Sphere 


The calculation of volume referred to above involves, for a given 
centre-to-centre distance, computation of the common area of the count- 
ing circle and each disk of the approximated distribution. An expression 
for this area was derived as follows. Let the diagram below represent 
two circles of radii Ry and Ro with centre-to-centre distance D on the 


surface of a sphere of unit radius. 


The angles A, B, and C are given by identities of spherical trigono- 


MeUrys 


cos A = COSR2 ~ cosRicosD 


sinRysinD 

cos B = COSR1 - cosR2cosD 
sinRgsinDd 

ees ie cosD - cosRjcoskR2 


sinR,sinRo 


The area of the triangle OPQ is then A+B+C-n where A, B and C are in 
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radians. The area of circle 1 is 2n(1-cosRj) and that of the segment 
with angle A is A(l-cosR;). Similarly, the area of the segment of 
circle 2 defined by the angle B is B(1l-cosRy),. The sum of the areas 
of these two segments exceeds the area of the triangle OPQ by the area 
of the shaded portion in the diagram, which is half the common area of 
thescurctes .asLhus 

‘(common area) = (area of segment A) + (area of segment B) 


- (area of triangle OPQ) 


iN 


A(1l-cosR}) + B(1-cosR>) - (A+B+C-1) 
= Spa mee AcosRy - BcosR9 


Hence 


il 


common area 2(m- C - AcosR; - BcosR») « 
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Mean Vector Program 


This program takes a set of N axes or poles to planes and, using 
a specified reference line and angular distance, selects those lying 
within the specified angle of the reference line. Using these n obser- 
vations, it calculates the direction and length (R) of the resultant, 
and estimate (k) of the precision parameter of the population, the 
95% and 99% confidence radii on the direction of the mean, and the 
vector magnitude (R/n), using the formulas given by Watson and Irving 
(1957). 

After the mean direction has been calculated, all the data are in 
effect rotated until the mean is vertical. For each observation the 
Computer then calculates (1) the complement of the new plunge and (2) 
the new trend referred to the original trend of the mean as zero azimuth. 
The trends are then grouped into 72 five-degree classes, and the plunge 
complements into 12 concentric or polar classes each equal to one 
twelfth of the radius of the limiting cone. Using the calculated pre- 
cision parameter k, the computer then derives the expected frequencies 
in the polar classes assuming a Fisher distribution. Output from this 
part of the program consists of the polar class limits, expected fre- 
quencies in polar classes, observed frequencies in polar classes, and 
the observed frequencies in azimuthal classes. 

The next stage in the program calculates chi-square values for 
the azimuth distribution, employing a procedure designed to overcome 
the problem posed by elongate and double groupings. Such groupings, 
while having non-random azimuthal distributions about their central 


direction, will give nearly equal frequencies in four quadrants about 
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a central point when the quadrants are in a certain position. The 
procedure is as follows. Using the beginning of the first azimuth class 
as the origin, the azimuth frequencies are grouped into eight 45-degree 
sectors, then into four pairs of opposite 45-degree sectors, then into 
four quadrants beginning at the origin, and finally into four quadrants 
beginning 45 degrees from the origin. For each of these groupings a 
value of chi-square is calculated, comparing the observed Eregtencie: 
with a uniform distribution. The origin is then moved to the beginning 
of the second ciass and the procedure is repeated. This process is 
continued until all nine positions of the quadrant have been used. The 
chi-square values are printed in the form of a table. 

Using the computer output, chi-square values were calculated for 
the polar distribution after regrouping the classes so that the expected 
frequency in each class exceeded five. The regrouping usually employed 
combined the three innermost classes into a single group, the three 
outermost classes into a second group, and adjacent pairs of the six 
remaining classes into three more groups. No less than five groups 
could be used, since the number of degrees of freedom of chi-square is 
one less than the number of groups diminished by the number of parameters 
of the model distribution that are estimated (Fraser, 1958, p.270). 
Since in this case three parameters are estimated (k and two coordinates 
of the mean) the number of degrees of freedom is C-1-3 or C-4 where C 
is the number of classes. Since the number of degrees of freedom cannot 
be less than one, C cannot be less than five for this comparison. In 
the case of the azimuth distribution, the number of estimated parameters 
1S two, since two parameters are required to define the position of the 


mean and k is not relevant. Hence in this case the minimum number of 


Classesias leltZze4. 
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COMPUTER PROGRAMS 


GC PROGRAM © -...POENT. DENSITY (PRB RL 


54 
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14 
15 


DIMENSTON COSA(333),COSB(333),C0SG{(333),ND{(333) 
DIMENSTGON A(1000), B(LO000), GI{1000) 
O) a Rive Scitiae Hew Dus A ie 
READ(5,500) COSA(T),COSB( 2) ,COSG(I},J 
ORMA EE p 7X 2, 
[Pid eNEsiie GUAaw cs 
CONTINUE 


OREAD (5. 5013 KRE Ay NMP) 0 8 eRe ia) <0, ees A areas ee 


FORMAT(1I354XyI11) 
IF(NM.NE.1) GO TO 24 
PCA=KPC 
TEST=1.-PCA/100. 
READ(5,502)NSTN,gNZ,NN 


IF(NN.NE»«2)G0 TO 24 
NCDS=0 
READ(5,503)N1,N2,NNyNDIR 
FORMAT(I4,13,11,11) 
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IF ONL L720) Ge Tae 
NCDS=NCDS#4]1 


cai ecenmenw ite anasto oierets 


eeceentiimeseartieaerehtrattennce renee trees r. 


IF(NDIRWEQ.0)G0 TO 26 i 
Ps SOs Ne ce seo tn eight bale tee eS ee 


IF(NDIR.NE~1)GO TO 8 
TR=N1+90 

0. TH 

TR=N1-90 

GO TO 11 

PL=90-N2 

TR=N1+180 
IF(TR.»LT.360.)G0 TO 1 
TR=TR-360 

eo° TO aa 

TRENI 

PL=N2 


fs 


P1=3.14159265*P1L/180. 


TR=3.14159265*7TR/180. 
A = ERS RRC OHS{PI 
BIT J=SIN(TR) *¥COS(PL) 
GI TJ=SIN(PL) 
60 10.6... 
CONTINUE 
CDS=NCOS 
IC 5 = 3 
DENS=0. 
CSA=COSA(I) 
hss ok BB NS ws 8 en 
CSG=COSG({(1) 
DO 14 J=1,NCDS 

=ABS(CSAXA +#CSRB*R 4CS6G4G 
IFC TEST.GT.~AVIGO TO 14 
DENS=DENS+1 


iS 
se 


CONTI AUE 2 scien ig a at ce OE i 


ND (1 )=DENS*100./CDS 
DO 60) K=2,'2 
Wi R J EF{6.4609 A NSTIN.NCDS 


CENTRAL STORES HF -585 


600 FORMAT(*1PERCENTAGE OF POINTS PER',I2,1 
LUMBER! ,1492X,y 
PeS5Xsl4>? OBSERVATIONS") 
GO TO (27,28 329),NZ 

27 WRITE(6,606) 


PERCENT AREA 


606 IRMAT(* LINEARS SPECIFIED BY TREND AND PLUNGE ® 


Sou Ti 2s 
28 WRITE(6,607) 


STATION 


N 


ee GOd Fame (Pores DF PUANES: PLANES SPECIFIED BY STRIKE AND DIP!) 


G60 10 25 
29 WRITE(6,608) 


608 FORMAT({' POLES OF PLANES: PLANES SPECIFIED BY DIP DIRECTION AND 


CE Aa 
25 WRITE(6,;601) %(ND{(1), 1=1,+333) 


f 


DOSS) EGRMAT(I50%. 15./39%115.165110.16//30K%, 91 5/7/20, VAISS // LOX LS ITSL Ss) 
POM, 1 799/77 10%5 1 715/775K21915/77/5%51915/7/5X%,1915/716,93X%,15//5X1 G1 
Poy Gs 95 et oF 72 L1D7 716393 %51D/5X%,1915//16293X,1575X,191S///5X,19 157 
SSIEKRe LTS POXs i PISS 7 1LOXS 1 TI5/7 715X531 515 20X%s1315 


49X%,15,16,110,16,/50X%,15) 
60 CONTINUE 
(ERS ee Tae 


eaiwehtetsse02) 


602 FORMAT(* ERRDR....2PQINT COUNTER OUT OF 


SEQUENCES / /* 


LLUMNS 58-60 MUST BE IN SEQUENCE FROM 1-70 333") 


Sauer 
24 WRITE(6,503) 


D60s FIRMA (SY PARAMETER CARD MISSING OR MISPLACEDS) 


S?OP 
26 WRITE(6,605}) 


605 FORMAT(50X.'WRONG DATA FORM CODE'.//59x 


1 MISSING ON STRIKE AND DIP_ DATA CARD) 
100. STGP 


Pete OM foo VOR 


onal oy 


DIMENSTON A{1L0900 Bf 
DIMENSION NOAF (72), NO} 
PT=3.14159 
eee Pe LALS SPOS TSO 


[le thes ie SLT ESO ES a ne PC 


BPEUT2 §, CSC) 


ea READ PARAMETER CARD 


2 READ(5, 502)NSTN,NN,NG 
502 FORMAT(14,2x%,12,513 
IFINSTN.LT.03G0 TO 100 
IF{NN.NE.2)G0 TO 24 


a) 


"READ (5,501) NN, NTRRaNPLRANRAD 
501 FORMATI7X,11,14, 13213) 
IF (NNsNE.3)G0 TO 24 
P=CEXNRAD/12 


oOo 
-) 
i 
ae 
Se) 
aes 
~ 
pd 
—a 
m 


Si boat ALS 
PLR=NPLR 
TRR=TRR*CF 

P} R=Pl R*XCE 
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BREAD PARAMETER CARD. 


DIRECTION COSINES OF 


/1130X.915 


NUMBERS IN 


TION 33 READS PARAMETERS AND 
PRELIMINARY CALCULATIONS 


REFERENCE 


p33 


Lit 


AREF=COS{TRR)*COS(PLR) 
BREF=SIN{(TRR) *COS(PLR) 
GREF=SIN(PLR) 
C CALCULATE COSINE OF SEMI-ANGLE OF CONE 
RAD=NRAD 
RAD=RAD* 
TEST=COS{RAD) 
C INITIALIZE SUMS AND COUNTERS 
SA=0. 
SB=0. 
SG=0. 
NTOT=0 
N=0 
oo 30" tei; 
NOPF(I)=0 ee ore 
3C CONTINUE 


DD. 33) tbs FZ 
7 NOAF(1)=0 
33 CONTINUE 


PROGRAM SECTION 6: READS DATA, COUNTS TOTAL 
__NUMBER OF OBSERVATIONS IN SAMPLE, NUMBER 


WITHIN CONE, AND COMPUTES SUMS OF DIRECTION 
COSINES GF OBSERVATIONS WITHIN CONE 


tHe ey 


502)STR,DIP,NDIR 

503 FORMAT(F4.0,F3.0;12) 
IF(STR«LT.0.)GO TO 14 
NIOT=NTOT+1. 


Gaui, CALCULATE DIRECTION COSINES OF OBSERVATION — 
C VECTOR (AXIS OR NORMAL TO PLANE) 


NDTR=ENDIR+ 
GO. 10) (103545) -NDIR 
5 pi=90,-DIP 
TEUNDIR «NE 2G TO a ee eee 
TR=STR+90. 
GO TG 14 
8 TR=STR-96 
60 Tavis 
10. FR=STR 
PLSD I Paes 
11 PL=PL*CF 
TR=TR*CF 
={Cs RY*COS (TP 
BB=SIN{(TR)*COS{PL) 
GG=SIN(PL) 


__ CALCULATE COSINE OF ANGLE BETWEEN REFERENCE 
LINE AND OBSERVATION VECTOR 
13 CALPH=AREF *AA+ BREF*RB+GREF*GG 
F OBSERVATION ES. JE CONE Gf 0 


a te 


G NEXT OBSERVATION 

ITF {ABS{CALPH).LT-TESTIGO TO 6 

DESI a a ee er 
REVERSE OBSERVATION VECTOR IF NECESSARY 

IF(CALPH.GE .0.}60 TO i2 

AA=-AA 

88=—-s5 

GG=-GG 


cs 


__... INCREMENT SUMS OF DIRECTION COSINES AND RET 
TO READ NEXT DATA CARD 


Cie 


12 SA=SA+AA 
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C 


SG=SG6+6G 
A(N)=AA 
BI N)=BB 
G(N)=GG 
GO TO 6 


PROGRAM 
SvaTISTICS 


SE 


Oh 4s CULATES AND PRIN 


14 WRITE(6,610) NSTNsNG, NTOTSNTRR yNPLRNRAD Ny SAySBySG 
— S10 FORMAUULIT/ (1/213 39X% 2 SSTATION NUMBER '+145'-'s12.5Xs "TOTAL SAMPLE: 4s 


1 
2 
a 
4 


ee 
[F(SG.GE.0.)G0 TO 15 


ee aie 
10X,tAXIS 
DEGREES 


OBSEI 


OBSERVATIONS! // 
OF CONE AT*,14,13;5X,"SEMI-ANGLE OF 
PATIONS WITHIN 


CONE:',13, 
ONE: "1 


20K, *SA=",F11.4,5K ,*SB=4,F11.455X, 'SG=',F11.4//) 
18 R=SORTISA*SA+SR*SB+SG*SG)} 


AN=N 


Vo 


f 
$B=-SB 
SG=-SG 


eM Ae CRC ASS 2 


16 


VPL=ARSIN(SG/R) 


IF (VIR.GE.0.)G0 TO 16 


ie R4+180 
IF{S8.GE.@.)G0 
VTR=VTR+180. 


re wee de 
Q5=20.**(EXP)-1 


QQ5=1.-Q5*{ AN-R 
AL PH5S=ARCL Q@5 
Q1=100.**{ EXP) —- 
Q91=1.-Q1%* { AN-R 


i Se 


TO 


Y/R 
CF 

1. 

)7R 


i7 


BNA at ANUS LAN. 6s ee 


600 
2 


320X%—9'95 PERCENT CONFIDENCE RADIUS ABOUT 
a -AZ0K%s "99 PERCENT CONFIDENCE RADIUS ABOUT MEAN 
ZOKs "ESTIMATED PRECISION 
Z2OX,* MAGNITUDE OF RESULTANT VECTOR 
N)*100 


5 
6 


AAO 


XK=(AN-1.)/ (AN- 


R) 


WRITE(6, 600) VTR»VPL,yALPH5,ALPHI,XK,R,RL 
MEAN VECTOR 
OF MEAN 


FORMAT(20X 
20X5'PLUNGE 


OX» *VECTOR MAG 


ee 


BR=SB/R 

pR=SG/R 
CVPL=SORT(1.—-GR 
SVTR=BR/CVPL 


‘TREND OF 


N 


*GR) 


DE 


i. 


VECTOR 


=(R 


CALCULATE 
AND 


BES ee Av Ne 8 ee ee 


BG Si ta riN 


CPH I=AR*A(T)4+BR*BCI)4+GR*G( 1) 


E(CPH a! 
ACI )=-A(T) 
B(T)=-B{T) 


C 


O° FG 


49 


Be 51 a et ee 


49 


CPHI=-CPHI 
PHT=ARCOS(CPHI ) 
% TR+R 


Re = fh 
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Bra 


TR 


PARAME 


RESULTANT, 
AG VERTICAL PLANE THROUGH 


=',F6. ‘DEGREES! £7 

ee Sag TEBE Gee otf ¢ 

MEAN =*,F5.1,' DEGREES? // 
= 7% st 5 1 ? i DEGREES*// 
TER OF POPULATION =",F7.25// 
(A =", FG6.1// 
PERCENT =" .F5 : 
POLAR ANGLES BETWEEN OBSERVATIONS 
AND AZIMUTH ANGLES REFERRED 
RESULTANT 


PERCENT! / 


CSC=(CBE-CPHI*CVPLI/(SIN( PHI )*GR) 
TE(CSC.LT..0.99991G0 10 S32 
CSC=0.9999 

32 TF(CSC.GT.-0.9999)G0 TO 34 
Coo=-U.6 799 

4 C=ARCOS(CS 
XYZ=BR¥*A(1)-AR*B(1) 
IF (XYZ.LT.0.j)G0 TO ef 
C=2 *PI-C 


- ee ARS GRUNTS a 


27 J=36.*C/PIt1. 
NOAE =NOAF(J)41] 
J=PHI/P+1. 
NOPF(J)=NOPF(J)+1 
3] CONTINUE ee waa nih lean a kee ee 
C “CALCULATE POLAR CLASS LIMITS, EXPECTED 
c FREQUENCY IN POLAR CLASSES 
OX 16 
Da 29 Pods ie 
F=I 
PLIM(L)=PxF/CE 
Q1=—XK +XK % KCOS (PF) 
Q=2.71828**Q1 
EPF {1 )=AN*(QX- 
QX= 
29 CONTINUE 
C ___ PRINT CLASS LIMITS AND OBSERVED AND EXPEC TET 
Cc FRFQUENCIES 
WRITE(6;648) (PLIM(I), T=1,12) 
8 FORMAT(* POl AR ASS MIT V1 ZFS 
WRITE(63647) (EPF(T), I=1,12) 
647 FORMAT(!' EXPECTED POLAR FREQ 7 12F8.0/) 
_-WRITE(6:646)  (NDPE(T)s Ielsl2) PRED EUE MEI NA he Bienes Bc RC 
646 FORMAT(* OBSERVED POLAR FREQ 1 L2neys2 i, 
WRITE(6,645) (NOAFCI), T=1e 72) 
45 FORMAT{* OBSERVED AZIMUTH FREQ *,1815 
a CALCULATE CHI-SQUARE ON AZIMUTH DISTRIBUTION 
C USING DIFFERENT PATTERNS OF CLASSES AND 
_NINE DIFFERENT ORIGINS AT 5-DEG INI ERVALS 


™ + 
oe 


WRITE(6,700) 
700 FORMAT(///10X%, * CHI-SQUARE VALUES FOR AZIMUTH DISTRIBUTION: '// 
22X.'180-DEGREE DISTRIBUTION’ Xs 1! Xs '360- REE DISTRIB! 
2ON'/19X%, #2 CLASSES’, 19X%9'1's5X9'4 CLASSES", 5X_"E%39Xs "4 CLASSES 
aX GLa 8 es se wren bee 45-D SECT',3Xe21'1%,9X, "QUADRANTS *%, 
T<0 
36 S1=0 
Sa=Q 


D0 37 1=1;9 
37 $1=S1#NOAr11) 
PE BeBe i ee a acc a 
38 §$2=S2+NOAF (1) 
DO 39 1=19; 2t 
Qq =S24NDA 
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DOD 40 1=28,36 
40 S4=S4+NOAF(T) 
DO 41 1=37,45 
41 S5=S5+NOAF (1) 
DON 42 1=46,54 


6= DOFNDA 


OB 43 1=55,63 
43 ST=ST+NOAF (1) 
Ye Mit fo Leak AS a 
44 S8=SB+NCAF{(1) 
46 IFIN.LT.10)G0 TO 48 


ay 


CS(1)=(1S14+S3+S5+S7-E)**24+(S24+544564+S58-F)**2)/5F 
CS(2)=((S14+S2+S54+S6-E) **24+( S34S5445 74+S58-E eK?) /E 
a esi BI SIUS24534S64+S7-F)**24+(S8+S1 +S44S55-F ) **2) /E. 


ICS$=3 


IF(N.LT.20)G60 TO 45 


==15 


CS(4)=((S1#S5-E )**24+(S2+S6—-E )**24(S 345 7-E) ®K24(5445B-F) RQ) /F 
CS(5)=((S14S2-E )** 2+(S34+S4-F ) **24( $54S56-F ) **24( 5745 8-E ) kk?) /E 
Be AO VENTS Ae SA SE) 2 +(S4FESS—E ) 2 41 S64S7—-EF) **24+(S84+51-F) *e2)/E 


ICS=6 


IF(N.LT.40)G0 TO 45 


2 


=H 


CS(7)=({ S1-E)** 24+(S2-F) e241 S3-E) R24 ( S4-E) KK24 (S5-F) KK24(S6—F) RR 


1 41S 7-£) *#24+(S8-E)**2)/7E 


Se 5. Sa 


45 WRITE(6,643) (CS(I)>I=1, ICS) 
643 PORMAT(1X+F11.2,7Xs20*T*,F12.2)F 8.2, 7X_ 1% 9F120257X)) 


=[T+ 


IF({7.GE.9)GO TO 3 
NOAFI1=NO AF (1) 
CU ee | ye Se EU a ee 
47 NOAF(IJ=NOAF(I+1) 
NOAF(72)=NOAF1 


a TO 3 


48 WRITE(6;,641) 
641 FORMAT(10X,"SAMPLE TOD SMALL TO CALCULATE CHI-SQUARE') 
SWS Nee 2 <r 


: 
oO} 
' 


ERROR MESSAGE 


24 WRITE(6,603) 


603 FORMAT 
100 STOP 
END 


C PROGRAM - 


PARAMETER CARD MISSING OR MISPLACED") 


VECTOR MAGNITUDES 1 


_WIMENSTON COSA1334).C00SB(334).C9SG6(334) NL 1334) 
DiIMENSTON A{1000), B{1000), G{1060) 
c READ AND STORE DIRECTION COSINES GF NORMALS TO REFERENCE PLANES, 


Pe eA i 


iA ARDS RE Nj FQUENCE 


PT=3.1415926 
DO 2 Y=1,333 


pee tos OOP COS AIT) -COSBIRVACOSG(T)s tS 
2 FORMAT (7X,3F11.8,17X_,13) 


IFCI.NE.J)GO TO 23 


q 


IN 
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oo) 


17 PL=P 


ca 


aa 


fon) 


14 


4h 


A 


READ PARAMETER CARD 
READ(5,502)NSTNgNZsNNGNG 
FORMAT(14,12;12,13) 
IF(NSTNeLT.G)IGO TO 100 
TF (NN.NE.2) TO 24 

=N 


READ 


{% 
GO 


DATA: COUNT CALCULATE DIRECTION COSIN 

ES AND STORE 

NCDS=C 

READ(5,503)STR,DIP,NDIR 

FORMAT(F4.0yF 340912) 
-{S1R (Fs GS ee 

NCDS=NCDS#1 

GO TO (9,10) ;NZ 

TF(NDIR~EQ,0)60 

PL=90.-DIP 

IF(NDTIR.NE-1)GO TO 
R=STR+90 


“fe oh 


OBSERVATIONS,» 


2 £ 


TO epee chaskomeks 


a2 


Nie 
GO 


TR=STR-90. 


GO 
TR=STR 
PL=DIP 


10.12 


TR=H=PI*xTR/I8C0. 
l=NCDS 


ALTI=ECOSLIRIFCOS (PL) ar ce 


B(T)=SIN(TR) *COS (PL) 

GCL)=SIN(PL) 

GO TO 4 

CDS=NCDS 

LOOP TO CONSIDER EACH REFERENCE PLANE 
DO 6 1=12333 


CSA=COSAIT3 

GCSB=COSB(LT9 

CSG=COSG 

INITIALIZE SUMS OF DIRECTION COSINES 

SA=0. 

SB=0. AN 


SG=0. 
LOOP TO CONSIDER & 
Quiz. = PNG 

CALCULATE COSINE 
OBSERVATION 


ACH OBSERVATION 


OF ANGLE BETWEEN NORMAL TO REFERENCE PLANE AND _ 
VECTOR 


CALPH=CSAXALII4CSB*B{ JI +CS OCU IP orice 


SIDE OF REFERENCE PLANE 


TO 


ADJUST OBSERVATION 
IF{CALPH.GE.0.)G0 
SA=SA-A 

SB=SRB-B(J) 
$G6=SG6-GtJ) 


ee 


VECTOR 
TO 8 


LOWER 


INCREMENT SUMS OF DIRECTION COSINES 

SA=SAtACJ) 

R= SBt+B 

SG=S6+G(J) 

CONTINUE 

CALCULATE MAGNITUDE GE RESUI TANT AND EXPRESS AS PERCENT 
R=SORT (SA**24+ SBR**24+SG%*%2 ) 

NL{I)=100.*R/CDS + 0.5 

ANTIN 
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sl PRINT RESULTS AND RETURN FOR NEXT DATA SET 
WRITE(6,6090)NSTNyNG,NCDS 
600 FORMAT{*1VECTOR MAGNITUDES (100*R/N) FOR 3233 POSITIONS OF REFERENC 


1— PLANE STATION NUMBER? ,14,°—',12+10X%,14," GBSERVATIONS',/* Ri 
2GHT) HAND DIGIT GS At PGWLE OF REFERENCE PLANE') 
WR 6260 {N =} 4333 


601 FORMAT(50X,15,/39X%.15,165110,16//30%,915///20X,1315/// 15X15 I5///]) 
Doty SI LONG ae Ose OAS ADIT SSF POX ROIS ///5Xo 1915 /16,93X%,154/5X,191 
25116,93X%s 1 5//21 15/1659 3K 2 LS/SXsILIIS//16 293K TOSS LOLSS//OX SLOT SS 
OUP Seg hol o/ VLOG TITS /7/10X%4 1 1157/7 15K,15157//20X%,1315///30X%s9T5 7/3 
49%y15,16,110,16,/50X,15) 

GO T0 3 
C ERROR MESSAGES 
23 WRITE{16;,602) 
B02 FORMATI* ERROR. ... COUNTER CARDS OUT DE SEQUENCE /' NUMBERS IN CO 
I2IMNS So-bO MUST BE IN SEQUENCE FROM 1 TC 3337) 


STOP 
4 WRITE(6. 
603 FORMAT{* PARAMETER CARD MISSING OR MISPLACED - SEE INSTRUCTIONS) 
STOP 


I ae ee le : 
665 FORMAT(50X,'" WRONG DATA FORM CONE’, //59X%,!OR',// 35X,'DIRECTION CODE 
1 MISSING ON STRIKE AND DIP DATA CARD?) 


00 STOP 
END 
MODLETED POINT DENSITY .4 MF PD4 


aaA 


PART 1. -— DIMENSION VARTABLES AND DEFINE. CONSTANTS. 
DIMENSION €0SA{333), COSB{(333), COSG(334), ND(334) 
DIMENSIGN A{1000), B(10090), G(1000), VOL{(100 
DIMENSION R{(45), DPRIN{45), AR(45) 

P1=3.1415926 
E=2.71828 

BP alsSe hie |.) 1 EL eee eee eae 
HD=9 .O1*CF 


C 
a 
3 
C PART 2 - 
_C READ AND STORE DIRECTION COSINES OF COUNTING LOCATIONS » CHECKING THAT 
C COUNTER CARDS ARE IN SEQUENCE 
pO 211=1,333 
READ(5,500)COSA(T) »COSB(T) »COSGI] 
500 FORMAT(7X,3F 11-81 7X 413) 
IF(I.NE.J)GO TO 23 
Bae ON Ee. 
C 
C 
C PART 3  -— READ PARAMETER CARD 1, PERFORM CHECKS, AND CHANGE 
C PARAMETERS TO FLOATING POINT 


_ READ(5,501) KPCAsNDISPs»NN2KP, KORY 
501 FORMAT(13,12,13,14,12) 
IFINN.NE~1)GO TO 24 
PCA=KPCA 
CENTRAL STORES HF-585 


CALCULATION OF PROBABILITY DISTRIBUTION 


SDPR=0. 


EXP1L=-—XK *QX 
EXP2=-XK*Q 
DPR=E**EXP]—ExXxX EXP 2 
SPPR=SDPR+DPR 
AR{T)=2.*P TO 
ART=AR UT} 
fFLECEO. 1L).GO 16 2¥ 
DA=AR CTI ~ARLT=—1) 
PRU=DPR/DA 
DPRU(I-1)=PRUX—PRU 
[FINDISP.EQ.0)60 TO 41 
5,610) 1, DPR, SDPR» ARI » DAs PRU 
610 FORMAT(5%,?1=',13,5X%,'DPR="_yF6.4,10K, SDPR=*%,F6.4;s 
} 10X%,* A=", E10. 3710X%_ *DA=",E10 -35910X, *PRU=',E10.3) 
41 JEL SDPRG1.0,.9993760 [0 tt 
GO FG 22 sei } a 
27 PRU=DPR/ AR (1) 
FC ND P,FD.0) C0 ieee 
WRITE(6,612)1,DPR,SDNPR,ZARI, PRU 
612 FORMAT(5X,_,*I=" 513, 5X_*DPR=" ,F6.4e10KX,* SDPR=",FO.4s 
1 10X%s* A=". F10.3233%s 1 PRUH" FIG e ee eee 
22 QX=9 
PRUX=PRU 
Tees ENE DE fa 
17 DPRU{T)=PRU 
[MAX=I 
TEIND ISP ER CO Ti ae 
WRITE(6,609) {(DPRU(T) y i 
609 FORMAT(/10X%,"*DPRU IN ORDER O 


=1, 1MAX) 
F INCREASING 1'/10X,10F10.4) 


SR 


PART 5 -— CALCINATION GE AREA AND RADIUS OF COUNTING Cigt. ee 
42 ACC=2.*PI*PCA/100. 

CSC=1.-PCA/100. 

=fAPo ns See 

CD=C/CF 

ICD=CD+0.5 

LE{NDISP.£90.0)1G0 TO 43. 

WRITE(6,608) CD 
608 FORMAT(/10X,"*RADIUS OF COUNTING CIRCLE=",F5.2///) 


oe ae Se ae Oe 


PART A. COMPUTATION OF TRE UARR Ar eri eee eee ee ee eee 
NERO=90 
43 1D02=I1CD+I MAX 
= 1[CO-IMAX 
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IF(IDDL.GF.~0)G0 TO 44 
TDD1=0 
GO TO 45 
44 00 46 J=1,10N01 
VOL(J)=1. 
46 CONTINUE 
45 DO 47 J=1DD1,1DD2 
D=J*CF 
CSD=COS(D) 
VOLUM=0. 
[=I MAX 
F=] 
48 RA=R(1) 
Tt. WESABSTC-RAI+HDIGG TC 50 
IFAD. GT «CtRA-HD)GO TO 49 
CSR=COS(RA) 
CAL=(CSR-CSC*CSDI/(SIN(UC}*SIN(D)) 
CBE=(CSC-CSR¥CSDISISIN(RA)*SINID 
CGA=(CSD-CSC#CSRI/(SIN(C)*SINIRA)) 
LF CABS( CAL V2GT.1.)60 TO 36 
LELABS(CBE) .G6751,)60 TO 36 
TPUABS(CGA).GT.1.)G0 TO 36 
33 ALPHA=ARCOS(CAL) 
BETA=ARCOS(CBE 
GAMMA=ARCOS{(CGA) 
AREA=2 .*{P1-GAMMA-—ALPHA*CSC-BETA*CSR ) 
DLS ROG 3. Sie Ea 


S SBR ITERG) 606)CD+JS, 1MAX, * »CAL »CBE,CGA 
606 FORMAT (10X, * #* eR KRS MEME ARK ai a Deak i as aie HE ak OK ack tf 


Rar ek KKK TNVAL TID COSINE 
Poke RADTUS 0 DUNT IA 2IRCLE == ' B46? ‘ 
3 bone (UP OIANM.© IN SIENTHS OF DEGREE =*,705/ 
4 15Xs*IMAX =*,14,20XK,'1=',14/ 
5 2 Oks CAL = "615.551 0X, "CBE=" ELS eas LOX, UGAS9 ELS eb L) 
AREA=0. 
NERO=1 
G0 T1052 
SO UIFLAR(I 7 2GE.ACTCIGO TO 51 
ARFA=AR{T) 
Mi Aloe 228 BNE ele a 
51 AREA=ACC 
52 VOLUM=VOLUM+AREA*DPRU(T) 
TF(1.GT.0)G0 TO 48 
49 MOA Sl p=V OL UM 
47 CONTINUE sd 
IF (NERO.EQ. 1)G0_ tO 100 
IDD3=1DD2+1 
FINDISP.EQ.0)G0 TO 3 
DO 55 I=1,1DD3 
IDTSt=1-1 
Mowe Pee Ol ONO S Te VOLT Bec tod. 
613 FORMAT (5X, * DISTANCE =',15," DEGREES',20X,"VOLUME =',F6.4) 
55 CONTINUE 


~~ 


ed 
re 
DEG 


PART 7. = READ PARAMETER CARD 2» READ DATA, CALCULATE DIRECT TON 
~ COSINES AND STORE, COUNT DATA CARDS 
3 READ(5,502) NSTN»NZ.NNING 
5 RMAI . , an 
CENTRAL STORES HF-585 


C 
C 
% 
G 


. 
3 


IF{NSTN.LT.0)GG TO 100 
IF(NN.NE.2)GO0 TO 24 
NCDS=6 


6 READ(5,503)N1yN2—NN,NDIR 


503 


8 


FORMAT(14,13,11,11) 
F(NN.N 0)GD {joe 
IF{N1.LT.OIGO TAG 12 
NCDS=NCDS+#+1 

I=NCDS 

GO TO (95547) 5NZ 
ITF(NDIR.EQ.0)GO TO 26 
P|} ~90-N7 
IF{NDIR.NE.1)GO TO 8 
TR=N14+90 

GO. 10-11 

TR=N1-90 

GO TO 11 

Pj =~9N-N 

TR=N14+180 
IF(TR.LT.360.)G60 TO 11 


_TR=TR-360, 


GO 70 ll 


9 TR=N1 


A PLN 


1} 


12 


PL=PL*CF 

TR=TR*CF 
ALVI=HCOS(TRI*COSTPL) 
BCTJ=SIN(TR)*COS(PL) 
G{T)=SIN(PL) 

G.7n 4 

COS=NCDS 


PART 8 - COUNTING PROCEDURE 


iit oe a prere cnc Pais le oe te le ek 


DDMAX= {) +0 5 
DMAX=DDMAX*CF 


COUNT=0. 
CSA=COSA(L) 

CSB C7) 
CSG=COSGIL) 

DG 54 J=1,NCDS 


SDE ABSTUSSRAI tC SB eB tC SGeG (I) eee 


IF(CSD.LE-CSDM)GO TO 54 
K=ARCOS(CSD)/CF+1.5 
QUNT=COUNT+ VOL (kK 


54 CONTINUE 


Be 


ND (L)=COUNT*100./CDS 


BOS MAC LEL AR SC Rn OE Ne See OT 


C PART 9 ..=—, WRITE RESULTS AND RETURN FOR NEXT OATA. SET 


DO 40 J=1,KOPY 
WRITE Gs: 600) KPCAsNSTN;: NGsNCDSsKP v55%5 i 


600 FORMAT(TIMODIFIED POLE CQUNT PER',12,' PERCENT AREA STATION NU 


IMBER',14,9!'-",I2,19,* OBSERVATIONS, PRECISION PARAMETER (K) =", 1457 


'N BER RUNCATED h ix NDED MEppD4t 


CENTRAL STORES HF-S5&5 


WRITE(6,601) (ND(1I), IT=1,333) 

601 FORMAT(50X,15,/39X%515,16,110,1T6//30Xs 915///20X%,1315///15Xs1515SSS/1 
THX) PITS SOs Ul OST (SMe LOTS /7T 15X91 915/775 X91 915/16,93X%,15//5X,191 
25/16 93% 119//2115//16293X%215/5X%s 191 5//16,93XyT5/5X,1DT5S/S//5X,1915/ 
SF /5Ko7 a1 Ses 1Oks PT IS /SP7 OKs 1 T1S77//15X%.1515/77/20%5,1315///30K%,915/7/3 

| 49 6210.1 65/50X% 515 

40 CONTINUE 

GO TO 3 


Lay he BL a 


PART 10 - ERROR MESSAGES 
23 WRITE(6,602) 
602 FORMAT(' ERROR...+-POINT COUNTER OUT GF SEQUENCE'//* NUMBERS IN CO 
) |  _sULUMNS 58-60 MUST BE IN SEQUENCE FROM 1.79 3334) 
STOP 
24 WRITE(6,603) 
2 FORMAT(! PARAMETER CARD MISSING OR MISPLACED -— SEE INSTRUCTIONS? 
sToP 
26 WRITE(6,605) 
) 1 605 FORMATI50X, "WRONG DATA FORM CODE',//59X%."OR'2//35Xs "DIRECTION CODE 
1 MISSING ON STRIKE AND DIP DATA CARD‘) 
100 STOP 
) END 


) 
) 
) ne fe 


_ CENTRAL STORES HF-585 


Sample Number Original Program Modified Program 


» 3] 6e-7e 
5e-6e 
4e-5e 
3e—4e 
2e-3e 
le-2e 
a o=tc 


; 


—} 


FIGUREMeE 


Comparison of diagrams 
produced by unmodified and 
modified point density programs. 


All diagrams prepared using a 3% 
counting circle. Contour interval = 
e= 3%, the expected density for 
a uniform distribution. N= number 
of axes plotted. North at the 


top in all cases. 


Size of Gyosune ti nieg GCirele 


5% 


FIGURE 4. 


Contour interval = 


j | P Key to Fabric Diagrams 
e, where e is the expected density 


| 
F : : in the case of a sample from a uniform distributj = 
Effect on point density diagrams of See Sthip pop. -  Bye-ize FE 7e-8e §~= (ED 3e -4e 
and is equal to the per cent size of the counting [EB] ice-tle  ce-7e [Ei 2e-3e 
North at the top in all cases. [I 2c-10e [ Se -6e ple 
(HBB) ce- oc [MB 4e-se [_] o-te 


change in size of counting circle. 


circle. 


Key to Fabric Diagrams 


ll ze - Be (| se 
i] ce - Te (2 
isc - 6e Lie 
l-- = SG L_]| 


All diagrams prepared using a 3% counting circle. 
Contour interval =e = 3%, the expected density for 
distribution (505-507) and a_ hypothetical a uniform distribution. N=number of axes plotted. 
girdle distribution (508-515). North at the topin all cases. 


FS WGWits ©, Diagrams of samples from a _ uniform 


| 
| 


Location A Location D Location E Location F 


TPIGUIRE 6. Fabric diagrams, Locations A, D,E and F, except. Key to Stratigraphic Diagrams 


upper till at Location F. Key to Fabric Diagrams Stratigraphic Units General 


ES Lake Edmonton Sediments Ground Surface 
EBB 2c 100 ae ES] upper Till LJ Lower Limit of Exposure 
ED e- 9 sirraals Tofield Sand 
Ese erie: =] Lower Till 
[E) 6e-7e [aa le—2e ed 


the expected density for a uniform distribution. N= number of axes plotted. (HB se-6 ori 
L: moe = lle 


| if Not Ground Surface 


All diagrams prepared using a 3% counting circle. Contour interval = e=3%, 


3 Upper Limit of Exposure 


Diagram Not Extending to 
Base of Exposure 


Saskatchewan Gravels 


Edmonton Formation 


North at the top in all cases. 
Vertical Scale for All Diagrams 


Location M 


FIGURE 7 


Fabric diagrams, 


Locations M, G and H. 


All diagrams prepared using a 3% counting circle. 


Contour interval =e = 3%, the expected density for 


Q uniform distribution. N= number of axes plotted. 


North at the top in all cases: 


Location G 


Key to Fabric Diagrams 


i] 9e-l0e 4e-Se 
[Mh] se - 9e se-4e 
[HBS 7-8 2e-3e 
WBlce-7e le -2e 
[i] se-6e [J O -le 


Location H 


Key to Stratigraphic Diagrams 


Stratigraphic Units General 


f=] Lake Edmonton Sediments Ground Surface 

[9] Upper Tit | al 

| Tofield Sand Upper Limit of Exposure 
[5] Lower Till if Not Ground Surface 
Saskatchewan Gravels 


Es) Edmonton Formation 


| 
Vertical Scale for All Diagrams 1": 20' 


Lower Limit of Exposure 


Diagram Not Extending to 
Base of Exposure 


Location 


FIGURE 9. ~— Fabric diagrams, 


Locations |, J and K. 


All diagrams prepared using a 3% counting 
circle. Contour interval = e= 3%, the 
expected density for a uniform distribution. 


N = number of axes plotted. North at 
the top in all cases. 


Key to Fabric Diagrams 


| 5e -6e 


ae = se 
(MM) Se —4e 
[Mil zee 
] le -2e 
_]| © = ie 


Location J 


Location K 


Key to Stratigraphic Diagrams 


Stratigraphic Units | General 


Ground Surface 


(9) Lake Edmonton Sediments 
= Lower Limit of Exposure 


Upper Limit of Exposure 
if Not Ground Surface 


[7] saskatchewan Gravels Diagram Not Extending to 
[9] Edmonton Formation Base of Exposure 


Vertical Scale for All Diagrams 1’: 20 


WNW ROAD LEVEL: ELEVATION DATUM 


SASKATCHEWAN 


All diagrams prepared using a 3% 
counting circle. Contour interval = 
e= 3%, the expected density for 
a uniform distribution. N= number 
of axes plotted. North at the 


top in all cases. 


Key to Stratigraphic Diagrams 


Stratigraphic Units 
eee Lake Edmonton Sediments 


Upper Till 


fal Tofield Sand 


‘| Lower Till 
Saskatchewan Gravels 
fe) Edmonton Formation 


Vertical Scale for 


General 
Ground Surface 

eee Lower Limit of Exposure 
Upper Limit of Exposure 
if Not Ground Surface 
Diagram Not Extending to 


Base of Exposure 


All Diagrams 


FIGURE 


10. 


Fabric diagrams, 


Location B. 


to Fabric Diagrams 


| 
[4 


lle-|2e 
10e- Ile 
9e -10e 
Be - 9e 
7e -8e 
6e- 7e 
5e- 6e 
4e-Se 


HGOSB00Es 


3e- 4e 


] 
LJ 


2e- 3e 


U 


le -2e 


U 


Key to Fabric 


300° 


FIGURE It. 


Diagrams 
Shanes Mahe Reproducibility of 50-axis samples from sites 34, 
i a37 35 and 37 at Location F. 
Per 


lLiewer Wall 


Lake 


All diagrams prepared using a 3% counting circle. 


| 
Sackaicwemere mGrevels ; Contour interval =e = 3%, the expected density for 


a uniform distribution. N= number of axes plotted. 
North at the top in all cases. 


Composites of 
Horizontal Groups 


East End 


Centre 


West End 


9000 05 


Middle 


Bottom 


39 Lake Edmonton Sediments 


Upper n+ 35 = 
' 


Saskatchewan Gravels 


Location of Sample Sites 


Key to Fabric Diagrams 


Trend and plunge 
024 10 of preferred 


D] 9e-l0e (PS ] 4e-5 
Cae. Sel 7 orientation 


055 o7 Dip direction and 
7 dip of preferred plane 


FIGURE 


Composites of 
Vertical 


Groups 


12 


Fabric variability in upper till at Location F. 


All diagrams prepared 

Contour interval =e= 3 
a uniform distribution. 
North at the top in all 


using Gg 3% counting circle. 
%, the expected density for 
N=number of axes plotted. 
cases. 


< indicates sample site 


several 


43 » 
— hundred 


yerds 
OAD 


466-| 466-8 


Sketch (not to scale) of sample locations showing 


line of stratigraphic section 


Key to Stratigraphic Diagrams 


Stratigraphic Units General FIGURE 13 Fabric diagrams, Location L. 
Key to Fabric Diagrams 


Lake Edmonton Sediments [| Ground Surface 
[ilee-loe ea [J] Upper Till 
[IM] 8e-9e = Tofield Sand Upper Limit of Exposure 


ia 9) Lower Till if Not Ground Surface All diagrams prepared using a 3% counting circle. 
7e-8e 


Lower Limit of Exposure 


[EBB] saskatchewan Gravels Bicapem Wey Extending te Contour interval = e=3%, the expected density for 
Edmonton Formation Base of Exposure a uniform distribution. N=number of axes plotted. 


North at the top in all cases. 


Vertical 


FIGURE 


14. 


Effect of a:b ratio on a-axis orientation. 
Data from sites 34, 35 and 36 at Location F. 


to Fabric Diagrams 


BR ce-7e 
fe 5e—G6e 
pa) 

= Se—4e 
2e—Se 
[4] la =2e 


— 


|) OH le 


po /% contour 


All diagrams prepared using a 3% counting circle. 
Contour interval =e= 3%, the expected density for 
a uniform distribution. N=number of axes plotted. 
North at the top in all cases. 


Effect of a-axis length 
on a-axis orientation in 
two shape categories. 

Data from sites 34, 35 
and 36 at Location F. 


a>25mm 


to Fabric Diagrams 


9e -|0e 


|} Be - 9e 


7e - 8e 


he] 6e- Te 


5e - 6e 


4e-5e 
3e- 4e 
2e- 3e 
le - 2e 
|| ©= te 


Poles af a:b planes 
from sites 34, 35 


and 36 at Location F. 


FIGURE |5 
All diagrams prepared using a 3% 
counting circle. Contour interval = 
e= 3%, the expected density for 
a uniform distribution. N= number 


of axes plotted. North at the 


top in all cases. 


FIGUIRIE 


Key 


to Fabric 

Sz] 6e- 7e 
[__] 5e-6e 
[) 4e-5e 
[0] se -4e 
2e -3e 
le -2e 
[4] ©=te 


16. 


Diagrams 


Sampling experiments using fabric data from 
upper till at Location F,. 

Top row: 50 measurements from 5 sites 

Middle row: 80 measurements from 9 sites 

Bottom row: 200 measurements from 9 sites 


All diagrams prepared using a 3% counting circle. 
Contour interval =e = 3%, the expected density for 
a uniform distribution. N= number of axes plotted. 
North at the top in all cases. 


